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ABSTRACT 


This  report  presents  a new  computer  program  that  can 
be  used  to  design  and  predict  the  performance  of  contrarotating 
propellers.  The  new  program  utilizes  the  latest  numerical 
computation  techniques  developed  for  the  design  of  contra- 
rotating propellers.  The  hydrodynamic  pitch  angle  distribution 
is  specified  as  input  and  the  design  calculations  are  made 
using  Lerbs’  moderately  loaded  single  screw  lifting  line 
propeller  theory.  The  propeller  interaction  effects  (the 
most  Important  new  feature  of  the  design  procedure)  are 
obtained  using  Kerwin's  field  point  velocity  program  developed 
using  finite  bladed  lifting  surface  theory. 

An  analysis  of  results  obtained  for  a sample  set  of 
contrarotating  propellers/  designed  using  the  nevf  method 
presented,  show  that  different  design  and  performance  pre- 
dictions are  obtained  for  the  aft  propeller  when  results 
are  compared  to  those  calculated  using  the  old  method. 
Calculations  made  using  the  new  design  method  are  considered 
more  accurate  due  to  the  improved  method  of  dei'errainlhg 
the  propeller  interaction  effects. 

A FORTRAN  listing  of  the  new  program,  developed  to 
run  on  the  computers  at  the  David  W,  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC)  is  presented  as  well  as 
input  and  output  obtained  for  a sample  set  of  contrarotating 
propeller  designs. 





ADMINISTRATIVE  INFORMATION 


This  work  was  sponsored  by  the  Naval  Sea  Systems  Command « 
SEA  034,  and  carried  out  under  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (DTNSRDC)  Nork  Unit 
1544-256,  Task  14438. 


INTRODUCTION 

The  David  H.  Taylor  Naval  Ship  Research  and  Development 
Center  (DTNSRDC) , Carderock  Laboratory,  was  requested  by 
the  Naval  Sea  Systems  Command  (NAVSEA)  to  develop  a computer 
program  that  can  be  used  to  better  design  and  predict  the 
propulsive  and  cavitation  performance  of  contrarotating 
propellers.  Contrarotating  propellers  have  been  used  on 
naval  vessels,  especially  torpedoes,  because  they  offer 
advantages  over  single-screw  propellers  by  being  more 
efficient,  having  smaller  optimum  diameters,  and  being  torque- 
balanced  resulting  in  better  stability.  Unfortunately, 
the  contrarotating  propeller  theory  and  design  procedure 
have  not  kept  pace  with  the  advancement  developed  for  single- 
screw propellers  where  the  design  and  perforihance  predictions 

can  be  made  with  a high  degree  of  accuracy. 


The  old  contrarotating  propeller  design  procedure  used 
at  DTKSRDC  is  based  on  Letbs*  theory  of  References  1,  2, 
and  3.  This  method  requires  that  the  same  hydrodynamic 
pitch  distribution  be  specified  as  input  for  the  forward 
and  aft  propellers  and  results  obtained  show  that  similar 

circulation  distributions  are  computed  for  both  propellers. 

% 

The  average  axial  and  tangential  propeller  induced  velocities 
used  to  determine  the  propeller  interaction  effects  are 
derived  using  the  uniformly  loaded  sink  disc  theory  in  this 
design  method. 

DTNSRDC  is  presently  evaluating  a package  of  contrarotating 
propeller  design  computer  programs  based  on  lifting  line 
and  lifting  surface  theories#  recently  developed  by  Nelson 
in  Reference  4.  The  propeller  lifting  line  theory  used 
requires  the  circulation  distribution  (which  must  be  the 
same  for  the  forward  and  aft  propellers)  rather  than  the 
hydrodjmamic  pitch  distribution,  to  be  specified  as  input. 

The  corresponding  hydrodynamic  pitch  distribution  is  calculated 
using  Lerbs'  moderately  loaded  single  screw  theory 


1.  Lerbs,  H.W.,  "Contra-Rotating  Optimum  Propellers  Operating 
in  a Radially  Non-Uniform  Wake,"  David  Taylor  Model 
Basin  Report  941,  May  1955 

2.  Morgan,  William  B.  and  Wrench,  J.W.,  Jr.,  "Some  Com- 
putational Aspects  of  Propeller  Design,"  Methods  in 
Computational  Physics,  Vol.  i,  academic  Press  Inc., 

New  York,  p 301-331,  19,65 

3.  Morgan,  W.B.,,  "The  Design  of  contraro^ting  Propellers 
Using  Lerbs’  Theory,"  Transactions- of  the  Society  of 
Naval  Architects  and  Marine  Engineers,  vol.  68,  p 6-38, 
I960 

4.  Nelson,  D.M. , "A  Computer  Program  Package  for  Designing 
Wake-Adapted  Coun ter rojtating' -Propellers:  A Users  Manual," 
Naval  Undersea  Center,  Fleet  Engineering  Department 
Report  NUC  TP  494,  December  1975 


of  Reference  5,  which  has  been  extended  to  account  for  finite 
circulation  values  at  the  propeller  hub.  Nelson's  method  for 
determining  propeller  interaction  effects  (average  axial  and 
tangential  propeller  induced  velocities)  are  determined  using 
the  procedure  developed  by  Hough  and  Ordway  in  Reference  6, 
corrected  to  account  for  finite  blade  number  effects.  Pre- 
liminary results  from  Nelson's  computer  program  for  making 
pitch  and  camber  calculations  based  on  contrarotating  pro- 
peller lifting  surface  theory  show  that  lifting  surface  effects 
due  to  the  forward  and  aft  propellers  are  small. 

The  new  contrarotating  propeller  design  ccxsputer  program 
based  on  Letbs  moderately  loaded  single  screw  propeller  theory 
of  Reference  5 is  similar  to  the  computer  program  developed 
for  single  screw  propellers  in  Reference  7 except  that 
additional  calculations  are  required  to  account  for  propeller 
interaction  effects  required  in  the  contrarotating  propeller 
design  program.  The  average  axial  and  tangential  propeller- 

induced  velocities  needed  to  determine  the  propeller  interaction 

5.  Lerbs,  H.W, , "Moderately  Loaded  Propellers  with  a Finite 
Number  of  Blades  and  an  Arbitrary  Distribution  of  Circulation, 
Transactions  of  the  Society  of  Naval  Architects  and  Marine 
Engineers,.  Vol.  60,  p 73-111,  1952 

6.  Hough,  G.R.  and  Ordway,  D.Ei,  "The  Generalized  Actuator 
Disk,"  Advanced  Research  Report  TAR-TR-6401,  Therm,  Inc., 
January  1964 

7.  Caster,  E.B.,  Diskin,  J.A. , and  LaFone,  T.A.,  "A  Lifting 

Line  Computer  Program  for  the  Preliminary  Design  of  Propellers 
David  W.  Taylor  Nayal  Ship  Research  and  Development  Center 
Report  SPp-!595-bl,  November  1975 


effects  are  derived  using  Kerwin's  field  point  velocity 
program  described  by  Denny  in  Reference  8.  These  calculations 
represent  the  most  important  new  feature  of  the  design  procedure 
presented.  Unlike  the  old  method  ' * t the  new  contra- 
rotating propeller  design  coaster  progrM  allows  different 
hydrodynamic  pitch  distributions,  wake  distributions  and  rpm 
values  to  be  specified  as  input  for  the  forward  and  aft  propellers. 
Nelson's  li^'ting  surface  theory^  for  contrarotating  propellers  should 
be  used  to  calculate  the  final  pitch  and  carsber  for  these 
propellers.  Preliminary  results  using  Nelson's  lifting  surface 
theory  show  that  the  lifting  surface  interaction  effects  due 
to  the  forward  and  aft  propellers  are  small.  As  a result, 
lifting  surface  theory  developed  for  single  screw  propellers 
(Reference  9)  may  be  used  to  determine  the  final  pitch  and 
camber  for  the  contrarotating  propellers  if  Nelson's  program 

is  not  available.  The  estimated  propeller  stresses  are- 

10 

calculated  based  on  simple  beam  theory  modified  to  account 
for  effect  of  rake  and  skew.  The  propeller  weight,  spacing 
between  propeller  blades,  chord;  lengths  for  lifting  surface 

pitch  and  camber  calculations  and-  blade  load  distributions 

8.  Denny,  Stephen  B.,  "Comparisons  of  Experimentally  Determined 
and  Theoretically  Predicted  Pressures  in.  the  Vicinity  of  a 
Marine  PropeUer,"  Naval  Ship  Research  and  Development 
Center  Report  23.49,  May  1967 

9.  Morgan,  W.B.,  Silovic,  Vladimir,  and  Denny,  S.B.,  "Pro- 
peller Lifting  Surfaqe  Corrections,^  Transactions  of  the 
Society  of  Naval  .Architects  and  Marine  Engineers,  Vpl.  76, 
p 3p9-347,  1968 

10.  Eckhardt,  M.K.  and  Morgan,  W.B.,  "A  jPropeller  Design  Method^" 
-Transactions  of  the  Society  of  Naval  Architects  and- Marine 
Engineers,  Vol>.  63,.  p 325,-;370,  .1955 
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for  finite  element  stress  calculations  are  also  included  in  the 

calculations  made  using  this  method.  Parameters  for  cal- 

culating  the  minimum  expanded  area  ratio  by  Xeller'^  and 

12  13 

the  methods  of  Burr ill  and  Brockett  are  used  to  check 
the  cavitation  performance  of  each  propeller.  An  input  option 
to  specify  the  hull  radial  induced  velocity  (not  available  in 
the  program  presented)  would  also  inorove  the  accuracy  of  the 
design  calculations. 


PROPELLER  LIFTING-LIHE  THK)RY 
As  mentioned  in  the  Introduction,  the  new  contrarotating 
design  procedure  follows  closely  Lerbs*  moderately  loaded, 
single-screw  lifting-line  propeller  theory  as  described  in 

O 

Reference  7.  Kerwin's  field  point  velocity  program  is  used 
to  determine  the  propeller  interaction  effects  (average  axial 
and  tangential  propeller-induced  velocities) . The  diameter 
of  the  aft  propeller  is  based  on  contraction  of  the  slip 
stream.  These  calculations  are  made  using  the  continuity  equation  as 

described  in  Reference  1,  once  Lerbs'  distance  factors  (g  ) 
values,  plotted  in  Figure  1,  are  specified.  More  exact  con- 

traction  of  the  slip'  streeun  calculations  can  be  made  using 

11.  Keller,  J.  Auf'm,  "Enige  aspecteri  bij  het  ontwerpen  van 
Scheepsschroeven,"  SChip  eh  Werf,  No.  24,  p 658-662,  1966 

12.  BurriU,  L.C.  and  Emerson,  A.;  "Propeller  Cavitation; 

Further  Tests  on  l6-Ihch  Propeller  Models  in  the-  Kings' 

College  Cavitation  Tunhel/i"  Transactions  of  the  North 
East  Coast  Institution  of  Engineers  and  Ship  Builders, 

•V61.  78,  p 295-320,  1963t64- 

13.  Brockett,  Terry,  "Minimum'' Pir'essure' Envelopes  for  Modified 
NACA  66  Sections  with  NACA  a»0.8:Camber  and'  BUSHIPS  Ty^  1 
and  Type  II  Sections,"  David  Taylor  Model  Basin  Report 
1780,  1966 
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the  propeller-induced  field  point  velocity  progcM  of 
Reference  8 if  velocities  are  calculated  at  a sufficient 
mmber  of  points.  This  approach  was  not  used  in  the  pro- 
gram presented  L^cause  it  would  result  in  a significant 
increase  in  the  core  size  and  running  tine  of  the 
computer  program. 

The  main  steps  in  the  new  contrarotating  propeller 
design  procedure  presented  are  as  follows: 

1.  Lerbs*  moderately  loaded  single-screw  lifting- 

5 

line  propeller  theory  is  used  to  make  design 
calculatons  for  the  forward  propeller.  The 
average  axial  and  tangential  velocities  induced 
by  the  aft  propeller  on  the  forward  propeller 
are  not  included  initially  in  these  calculations^ 
but  must  be  included  in  subsequent  calculations. 

2.  Lerbs'  contrarotating  (equivalent)  propeller 
design  procedure^  is  used  to  calculate  the 
aft  propeller  diameter. 

3.  Kerwin's  field  point  velocity  program  described 
in  Reference  8 is  used  to  compute  the  average 
axial  and  tangential  velocities  induced  by  the 
forward  propeller  on  the  aft  propeller. 
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The  aft  propeller  design  calculations  are  atade 
using  Lerbs*  moderately  loaded  single-screw 
lifting-line  propeller  theory^  where  the  average 
axial  and  tangential  velocities  Induced  by  the 
forward  propeller  on  the  aft  propeller  (step  3) 
are  included  in  the  aft  propeller  calculations. 
These  steps  (1  through  4)  are  repeated  until 
the  propeller-induced  velocities  converge. 


DESCRIPTION  OP  INPUT  DATA 

Appendix  A presents  a list  of  dimensioned  propeller 
design  parameters  required  for  the  new  computer  program 
developed  which  must  be  in  the  international  system  of 
units  (SI  units) . This  list  also  contains  conversion 
factors  (Kgj.)  changing  dimensioned  parameters  from  SI 
units  to  English  units.  Effective  power  (P^) / speed  (V) , number 

of  blades  (Z) , diameter  (D) , propeller  wake  / the 

hydrodynamic  flow  angle  distribution  (Bj)  and  distance 
factors  (g  ) from  Figure  1 are  required  input  parameters 
in  order  ‘to  make  nonviscous  propeller  design  calculations 
based  on  lifting-line  theory.  The  radial  distribution  of 
blade  chord  lengths,  npndimensionalized  on  diameter,  (c/D) 
and  section  drag  coefficients  (C^)  must  be  specified  as 
input  if  design  calculations  are  to  account  for  the  viscous 
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drag  effects  on  the  blades.  The  radial  distribution  of 
naxinum  thickness  nondimensionalized  on  chord  length  [t/c) , 
blade  rake  angle  and  nke«#  angle  (Qg)  parameters  are 

input  so  propeller  stresses  based  on  beam  theory  of 
Reference  10  can  be  calculated.  The  static  head  (H)  is 
input  so  the  blade  section  cavitation  number  (o)  can  be 
computed.  Appendix  A also  gives  the  complete  input  fortaat 
(card  numbers/  format  and  description  of  input  pareuneters) 
for  the  computer  program.  A brief  description  of  how  some 
of  these  parameters  can  be  determined  will  be  discussed 
next. 


Effective  Power#  Speed,  and  Shaft  Power 

Effective  power  and  speed  are  normally  obtained  from 

model  self-propulsion  experiments.  Input  effective  power 

(P_)  and  shaft  power  (P_)  are  defined  as  follows: 

£ S 


Pg=VT(l-t) 

(1) 

Pg=«2TtnQ 

(2) 

where  n =*  propeller  revolutions  i 

P_  = shaft  power,  5 

f* 

P_  » effective  power,  \ 

Cl  s 

Q » propeller  torque,  j 

T = propeller  thrust, 

V * ship  speed,  and 

(1-t)  » thrust  deduction,  which  may  vary  with 
propeller  diameter  and  speed. 


Nondimensional  Radial  Distance  (x) 


This  is  a reference  set  of  eleven  nondimensional  radial 
distances  at  which  all  other  distributions,  either  input 
or  calculated  by  the  computer,  are  defined  as  existing. 

In  general,  x.  * t ^/R,  with  the  restrictions 

*1  ■ "h''® 

= R/R=l, 

where  r^^  = the  distance  along  the  propeller  reference  line 
from  the  shaft  axis  to  the  ith  section, 
rj^  = propeller  hub  radius,  and 
R = propeller  tip  radius. 

Propeller  Wake 

The  radial  distribution  of  the  axial  wake  (1-w  ) which 
varies  with  propeller  diameter  is  also  required  input 
data.  The  circumferential  mean  of  the  axial  velocity 
distribution  d~w^)  is  obtained  from  a wake  survey 
without  the  propeller  operating.  However,  the  (l-^w  ) 
wake  distribution  must  be  corrected  for  the  propeller 
action.  No  completely  satisfactory  method  is  presently 
available  to  obtain  this  correction,  but  an  approximation 
of  the  radial  distribution  of  the  wake  with' 

the  propeller  operating  is  obtained  as  follows: 


Wake  distcibution 


a-\)  - ^l'W^)(l-w^>)/(l-V  (3) 

where  (1-w  ) » radial  distribution  of  the  circumferential 
c 

mean  wake  from  wake  survey  data, 

(1-w^)  * effective  wake  from  self  propulsion  data 

2 3. 

(l-w^)  » volume  mean  wake,  [2/(l-Xj^) 

h 

R * propeller  radius, 

r « propeller  local  radius, 

r^  » propeller  hub  radius, 

X « nondimensional  radial  distance  (r/R) , and 

X.  * nondimensional  hub  radius  (r. /R) . 
h h 

The  propeller  wake  distribution  may  also  vary  with  propeller 
diameter  depending  on  the  hull  characteristics  of  the  vessel. 

Advance  Angle  Distribution  Option 

The  advance  angle  distribution  (tan3)  defined  as 
V(l-w^)/(TrnDx)  is  normally  calculated  on  the  computer  for 
the  case  where  the  propeller  wake  fro*"  Equation  (3) 

sufficiently  represents  wake  in  the  plane  of.  the  propeller 
being  designed.  For  most  single  screw  propeller  designs 
this  approach  gives  good,  performance  predictions.  If  a 
propeller  operates  inside  a duct  or  in  the  vicinity  of  another 
propeller  as  in  the  case  of  tandem  or  contrarotating  propellers, 
the  axial  (w  /V)  and  tangential  (w. /V) . velocities  induced 
by  these  additional  sources  can  be  accounted  for  using 
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different  methods.  For  contrarotating  propellers,  the 
propeller  interference  velocities  are  computed  using  Kerwin's 
field  point  velocity  program  described  in  Reference  F in 
the  following  manner: 

tanP=  [{1-w J+w^/v]/r(x/Xj-w./v)  (4> 

where  w^/V  = axial  velocity  induced  by  forward  and  aft 
propellers  on  each  other, 
w^A  = tangential  velocity  induced  by  forward  and 
aft  propeller  on  each  other, 

= advance  ratio  based  on  ship  speed,  V/(imD) 

V = ship  speed,  and 

D * propeller  diameter. 

It  can  be  seen  from  Equation  (4)  that  for  the  case  where 
(w  A)  and  (w  A)  values  are  specified  as  zero,  the  advance 
angle  tan0  is  calculated  in  the  usual  manner  when  designing 
single  screw  propellers. 

Hydrodynamic  Flow  Angle 

The  hydrodynamic  flow  angle  distribution  (tanB^)  can 
be  specified  as  input.  An  option  is  included  so  Lerbs* 
optimum  tanBj  distribution^®  can  be  calculated  by  the  ccxn- 


puter  as  follows: 


(5) 


tan$j  » (tan8/n^)  ( 

where  * propeller  ideal  efficiency 

tanB  = advance  angle  distribution. 

Lerbs'  optimum  tauB^  distribution  usually  results  in  optimum 
propeller  efficiency.  If  other  factors  such  as  cavitation, 
strength  and  vibration  are  considered,  the  input  of  an 
alternate  tanBj  distribution  may  be  desired. 

Static  Head 

The  static  head  (H)  at  the  shaft  centerline 

is  required  input.  This  parameter  (H)  is  defined  as 

H +H  -H  , where  H is  the  shaft  submergence,  H is  the 
s a V s ^ a 

atmospheric  pressure,  and  is  the  vapor  pressure  of  fluid 
which  is  normally  small  compared  with  H and  may  be  neglected. 
The  static  head  (H)  is  used  to  calculate  the  section  cavitation 
number  (o)  in  Equation  (25)  and  the  Burr  ill  cavitation  number 
^ of  Equation  (30) . 

Blade  Outline  and  Expanded  Area  Ratio 

The  blade  outline  (c/D)  and  expanded  area  ratio  (A^/A^) 
must  be  input  for  the  design.  An  expanded  area  ratio  (A^/A^) 
is  calculated  on  the  computer  according  to; 


Expanded  Area  Ratio: 

V^O  ^ ^2Z/ti)/^  c/D  dx  (6) 

h 

where  c = chord  length, 

c/D  = nondlmensional  chord  length 
Z = number  of  blades 

The  final  blade  outline  and  expanded  area  ratio  should 
be  chosen  to  give  satisfactory  propeller  strength  and  cav- 
itation characteristics. 

Blade  Thickness  to  Chord  Ratio 

The  input  of  maximum  thickness  to  chord  ratio  (t/c) 
values  allow  an  estimate  of  the  propeller  principal  stresses 
(see  The  Propeller  Stress  Calculations  Using  Beam  Theory 
section  discussed  later)  based  on  beam  theory^®  to  be 
calculated  during  the  preliminary  design  stage  of  the  pro- 
pellers. From  a rough  estimate  of  the  blade  outline  (c/D) 
for  the  final  design  and  an  estimate  of  the  radial  distribution 
of  thickness  (t/D)  based  on  fatigue  strength^^,  the  following 
equation  can  be  used  to  obtain  initial  (t/c)  input  values: 

Blade  Thickness  Ratio: 

t/c  * (t/D) /(c/D)  (7) 

where  t/D  = radial  distribution  of  thickness  (can  be  estimated 
from  Reference  10. 
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Rake  and  Skew 


The  rake  angle  at  the  blade  tip  (0^^)  and  the  skew  angles  (6^) 
for  a design  are  specified  to  permit  adequate  predictions 
of  principal  propeller  stresses  using  the  beeun  theory  method 
described  in  Reference  10  and  discussed  later  in  the  Propeller 
Stress  section  of  this  report. 

The  rake  (6.)  is  defined  consistent  with  Reference  14 
as  the  distance  from  the  propeller  plane  to  the  generator 
line  in  the  direction  of  the  shaft  axis.  Aft  displacement 
is  considered  positive  rake. 

Since  the  skew  angles  {9g)  significantly  affect  propeller 
unsteady  forces,  a computer  program  based  on  the  unsteady 
contrarotating  propeller  lifting  surface  theory  of  Reference  15 
can  be  used  to  select  the  skew  angles  (3g)  for  the  design. 

The  input  skew  angles  (0g)  in  degrees  are  defined  as  the 
angular  displacement  of  points  on  the  blade  reference  line 
from  the  propeller  reference  line  in  the  projected  view. 


14.  Cumming,  R.A.,  Dictionary  of  Ship  Hydrodynamics  - Propeller 
Section,  14th  International  Towing  Tank  Conference 

1975,  Report  of  Presentation  Committee,  Appendix  VII,  1975 

15.  Tsakonas,  J.  and  Jacobs,  W.R.,  "Counterrotating  and 
Tandem  Propellers  Operating  in  Specially  Varying,  Three- 
Dimensional  Plow  Fields,"  Davidson  Laboratory,  Stevens 
Institute  of  Technology  Report  1335,  September  1968 
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Section  Drag  Coefficient 


In  order  to  account  for  viscous  effects  when  predicting 
the  performance  of  a propeller,  the  section  drag  coefficient 
(C^)  must  be  specified  as  input.  A section  drag  coefficient 
(C^)  value  of  0.0085  usually  gives  reasonable  estimates 
of  model  propeller  drag  for  propeller  shapes  normally  used  at 
DTNSEIDC  in  the  past.  For  propellers  having  very  thick  blades, 
the  following  equation,  available  as  an  input  option  on 
the  computer,  and  derived  as  a function  of  maximum  thickness 
(t/c)  values  using  experimental  data -from  NACA  66  type 
section^^'^' , will  give  a better  estimate  of  the  section 
drag  coefficient  (C^) ; 

Section  Drag  Coefficient: 

= CpQ  [l+1.25(t/0  + 125  (t/c)^j  'B) 

where  C„  is  the  frictional  resistance  of  the  section,  e.g., 
FO 

C - « 0.008  for  Reynolds  number  of  approximately  10^  and 
r O 

g 

0.004  for  Reynolds  number  of  approximately  10  . 

Options  for  using  alternate  nonlinear  distributions,  or 
a constant  distributions  are  also  available. 

Lerbs  Axial  Distance  Factors 

It  was  noted  earlier  that  Lerbs'  axial  distance  factors 
(9g) / rather  than  the  use  of  the  more  correct  propeller- 

induced  velocities  from  Reference  8,  were  used  to  obtain 

16,  Abbot,  Ira  H.  and  Von  Doenhpff,  Albert  E.,  "Theory 

of  Wing  Sections  Including  a Summary  of  Airfoil  Data," 
Dover  Publication  Inc.,  New  York,  Library  of  Congress 
Catalog  No.:  60-1601,  1949 

17.  Hoerner,  S.F.,  "Fluid-Dynamic  Drag,"  Published  by  the 
author.  Midland  Park,  New  Jersey,  1965 


contraction  of  the  slip  streaai  and  the  aft  propeller  diameter 
in  order  to  nsinimize  core  size  and  running  time  of  the  con^uter 
program.  The  g values  from  Reference  1,  plotted  as  a function 

A 

of  propeller  radius  and  spacing  between  propeller  blade 
center  lines,  are  presented  in  Figure  3. 


DESCRIPTION  OF  OOTPOT  DATA 
Lerbs*  lifting  line  theory  is  used  to  calculate  the 
propeller  lift  coefficient  (C_),  nondiraensional  circulation 
(G) , hydrodynamic  flow  angle  (6j),  axial  induced  velocity 
(U^/2V) , and  tangential  induced  velocity  (U^^2V) . These 
lifting-line  calculations  take  into  account  viscous  drag 
effects  on  the  propeller  by  specifying  as  input  in  the 
computer  program  the  propeller  section  nondimensicnal  chord 
length  (c/D)  and  section  drag  coefficient  (C^) . A method 
for  obtaining  values  for  c/D  and  is  discussed  in  the 
Description  of  Input  Data  section  of  this  report.  With 
these  parameters  available,  the  necessary  design  and  per- 
formance prediction  parameters  for  contrarotating  propellers 
can  be  obtained. 

Thrust  and  Power  Loading  Coefficients,  and  Propulsive  Efficiency 
The  new  contrarotating  propeller  design  computer  program 
calculates  the  thrust  (C^g)  and  power  (Cpg) / loading  coefficients 
based  on  ship  speed  and  the  estimated  propulsive  efficiency  (n^) 
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foe  these  propellers  In  the  following  manner: 
Thrust  I«oadiMg  Coefficient: 


^fwd,aft 

e 


(l-etanSj)/(dC^g^/dx)dx»T/(  (p/2)tirV) 

h 


"=T=>Cli  - '=Ts>f»a  " ""Ts'att 

Power  Loading  Coefficient; 

(Cps  aft*-^x  (1+e/tanSj)  (dC  /dx)dx=Pg/(  (p/2)1!rV) 

e ' h 

^Ss^CR~^Ss^fwd  ^^^aft^^fwd^  ^Ss^aft 
Propulsive  Efficiency: 

>o.=a-t)  VCp3=V^„ 

where  based  on  ship  speed 

‘Ss>desigh'^^C<0^^>”“''^’ 


fwd,aft  = subscripts  for  forward  and.  aft  propellers» 
respectively 

CR  = subscript  for  the  set  of  contrarotating  propellers 


C . = nondimens ional  inviscid  local  thrust  loading 

coefficient,  4ZG(x/Xg-U^/2V) 

C . = nondiraens ional  inviscid  local  power  loading 

coefficient,  (4Z/Xg)xG(  (1-Wj^)+U,j,/2V) ) 

G = nondimens ional  circulation  from  lifting  line  theory 

U^/2V  = axial  Induced  velocity  from  lifting  line  theory 

U^2V  = tangential  induced  velocity  from  lifting  line  theory 

V = ship  speed 

p * density  of  fluid 


(9) 

(10) 

(11) 

(12) 

(13) 
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The  calculated  thrust  (T)  is  obtained  fro*  Squation  (9) 
and  the  shaft  power  (Pg)  is  calculated  using  Bquatlon  (11) 
for  each  propeller. 

Other  paraneters  useful  in  designing  and  evauating 
the  performance  of  propellers  include  the  advance  coefficient 
(J) , ship  speed  advance  coefficient  (J^) , thrust  coefficient 
(K^) , torque  coefficient  (K^) , leoisent  due  to  thrust  (M^) , 
moment  due  to  torque  , moment  parallel  to  section  nose- 
tail  line  (M^^) , Boment  perpendicular  to  the  nose-tail  line 
(My©) * the  blade  loading  distribution  (LI) . These 
parameters  are  calculated  as  follows i 
Advance  Coefficientt 

J*V(l-w^)/(nD)«V^/(nD) 

Ship  Speed  Advance  Coefficient: 

Jy»V/(nD) 

Thrust  Coefficient: 

Ky»T/(pnV).(TtC^g/8)Jy^ 

Torque  Coefficient: 

Kg-Q/{pn^D^) «(Cpg/16) 

Manent  Due  to  Thrust: 

Mjjj(x)  «(pi|R^y^/.(2Z)  (x-x^),  (^-etanBj)  (dC^g^/dx)dx 

h 


(14) 


(15) 


(16) 


(17) 


(18) 


Moment  Due  to  Torque: 


Hqjj(x)-(p7IrV/(2Z)  {tanBj+e)[dC^^/dx]dx 


(19) 


Moments  Parallel  to  Section  Nose-tail  Line: 

(20) 

Moment  Perpendicular  to  Section  Nose-tail  Line: 

Blade  Loading  Distribution: 

LI(x)=JjPcV^^  C,  (22) 

L ii 

where  x,x^=propeller  nondimensional  radial  stations,  r/R  and 
4)  = propeller  pitch  angle  (<j)*6j  in  computer  program) 

s=  section  inflow  velocity,  /•.  (1"W^)  +U^/2v)'^+(  (x/Xg)  -U^2v]^ 
PROPELLER  STRESS  CALCULATIONS  USING  BEAM  THEORY 

A propeller  blade  must  contain  enough  material  to  keep  the 
stresses  within  a blade  below  a certain  predetermined  level. 

This  level  depends  on  the  material  properties  with  regard 
to  both  steady-state  and  fatigue  strength  and  to  both  mean 
and  unsteady  blade  loading.  The  material  selection  controls 
the  allowable  stress  level  and  the  blade  chord,  thickness, 
rake  and  skew  are  the  main  parameters  which  control  the 
blade  stress  for  a given  blade  loading.  Stresses  based 
on  beam  theory^^  in  the  propeller  blade  are  computed  for 
each  propeller.  Both  hydrodynamic  and  centrifugal  loadings 


are  considered.  Effects  of  rake  and  skew  are  included 


In  this  stress  calculation  procedure/  the  propeller  blade 
is  represented  as  a straight  cantilever  beas  of  variable 
cross-section  without  caaber.  Experimental  results  show 
that  the  neutral  axis  of  an  airfoil  section  lies  approxiisately 
along  the  nean  line#  so  caaber  is  not  considered  in  tlie 
stress  calculations  presented.  Only  the  maximum  principal 
stresses  calculated  at  the  mid-chord  of  each  section  ace 
printed  as  output  in  the  cocoputer  program.  Stcessej;  for 
the  final  design  should  be  calculated  by  finite  element 
techniques  if  rake  and  skew  for  the  propeller  differ 
from  usual  propeller  shapes. 

PARAMETERS  FOR  MAKING  BLADE  SURFACE  CAVITATION  CHECKS 

Brockett's  theoretically  derived  incipient  cavitation 

charts  of  Reference  13  can  be  used  to  predict  the  blade 

surface  cavitation  characteristics  of  each  propeller  once 

the  lifting-line  calculations  have  been  completed.  The 

two-dimensional  caraber-to-chord  ratio  (fjj/c) » ideal  angle 

of  attack  in  degrees  (a^) t section  cavitation  number  (0) 

nondiraensionalized  with  the  section  inflow  velocity  (V^) , 

and  the  .maximum  and  minimum  fluctuating  angles  of  attack 

(a  fOi  . ) in  degrees  are  parameters  that  must  be  determine^ 
max  min 

before  Brockett's  'incipient  cavitation  charts  can  be  usedi 
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These  parameters  are  calculated  as  follows: 

Section  r^axiiaum  Camber  to  Chord  Ratio  for  NACA  a»0.8  Heanline: 

f^c  * O.OSTSCj^  (23) 

Section  Ideal  Angle  of  Attack  in  Degrees  for  NACA  a>0.S  Meanline: 

* 1.54  C^  (24) 

Section  Cavitation  Number: 

O = 2g(H-xR)A^^  (25) 

where  g = acceleration  due  to  gravity 

H » static  head  at  shaft  centerline  (see  section  on 
Static  Head) 

= inflow  velocity  of  each  propeller  section. 

The  maximum  and  minimum  fluctuating  angles  of  attack 

(a  . ) in  degrees  are  calculated  using  the  method  derived 
in3x  tnin 

by  Lerbs  and  Rader  in  Reference  18.  These  calculations 
can  be  made  using  the  following  equations: 

Maximum  Fluctuating  Angles  of  Attack: 

'W  = 'i  - 

Minimum  Fluctuating  Angles  of  Attack: 

\in  “ ®i  ■ (27) 

where  -AB  = maximum  effective  angle  of  attack  in  degrees 
(from  wake  survey  data) 

+A6  minimum  effective  angle  of  attack  in  degrees 
(from  wake  survey  data) 

18.  Lerbs,  H.Wi  and  JUtder,  H.P.-r  "Uber  det  Auf.tr.iebsgradienten 
von  Profiien  im  Propeller  Verbahdj"  Schiffstechnik^ 

Vol.  9,  No.  48,  p i78-180,  1962 


22. 


The  paraneter  F(x)  in  Equations  (26)  and  (27)  is 
dependent  on  the  hydrodynamic  flow  angle  (B^) , the  advance 
angle  (8)  and  the  lift  coefficient  (C. ),  and  is  calculated 
on  the  computer  using  the  following  equation: 

P(x)  - l/(l+2irtan{8j-B)/Cj^)  (28) 

12 

The  Burrill  Cavitation  Charts  can  also  be  used  to  give 
an  approximate  check  on  the  cavitation  performance  of  pro- 
pellers. Burrill *s  thrust  loading  coefficient  (t^)  and 
cavitation  number  (0^  ^)  at  0.7  radius,  defined  as  follows, 
are  parameters  that  must  be  known  in  order  to  use  these 
cavitation  charts: 

Burrill  Loading  Coefficient: 

* T/{(p/2)Ap(v(l-Wj^,p  .^))^+(0.7TTnD)^}  (29) 

Burrill  Cavitation  Number  at  0.7  Radius: 

« 2gH/{(v(l-w^^jj  .y))2+(0.7TrnD)2}  (30) 

where  A_  » propeller  expanded  area,  Z/  c dr 


Aq  * propeller 
Ap  = propeller 

(P/D)  j.  * estimated 
estimated 

Keller’s  method^^  of 

ratio  of  the  propeller  is 


2 

disc  area,  TiD  /4 

projected  area,  (l.067  - 0.229(P/D)  )Ag 

propeller  pitch  ratio  at  0.7  radius, 
as  0.7Trtan8j 

predicting  the  minimum  expanded  area 
also  calculated  on  the  computer. 
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.Copy  cvcdloble  to  DDC  'does  not 
pen^t JuUr  reproduction 


The  minimum  expanded  area  ratio,  based  on  eliminating  back 
bubble  type  cavitation,  is  computed  as  follows: 

‘Wk  ^ 0.6Z)K^{Oq  .^(j^  + (0.7lT)^)}+K  (31) 

where  K = 0.2  is  used  for  single  screw  ships  with  manganese- 
bronze  propellers  having  rake  of  approximately  10  degrees, 

K = O.IO  is  used  for  twin-screw  ships  with  copper-aluminum- 
bronze  propellers,  K = 0.15  is  used  for  twin-screw  ships 
with  manganese-bronze  propellers,  and  K = 0 to  K - 0.05 
is  used  for  propellers  for  fast  ships  such  as  destroyers 
and  frigates.  The  constant,  K = 0.15,  was  used  in  developing 
this  computer  program.  Tf  a different  value  of  K is  desired, 
the  expanded  area  ratio  calculated  in  Equation  (31)  should 
be  adjusted  to  account  for  changes  in  the  value  of  K. 


Chord  Lengths  for  Lifting  Surface  Pitch  and  Camber  Calculations 
The  final  pitch  and  camber  for  each  propeller  can  be 
calculated  using  Nelson's  computer  program,  presented  in 
Reference  4,  based  on  lifting  surface  theory  for  contrarotating 
propellers.  However,  Nelson  has  shown  that  in  some  cases 
lifting-surface  interaction  effects  of  contrarotating  propellers 
are  small.  If  this  is  the  case,  single  propeller  lifting 
surface  theory  from  References  19,  20,  21,  and 

19.  Kerwin,  J.E.,  "The  Solution  of  Propeller  Lifting-Surface  Problems 
by  Vortex  Lattice  Methods,"  Department  of  Naval  Architecture  and 
Marine  Engineering,  Massachusetts  institute  of  Technology, 
Cambridge,  Massachusetts,  1961 

20.  Kerwin,  J.E.  and  Leopold,  R. , "Propeller  Incidence  Correction 
Due  to  Blade  Thickness,"  Journal  of  Ship  Research,  Vol.  7,  No.  2, 
1963,  p 1,  6 

21.  Cheng,  H.M.,  "Hydrodynamic  Aspects  of  Propeller  Design  Based  On 
Lifting-Surface  Theory:  Part  I - Uniform  Chordwise  Load  Dis- 
tribution," David  Taylor  Model  Basin  Report  1802,  1964 
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22  can  be  used  to  calculate  the  final  pitch  and  camber  for 
each  propeller.  These  progreuns  require  as  input  the  section 
chord  length/  (c/R),_  and  (c/R)__  nondimensionaiized  on 
propeller  radius,  in  terms  of  the  skew  angle  (0^) , hydrodynamic 
flow  angle  (6j)  and  blade  outline  (c/D) . The  parameter 
(c/R)j^g  and  (c/R)^j,  measured  from  the  leading  and  trailing 
edges  to  its  reference  line,  respectively,  are  calculated 
in  the  following  manner  on  the  computer: 

Chord  Lengths  Measured  from  Blade  Leading  Edge: 

(c/R),„  a xO  /{57.296  cosB.)  - c/D  (32) 

Xic*  S Jl 

Chord  Lengths  Measured  from  Blade  Trailing  Edge: 

(c/R)__  * x0  /(57.293  cosB.)  + c/D  (33) 

Tc  S X 

SPACING  BETWEcaJ  BLADES  AND  FILLETS 

Propeller  designs  should  have  enough  clearance  between 

23 

blades  at  the  hub  so  fillets  can  be  properly  applied.  Hill 
derived  the  following  equation  which  is  used  in  the  program 
to  estimate  spacing  between  blades  at  the  hub  without 
fillets: 

Gjj  » (2iirj^)/Z  - (tj^/sin  (J>)  (34) 

where  (|)  is  the  propeller  pitch  angle  (4»  * 8^  in  computer  program) . 
Based  on  a number  of  full-scale  propellers  built  with  standard 
fillets.  Hill's  blade  clearance  equation  was  modified 

the  following  manner  to  estimate  spacing  G^  between  fillets 

22.  Cheng,  H.M.,  "Hydrodynamic  Aspects  of  Propeller  Design 

Based  On  Lifting-Surface  Theory:  Part  li  - Arbitrary  Chordwise 
Load  Distribution,"  David  Taylor  Model  Basin  report  1803,  1965 

23.  Hill,  J.G.,  "The  Design  of  Propellers,"  Transactions 

of  the  Society  of  Naval  Architects  and  Marine  engineers, 

Vol.  57,  1949 
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at  the  hub  during  the  preliminary  stage  of  the  design. 

Gp  = {2nrj^)/Z  - (1.9tj^/sin4»)  (35) 

A layout  of  blade  sections  is  recommended  as  a final  fillet 
clearance  check. 

PROPELLER  WEIGHT  AND  CENTER  OF  GRAVITY 

The  approximate  propeller  weight  (Wp)  and  location 
of  center  of  gravity  (CG)  from  the  propeller  center  line 
is  also  calculated  for  each  design.  In  order  to  make 
these  calculations,  the  specific  weight  of  the  material 
(Dp)  must  be  specified  as  input.  The  propeller  hub  is 
assumed  to  be  a circular  cylinder  of  equal  length  and 
diameter.  The  propeller  center  line  is  located  at  the  mid 
length  of  the  propeller  hub  and  no  allowance  is  made  for 
the  propeller  bore  or  blade  root  fillets  in  these  calculations. 

The  weight  (Wp)  for  each  propeller  is  calculated  as  follows; 

Propeller  Weight; 

Wp  = Wg  + W„  (36) 

where  W 

n 

Wjj  = weight  of  hub,  tv'4(P^Djj) 

Djj  = hub  diameter 

a(x)  = area  of  section  at  radius  x,  2c(x)  t(x)T  Jt(x,x^)  dx 

c(x)  • chord  length  of  section  at  radius  x 

t(x)  = maxlmuj(»  thickness  of  section  thickness 


weight  of  blades,  fuz/"  a(x)dx 
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t(x,Xj^)  » chocdwise  distribution  of  section  thickness 

(NACA  66  nodified  thickness  form  is  used  in  program) 
Xj^  * nondiraensional  coordinate  along  the  section  chord 
(rj/c) 

Pp  » specific  weight  of  material  considered 

The  propeller  center  of  gravity  (CG)  with  respect  to  the 
blade  center  line,  where  plus  values  represent  the  distance 
ahead  of  the  center  line  and  negative  values  aft  of  the 
center  line,  is  calculated  in  the  following  manner: 

Center  of  Gravity: 

CG  = Mp/Mp  (37) 

where  Mp  = moment  of  the  propeller,  Zf^  a(x)  B(x)  dx 

h 

B(x)  * distance  of  center  of  gravity  from  propeller 

center  line,  y cos4»  + x sin^i  - 9^  x R tan<|)  - 9„  x + d„/2 

X * longitudinal  center  of  gravity  along  chord  line 

y = vertical  center  of  gravity  perpendicular  to  chord  line 

9p  = rake  angle  in  radians 

6 skew  angle  in  radians 

s 

({)  » pitch  angle  in  radians,  is  used  in  the  program) 
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COMPUTER  PROGRAM 

A new  computer  program  has  been  derived  for  the  preliminary 
design  of  contrarotating  propellers  using  the  CDC  6400,  6600  and 
6700  ccxnputers  at  DTNSRDC.  A core  storage  of  approximately 


27 


128,000  octal  is  required  for  the  computer  program,  and  it  takes 

approximately  4 minutes  to  compile  the  computer  program. 

The  actual  running  time  per  case  using  the  design  thrust 

\ 

option  is  approximately  6 minutes  and  when  the  design  shaft 
horsepower  option  is  used,  the  running  time  is  approximately 
9 minutes.  A detailed  description  of  the  input  and  output 
formats  for  the  computer  program  is  presented  in  Appendix 
A,  and  a FORTRAN  listing  of  the  new  computer  program  is 
shown  in  Appendix  B. 

PROPELLER  DESIGN  THRUST  AND  POWER  OPTIONS 
The  thrust  option  can  be  used  to  make  lifting  line 
calculations  for  propellers  required  to  produce  a given 
thrust  at  specified  values  of  speed  and  rpm  (this  is  accomplished 

by  adjusting  the  input  tanSj  distribution),  or  the  power 
option  can  be  used  if  the  propeller  is  required  to  absorb 
a specified  power  at  a given  rpm  (in  which  case  the  speed 
is  determined) . 

From  each  calculated  power,  a new  value  of  speed  (assumed 
to  vary  as  the  cube  toot  of  the  ratio  of  the  design  and 
calculated  pover)  is  obtained  and  its  corresponding  effective 
power  is  obtained  from  the  effective  power  input  curve. 

Design  calculations  again  produce  a new  calculated  power, 
and  the  process  continues  until  the  closeness  criteria  of 
design  and  calculated  power  is  satisfied  (two  iterations 
are  normally  sufficient) . Smaller  increments  of  input  speeds 
in  general  cause  faster  convergence. 
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Once  the  basic  shape  of  the  distribution  is  defined 
(see  the  Hydrodynamic  Flow  Angle  Distribution  section)  the 
final  tanB^  distribution  tanBj  is  determined  using  the 
thrust  or  power  options,  by  malcing  lifting  line  calculations 
of  three  nondimensional  thrust  loading  coefficients 
^^Ts^2'  ^^Ts^3  three  hydrodynamic 

pitch  distributions,  Kj^  tanB^,  K2  *^3  tanBj  where 

Kj^  = 0.975,  K2  * 1.0,  and  “ 1.025.  Once  these  calculations 
are  obtained,  the  following  system  of  equations  are  set  up: 

"^Ts'l  = » ♦ “1  * CkJ 

(Ct,1j  = a ♦ BKj  . CK^ 

= a * 3K3  . 

from  which  values  of  A,  B,  and  C are  obtained.  Then,  values 

of  A,  B,  and  C are  substituted  in  the  following  equation 

to  obtain  the  value  of  K.. 

4 

Ss  = ^ ^ ^*^4 

where  Ss  = design  thrust  loading  coefficient  (Equation  9). 

DESIGN  CALCULATIONS  OF  A SAMPLE  SET  OF  CONTRAROTATING  PROPELLERS 
Design  calculations  are  shown  in  Table  1 for  a sample 
set  of  contrarotating  propellers  with  7 blades  and  an  expanded 
area  ratio  (Ag/A^)  of  0.293  for  the  forward  propeller,  and 
6 blades  with  A^/Aq  = 0.365  for  the  aft  propeller.  These 
propellers  were  designed  to  operate  in  a wake  with  a thrust 
loading  coefficient  based  on  ship  speed  (C^g)  of  0.265  and 
an  advance  coefficient  (J)  of  1,235.  Both  propellers  had 
a 100  percent  linear  skew  distribution  resulting  in  51.4 


and  60  degree  blade  tip  skew  angles  for  the  forward  and 

aft  propellers,  respectively.  Solid  and  dashed  lines  in 

Figures  2 through  8 show  calculations  obtained  using  the 
12  3 

old  design  method  ' ' and  the  new  computer  program, 
respectively.  These  figures  show  that  the  radial  distribution 
of  the  hydrodynamic  pitch  (nxtanS^) , nondimens ional  circulation 
(G) , axial  (U^2V)  and  tangential  (U^2V)  induced  velocities 
calculated  using  both  methods  differed  significantly, 
especially  for  the  aft  propeller.  The  thrust  loading  coefficient 
(C^g) , power  loading  coefficient  (Cp^)  and  propulsive 
efficiency  (n^)  calculated  for  the  sample  set  of  contra- 
rotating propellers  using  the  old  and  new  methods  are  printed 
in  Table  2.  Results  in  Table  2 show  that  the  old  design 
method  predicts  similar  thrust  (C^g)  and  power  (Cp^)  loading 
coefficients  for  the  forward  and  aft  propellers  where  as 
the  new  design  method  predicts  that  the  aft  propeller 
and  Cpg  values  are  significantly  higher  than  those  calculated 
for  the  forward  propeller.  The  propulsive  efficiency  (n^) 
of  0.886  predicted  for  the  sample  set  of  contrarotating  propellers 
using  the  new  method  is  approximately  2 percentage  points 
lower  than  the  ri^  value  of  0.912  calculated  using  the  old 
design  method.  Performance  predictions  using  the  new  design 
method  should  be  more  accurate  than  the  predictions  made 
using  the  old  design  method  because  of  the  improved  methods 
of  accounting  for  the  propeller  interaction  effects. 
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RECOMMENDATIONS 

The  following  investigations  should  be  carried  out  in  the 
future  in  order  to  better  design  and  predict  the  performance 
of  contrarotating  propeller*: 

1.  Utilize  induced  velocities  derived  by  Kerwin  in 
Reference  8 to  obtain  more  exact  calculations  for 
contraction  of  the  slip  stream  and  the  aft  propeller 
diameter.  This  subroutine  should  be  reprogramed  to 
significantly  reduce  the  running  time  presently 
required  to  make  it  practical  for  use  in  the  contra- 
rotating propeller  design  computer  program. 

2.  Include  the  effects  of  the  hull  radial  induced 
velocity  as  an  input  option  in  the  contrarotating 
propeller  design  procedure. 

3.  Conduct  model  propulsive  and  cavitation  experiments 
on  a set  of  contrarotating  propellers  to  validate  the 
design  procedure. 

4.  Options  for  specifying  noncylindrical  hubs,  arbitrary 
propeller  locations  and  alloiwances  for  the  propeller 
bore  should  be  added  to  the  program  presented  to 
improve  the  predictions  of  propeller  weight  and  center 
of  gravity. 
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Figure  2 - Hydrodynamic  Pitch  Angle  Distribution  for  Sample  Set  of  Propellers 
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Figure  4 - Circulation  Distributions  Calculated  for  the  Aft  Propeller 
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BURItlLL  THUOST  COEFF  TC«  .8«i*«*t<l 
SUR^ILL  CkVITITtOH  COEFF  SI6HRC*.T»«  .1Z««C*«I 
CIEMINCE  IT  K0»  BETMEEM  BL*0ES/0>  -.8285B3B5 

CLEIIIHCE  IT  HOB  BETHEEH  Fl£.UTS/0»  -.81491592 


TABLE  2 


THRUST  AND  POWER  LOADING  COEFFICIENTS  AND  PROPULSIVE  BF7ICIOICY 
CALCULATED  USING  THE  OLD  AND  NEW  DESIGN  METHODS 


Perfocnance  Predictions 

Old  Method^' 

New  Method 

Thrust  Loading  Coefficient 
^Ss^  fwd 

0.137 

0.133 

Thrust  Loading  Coefficient 
aft 

0.137 

0.141 

■Thrust  Loading  Coefficient 
^^Ts^  CR 

0.265 

0.265 

Power  Loading  Coefficient 
^^Ps^  fwd 

0.113 

0.108 

Power  Loading  Coefficient 
<^Ps^aft 

0.116 

0.128 

Power  Loading  Coefficient 
^^Ps^  CR 

0.222 

0.229 

Propulsive  Efficiency 

0.912 

0.886 

45 


TABLE  1 


List  o£  dimensioned  i^put  and  output  parameters  used  by 
computer  progreun  based  on  English  and  SI  units 


Parameter 

English  Units 

SI  Units 

KSI^^^ 

Shaft  power  (Pg) 

hp 

KH 

0.7457 

Effective  power  (P  ) 
£« 

hp 

KW 

0.7457 

% 

Ibm/ft^ 

kg/m^ 

16.01846 

V 

knots 

m/^ 

0.514444 

V 

ft/sec 

knots 

0.5924 

D 

ft 

m 

0.3048 

H 

ft 

m 

0.3048 

P 

Ibf  sec^/ft^ 

kg/m^ 

515.3788 

n 

rev/rain 

(3) 

rev/min  1.0 

T,  weight 

Ibf 

N 

4.44822 

V 

r 

ft/sec 

ra/s 

0.3048 

LI 

Ibf /ft 

N/m 

14.5939 

^b 

in  Ibf 

Nn 

0.112985 

Qb 

in  Ibf 

Nn 

0.112985 

M 

xo 

in  Ibf 

Nn 

0.112985 

M 

yo 

in  Ibf 

Nn 

0.112985 

Max  Stress 

Ibf/in^ 

Pa 

6894.757 

SKEW 

deg 

deg<3> 

1.0 

RAKE 

deg 

deg^^^ 

1.0 

Mass  polar  moment 
of  INERTIA 

Ibm  in^ 

2 

kgm 

0.00029264 

(1)  Multiply  English  Units  by  KSI  to  get  SI  Units, 

(2)  Computer  program  uses  knots  in  both  English  Units  option 
and  SI  Units  option. 

(3)  These  are  not  SI  Units  but  are  permitted  to  be  used  in  the 

SI  system  according  to  International  Standards  Organization  (ISO) 
Standard  No.  1000. 


APPENDIX  A 
INPUT  FORMAT  (FORNARD 
Card  No.  Format  Parameter 

1 P8.6  d/R 

2 12A6  

3 F8.4  Pg  or  0.0 

4 F8.4  1.0  or  5.0 

Pp 

5 F8.4  V 

6 F8.4  Pg 


PROPELLER} 

Description  of  Input 

Axial  distance  between  the 
forward  and  aft  propeller 
planes  nondime ns ionali zed 
on  forward  propeller  radius. 

Identification  of  design 
input  data. 

Input  the  design  shaft  power 
if  the  power  option  is  used. 
Input  0.0  if  the  design 
thrust  option  is  used  (see 
the  section  on  the  propeller 
design  thrust  and  power 
options) . 

Number  of  different  speeds 
(V)  and  effective  power  values 
to  be  punched  on  Cards  5 
and  6,  respectively.  Input 
1.0  for  the  thrust  option 
and  5.0  for  the  power  option. 

Specific  weight  of  propeller 
material  (Pp)  for  weight 
calculations. 

Design  r;peed  (V) . For  the 
power  option,  input  five 
different  speeds  in  increasing 
order  with  the  third  speed 
being  the  best  estimate 
of  final  speed.  Input  one 
speed  for  the  thrust  option. 

Design  effective  power  values 
corresponding  to  input  speeds 
on  Card  5.  (See  Equation 
(1) ). . 
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mm 


Card  No. 
7 


Format 

P8.6 


Parameter 

D 


(1-V 


(1-t) 

H 


Description  o£  Input 
Propeller  diameter  (D) 

Effective  wake  (1-w^) 

Thrust  deduction  (1-t) 

Static  head  (H)  at  shaft 
centerline.  (See  section 
on  Static  Head) . 


Density  of  water  (p) 


Skew  Angle 
(8^)  Option 


For  a linear  skew  angle 
distribution  input  percent 
of  sicew  in  decimal  form. 

For  nonlinear  skew  input 
-1.0  and  the  skew  angles 
in  degrees  (8  ) must  be 
punched  on  Cards  22  and 
23. 


F8.6 


1.0 

C„ 


Rake  angle  at  blade  tip 
in  degrees.  Use  +8„  value 
for  blades  raked  aft  and 
-8_  for  forward  rake. 

K 

Input  1.0 

Use  CD>10  to  input  the  radial 
distribution  of  drag  coefficients 
(CD) ; 0<CD<10  to  input  a constant 
drag  Cp^CD  at  all  radial  stations; 
CO^O  causes  the  computer  program 
to  calculate  the  radial  distribu- 
tion. of  drag  coefficients  using 
the  equation  Cp»0.008(l+1.25(t/c) 
+125(t/c) '*);  -10<CD<0  causes  the 
computer  program  to  use  a constant 
frictional  resistance  Cp<y=ABS (Cpo) 
in  Equation  ^ C0»Cpo(l+1.25(t/c) 
+125(t/c) '*);  CEK-IO  to  input  the 
radial  distribution  of  frictional 
resistance  (CFO),  values  to  be  used 
in  Equation  8. 


tahS  Option  Input  0.0  if  tan8  is  calculated 

0.0  or  1.0  on  the  computer  in  the  normal 

manner.  Input  1.0  if  tan3 
d^tribution  is  arbitrary 
and  punched  on  Cards  26 
and  27. 


•«9, 


^ '•  ‘ ' • V. » J ' c-''-',-.-  ♦ . 


Card  No. 

Format 

Parameter 

Description  of  Input 

9 

F8.6 

Z 

Number  of  blades 

10 

F8.6 

V‘o 

Expanded  area  ratio 

11 

F6.6 

n 

Revolutions  per  unit  time 

12-13 

F8.6 

X 

Nondiraensional  radial  station 
(eleven  values  arbitrary 
spaced  from  propeller  hub 
to  tip) 

14-15 

F8,6 

(1-w^) 

Propeller  wake  (l“Wj^)  cor- 
responding to  input  X values 

16-17 

F8.6 

tanBj 

Tangent  of  hydrodynamic 
flow  angle  (tanSj)  corresponding 
to  input  X values.  If  Lerbs 
optimum  tanS  distribution 
is  desired,  0.0  values  must 
be  punched  in  these  cards 

18-19 

F8.6 

c/D 

Section  chord  lengths  correspond 
ing  to  input  x values 

20-21 

F8.6 

t/c 

Thickness-chord  ratios 
corresponding  to  input  x 
values 

22-23 

F8.6 

Optional 

Blade  skew  angles  corresponding 
to  input  X values  to  be 
input  only  if  -1.0  for  skew 
option  is  punched  in  Card 
7.  Otherwise  Cards  22  and 
23  must  not  be  used  as  input 
cards. 

24-25 

F8.6 

Cjj  Optional 

Section  drag  coefficient  cor- 
responding to  input  values  if 

O)>10  on  Card  8;  frictional  re- 
sistance of  section  if  Cp<10. 


5 b 


I 


Car*i  Mo.  Format 
26-27  F8.6 


28-29  F8.6 


Parameter  Deacription  of  Input 

tanS  Optional  Tangent  of  advance  angle 

corresponding  to  input  x 
values  to  be  input  only 
If  1.0  for  tanfl  option  is 
punched  on  Card  8.  Other- 
wise Cards  26  and  27  must 
not  be  used  as  input. 

q Distance  factors  (g  ) 

* corresponding  to  input  x 

values  from  Reference  1 
and  Figure  1. 


INPUT  FORMAT  {AFT  PROPELLER) 


Card  No.  Format  Parameter 

30  12A6 

31  F8.4  Pg  or  0.0 


32  F8.4  1.0  or  5.0 


Pp 


33  F8.4  V 


34  F8.4  P_ 

c 

35  P8.6  0.0 


(1-Wj.) 


Description  of  Input 

Identification  of  design 
input  data. 

Input  the  design  shaft  power 
delivered  at  the  propeller 
if  the  power  option  is  used. 
Input  O.Oif  the  design  thrust 
option  is  used  (see  the 
section  on  the  propeller 
design  thrust  and  power 
options) . 

Number  of  different  speeds 
(V)  and  effective  power 
values  to  be  punched  on 
Cards  33  and  34,  respectively. 
Input  1.0  for  the  thrust 
option  and  5.0  for  the  power 
option. 

Specific  weight  of  propeller 
material  ip^)  for  weight 
calculations. 

Design  speed  (V) . For  the 
power  option.  Input  five 
different  speeds  in  increasing 
order  with  the  third  speed 
being  the  best  estimate 
of  final  speed.  Input  one 
speed  for  the  thrust  option. 

Design  effective  power  values 
corresponding  to  input  speeds 
on  Card  33.  (See  Equation  (1)). 

Input  0.0  for  aft  propeller 
diameter  which  is  calculated 
on  the  computer. 

Effective  wake  (1-w^) 

Thrust  deduction  (1-t) 
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Foraat 


F8.6 


F6.6 


P«r«Mt«r 

R 


P 

Skew  Angle 

(0  ) Option 

8 


9 


R 


1.0 


Description  of  Input 

Static  Read  (H)  at  shaft 
centerline.  (See  section 
on  Static  Head) . 

Density  of  water  (o) . 

For  a linear  skew  angle 
distribution  input  percent 
of  skew  in  decieal  fore. 
For  nonlinear  skew  input 
-I.O  and  the  skew  angles 
in  degrees  (9  ) must  be 
punched  in  Cards  50  and 
51. 

Rake  angle  at  blade  tip 
in  degrees.  Use  +9_  value 
for  blades  raked  att  and 
-0j^  for  forward  rake. 

Input  1.0 


Cjj  Use  CD^IO  to  input  the  radial 

distribution  of  drag  coefficients 
(CD) ; 0<CD<10  to  input  a constant 
drag  Co^CD  at  all  radial  stations; 
CD»0  causes  the  computer  program 
to  calculate  the  radial  distribu- 
tion of  drag  coefficients  using 
the  equation  Cd=0.00S[i+1.25( (t/c) 
+125(t/c)  "*) ; -10<CD<0  causes  the 
computer  program  to  use  a constant 
frictional  resistance  Cpo*ABS(Cpo) 
in  Equation  8,  CD“CFo(ltl.25(t/c) 
+125  (t/c)**);  CD£-10  to  input  the 
radial  distribution  of  frictional 
resistance  (CFO)  values  ‘o  be  used 
in  Equation  8. 


tang' 'Opt ion 
0.0' or  1.0 


Z 


Input  6.0  if  tang  is  calculated 
on  the  computer  in  the  normal 
manner.  Input  1.0  if  tang 
distribution  is  arbitrary 
and;^puhched  on  Cards  54 
and  55. 

Nuiber  of  blades 


Expanded  area  ratio 


PS. 6 


Card  No. 

Format 

Parameters 

Description  of  Input 

39 

P8.6 

n 

Revolutions  per  unit  time 

40-41 

P8.6 

X 

Nondimens ional  radial  station 
(eleven  values  arbitrarily 
spaced  from  propeller  hub 
to  tip) . 

42-43 

P8.6 

(1-W  ) 

X 

Propeller  wake  (1-w  ) cor- 
responding to  input^x  values. 

44-45 

P8.6 

tanSj 

Tangent  of  hydrodynamic 
flow  angle  (tar.*^^)  corresponding 
to  input  X values.  If  Lerbs 
optimum  tanS  distribution 
i;i  desired,  6.0  values  must 
bi.'  punched  in  these  cards. 

46-47 

P8.6 

c/D 

Section  chord  lengths  corresponding 
to  input  X values. 

48-49 

P8.6 

t/c 

Tl.ickness-chord  ratios 
corresponding  to  input  x 
values. 

50-51 

F8.6 

6^  Optional 

Blade  skew  angles  corresponding 
to  input  X values  to  be 
input  only  if  -1.0  for  skew 
option  is  punched  in  Card 
35.  Otherwise  Cards  50  and 
51  must  not  be  used  as  input 
cards. 

52-53 

P8.6 

Cp  Optional 

Section  drag  coefficient  cor- 
responding to  input  values  if 
CD>10  on  Card  8;  frictional  re- 
sistance of  section  if  CD<lp. 

54-55 

P8.6 

V 

tanB  Optional 

Tangent  of  advance  angle 
corresponding  to  input  x 
values  to  be  input  only 
if  1.0  for  tan$  option  is 
punched  on  Card  8.  Otherwise 
Cards  54  and  55  must  not 
be  used  as  input. 

TABLE  I 


A DESCRIPTION  OF  THE  OUTPUT  PARAMETTERS 


Computer 

Parameter 

Parameter 

Description  of  Output 

V 

V 

Speed  (input) 

PE 

Effective  power  (input) 

D 

D 

Diameter  (input) 

Z 

Z 

Number  of  blades  (input) 

H 

II 

Static  head  (input) 

RHO 

P 

Density  of  fluid  (input) 

AE/AO 

V'‘o 

Expanded  area  ratio  (input  or 
Equation  (6)) 

CPS  I 

Nonviscous  power  coefficient 
(Equation  (11)  when  e=«0) 

CTSI 

Nonviscous  thrust  coefficient 
(Equation  (9)  when  e»0) 

XI 

X 

Nondimens ional  radii  (x) 

DCTSI 

dC 

Tsi 

Local  nonviscous  thrust  coefficient 
(Defined  in  Equation  (9)) 

DCPSI 

"Ssi 

Local  nonviscous  power  coefficient 
(Defined  in  Equation  (11)) 

CPS 

^PS 

Power  loading  coefficient  (Equation  (11)) 

CTS 

Ss 

Thrust  loading,  coefficient  (Equation  (9)) 

TANBI 

tanBj 

Tangent  of  hydrodynamic  flow  angle 

TANB 

tan8 

Tangent  of  advance  angle 

G 

G 

Nondimensional  circulation 

UT/2V 

U/2V 

Tangential  velocity  induced,  at 
lifting  line 

VA/2V 

V2V 

Axial  velocity  induced  at  lifting 
1 
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Computer 

Parameter 


Parameter 


Description  of  Output 


VR 


V 

r 


Section  inflow  velocity,  Equation  (22) 


SIGMAX  0 


Section  cavitation  number,  Equation  (25) 


CL 

ALI 


a. 


FM/C 


CD/CL 

F(X) 


C 

F(x) 


LI 


LI 


SKEW  ANGLE  S 
(C/R) LE 


s 

(c/R) 


LE 


(C/R) TE 


(c/R) 


TE 


T/R 

PC 


t/R 

n,. 


PS 


l-T 


(1-t) 


Section  lift  coefficient 

Section  two-dimensional  ideal  angle 
of  attaclr  in  degrees  for  NACA  a=0.8 
meanline,  Equation  (24) 

Section  two-dimensional  maximum  camber 
ratio  for  NACA  a=»0.8  raeanline. 

Equation  (23) 

Section  drag-lift  ratio  ( "Cp/C^^) 

Parameter  for  calculations  section 
fluctuating  angles  of  attack. 

Equation  (28) 

Propeller  load  distribution  for  finite 
element  stress  calculations.  Equation  (22) 

Propeller  skew  angle  in  degrees 

Chord  lengths  for  lifting  surface 
pitch  and  camber  calculations, 

Equation  (3  2) 

Chord  lengths  for  lifting  surface 
pitch  and  camber  calculations, 

Equation  (3  3) 

Ratio  of  section  thickness  to  radius 

Estimated  propulsive  efficiency, 

Equation  (13) 

calculated  shaft  power  delivered 
at  the  propeller,  Equation  (13) 

Propeller  thrust  deduction 

Propeller  effective  wake. 


1-WT 


Computer 

Parameter 

Parameter 

Description  of  Output 

TH{N) 

T 

Design  thrust/  Equation  (9) 

TH(N) 

T 

Calculated  thrust.  Equation  (9) 

A 

Area  of  section 

XBAR 

Longitudinal  position  about  x axis 
parallel  to  nose-tail  line  from 
centroid 

YBAR 

Vertical  distance  about  y axis 
perpendicular  to  nose-tail  line 
from  centroid 

IXO 

Moment  of  inertia  about  x axis 
parrallel  to  nose-tail  line 

lYO 

Moment  of  inertia  about  y axis 
perpendicular  to  nose-tail  line 

MTB 

«Tb 

Bending  moment  due  to  thrust. 
Equation  (18) 

MQB 

Bending  moment  due  to  torque, 
Equation  (19) 

MXO 

M 

xo 

Bending  moment  about  the  x axis, 
Equation  (20) 

MYO 

M 

yo 

Bending  moment  about  the  y axis, 
Equation  (21) 

MAX  STRESS 

Maximum  stress,  Reference  10 

RAKE  ANGLE 

0R 

Rake  angle  in  degrees 

PIXTANBX 

TPCtanPj 

Hydrodynamic  pitch  angle 

PXXTANB 

irxtan3 

Advance  angle 

HBXGHT  OF 
BLADES 

Equation  (36) , 
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Equation  (36) 


WEIGHT  OF  (W  ) 

PROP  (BLADES  ^ 

+ CYLINDRICAL 
HUB) 

CENTER  OF  CG  Equation  (37) 

GRAVITY  OF 

PROP 

CENTER  OF  Equation  (37),  m -0 

GRAVITY  OF 
BLADES 

KELLERS  Equation  (31) 

MINIMUM  EAR  » K 

SPEED  J Equation  (15) 

COEPF(JV) 

ADVANCE  J Equation  (14) 

COEFF  (JA) 

THRUST  K_  Equation  (16) 

COEFF  (KT) 

TORQUE  K Equation  (17) 

COEFF  (KQ)  ® 

PROPULSIVE  n Equation  (13) 

EFFICIENCY  ° 

ETAD 

BURRILL  T Equation  (29) 

THRUST  ^ 

COEFF  (TC) 

BURRILL  a Equation  (30) 

CAVITATION 

COEFF 

CLEARANCE  G Equation  (34) 

AT  HUB 

BETWEEN 

BLADES 


Conputer 

Paraaeter 

Parameter 

Description  of  Output 

CLBAMNCE  AT 
HUB  BSnfBHN 
FILLBTS 

Equation  C)***) 

CTS 

CR 

Thrust  loadlnq  coefficient  for  aet  of  conbrs- 
rotating  propellers  (Equation  10} 

CPS 

(C.,g)cR 

Power  loading  coefficient  for  aet  of  contra- 
rotating propellers  (Equation  12} 

BTAO 

^VcR 

Propulsive  efficiency  for  set  of  contra- 
rotating propellers  (Equation  12) 

PROGRAt-  GMAJN  (INPUT, OUTPUT, TAPE5sIKPUT,TAPFG  = OUTPUTI 
COHHON/CHEIGHT/X,CHORD,ThIC<MS,CAHBE-?,»»ITCH,S<FMR,OIAM,ZZ,OEM,tA«e 
l.ST 

lfPI,PP7,PPe,PP9,PPll,EHAKt , V$,FPS, SIGMA, £AR,3T 
DIHENSICN  9<3d,3«J)  ,01  (181)  ,92nai)  ,0T(  181)  ,o«»  (36?)  , 35(  362)  , 3G  (362) 
1 ,3F(  161)  ,08(181)  ,09(181)  ,810(181)  ,011  (181)  ,012(181)  , 013(181)  ,B1<*( 
2191) ,01 6(161) ,aZ (38,^8), OH (33,18) ,C(  58, 38),CC( 38,38) , INDEX  138,3) 
3,A(38, Z2) ,P9(38, 72) 

COMMON  B,ei,O2,0  3,Pt,,86,86,9F,O8,n9,OlO,dll,B12,O13,914,B16,AZ,0H, 
1C,CC, INDEX, A, 89 
COMMON  lOjJB, JC, JO, JOD, JEf 
COMMON  CL  1(11) 

01 men  SIGN  CHORC( 59 ) , TH IC XNS ( T6 ) , C AMBER ( Jd) , PI T CH ( 38)  ,S<EMR(38) 
1,X(^6) ,BrTAI(39)  ,91(11) 

0IHEN9ICN  X3(  38)  ,y«4(3fc)  ,X5(38) , X6133)  , V£L(  •))  ,EHP(9)  ,6LA(9)  ,EXX(9) 
OI-EMSICrj  AZZ(39,-<e) 
riMENSICN  AbHP(9) 

OIMFNSICN  XMH(a) 

OIMEN'^ICK-  CAtf  (,)  ,CAf  (9) 

OIMENSICU  FX(  11)  ,9?  J( 11) 

•DIM«-NSICN  VFLl  (9)  ,f  hoi  (9) 

OIMENSICN  L(11,14,?) 

OIMFNSTCN  CONST ( 11) ,o ( 11) ,H( 11) 
dimension  SVL («) ,S9£ (9) 

OIMENSICN  UA  ( 11)  ,UT  (1 1 ) ,><V(1!  ) ,UR  ( 11)  ,UAl  { 11 ,2)  ,UT1(  11 ,2)  ,UR1  ( 11,2 
1 ) 

OTMFNSIC“  AtX  (11,9) 

OIMcNSIGN  OEX  (11 ,3,2)  ,COM(  12,2) 

CIMINSION  FH0SHP13)  ,AFTchO(3) 

OIMENSICN  GZ(lll) 

DIHENSICN  AX(  11)  , PARI  11) 

DIMENSION  VSUOr.SOOl)  ,VSUuf-(il)  , »LPT  ( 1 1 ) , AV  ( 1 1 ) , 3V  ( 1 1 ) 

COMMON  CCONt(ll)  fCCTHOdl)  ,G:THP(11)  ,CCF0R(11) 

OIMe-NSION  AREA  (7)  , XJAP  (7)  , Y3AR  ( 7)  , A YEXO  ( 7)  , AYFY0(7)  , EMXC  ( 7)  ,EMY0(7 
1)  ,<^MT9(7)  ,EM0B(7)  ,<:T5'NAX(7) 

DIHENSICN  PXT3K  11)  ,PX-^R(11) 

01‘^ENSICN  SX(7) 

COMMON  FP1,PP2  ,PP3,PPA,FP5  ,PF6,Pr’10 
OIMFNSTCN  CAlKll) 

DIMENSION  HU0niM(6,2) 

DIHENSICN  CPT  (2)  ,rPS(2)  ,'*TS(2) 

COMMON  /HRT/  JPR 
9EAO(G,A06)  0EX(1,2,1) 

AOF  F0RMAT(r6.A) 

J3R=') 

IP0=1 
IOCR  = (' 

XDO=A. 0 
IDCaXOO 

PI  = 3.  IAI‘59265’50979 
OSI=16.018A6 
PWR=. 7656999 
VSI=. 5164664 
ELF=.3048 
€LI=.0254 
EL2=FLI*ELI 
EL4=EL2*FL? 


SIKX.112984S 
SMX=68<J'>.757 
SHT*U. 44822 
RH0SI=515.3788 
IPR*C 
JG*C 

00  813  K«l,2 

XRPHxl.O 

XZ7*i.C 

XE*=i,C 

JGxJG»l 

REA0(S«R3C)  (COHd, XI, 1x1,12) 
933  F0PMftT(1286) 

READ(E, ICOOC)  SPHR 

SHPsSPMR/PWR 

REAOJ^, lOaOO)  XVV,OSN  ■ 

HUB=C. 

OEN*OSN/OSI 

IRPHxXPPH 


7 


IVV=XVV 

REA0(5, 10009)  <SVL(I) ,1=1, IVV) 

PEAO(£,10009)  (SPE(I) ,I=1,IVV) 

REA0(5,  10009)  SOM, EHAKE .ETHRUS^SHOtSRO, XPS, RAKE 
00  409  i=l,IVV 
VEL(I)  = SVL(I)/VSI 
4C9  EHP(I)=SPE(I)/PHP 
OIAsSOM/ELF 
HEAO=SHO/ELF 


RHO=SPO/RHOSI 
IF(SHP)  103,102,103 
103  $HP»SHF/2.a 


NXVV=XVV 
1C?  CONTINUE 

PEAO(5,10C09) 

3<2,1)=C0 

JC=11 

IZZ=XZ7 

lEAxXEA 

PEAO(5, 10009) 

9LAS=3LA<1) 

PEAO<5,10C39) 

•:xxs=exx(i) 

PE  AO  (5, 10009) 
P£A0(5, 10009) 
PEA015, 10009) 
PEA0<5,1CC39) 
PEAO(5, 10009) 
PEAO<5,10CC9) 
IF(XPS)  6,7,7 
6 PEAO(5, 10009) 
00  26  1=1, JC 
AZZ(I,3fl)=AZZ 


TANBI,C0,TANB 


BLACI),I*1,IZZ) 

EXXCI),I=1,IEA) 

(XHM(I),I*l,IRPH) 
fX3fI),I=l,  JC) 
(X4(I),Izl,JC) 

(X5  a),i=i,jc) 

(XF(I),I*l,JC) 
CAZZfI,29) ,I=1,JC) 

(A7Z(I,24) ,I=1,JC) 

1,24) 


26  CONTINUE 
CONTINUE 

IF(ABS(C0).GE*1C.)  READ (9, 10089) 
IF«CO.r.E.lfl.  ) GO  TO  lOOiC 
TM*P(2,1) 


(Bn.,7)  ,I«1,.JC) 


1 


i 

i 

I 

i 

i 

i 


i 

! 


cro*.oo8 

IFCCO.GT.O.I  CFO»l. 

IF{CO.CT.>10.  .ANO«CO.LT.O«)  CFO«*CO 
00  10007  1*1, JC 

IF(CO.LE.O.I  TH*l.»i.?5*A2Za,?5l>l??.*AZ7II,25l**% 
IF(CO.LE.-10.  I CFP*8II,7I 
10007  B(I,7I«CF0*TH 
lOCC ' FOPMAT  C72H 

1 ) 
tOOOQ  F0PHAT(P8.U) 

10009  FOPHAT  C9F8.4) 

10010  IFCTANB.LE.O.)  PEAOI9,10009I  (BIZ, 8) , 1*1, JO 
IF(IOCP.NE.O)  GO  TO  57 

IFCK.NL.l)  GO  TO  57 
oFAn(5,10C09»  (GAIK  I)  ,r*t,ilJ 
5''  CONTINUE 

IFIHUR.NE.C.  » PFAO(5, 100091  ( HU«OIH ( I ,K ) , T= 1, 6» 

01 10,1,<1 =XVV 
P( 11,1 ,K) sXRPH 
00  «:2  1=1, IVV 
niI,l,tO=VEL(I) 

902  0( I,?,K» =EHP(I» 

0 (11 ,2,K)=EHAKF 
0ll,3,‘O=ETHRIJS 
0 (6,3,K>=XPS 
0(7,3,O=SHP 
0 (9,3,K»*CO 
0(1C,3,X)=TANB 
0( 1C,2,K)*0IA 
0(2,3,K»=HrAC 
0( 3,3, X) *PHO 
OI9,3,IO=XZ7 
C(E,3,K)»xeA 
0( «,3,X) *TANPI 
0(ll,3,*<»sQ£N 
n(ll,«,,K)=HUR 
00  e:3  1=1, izz 

903  0(  I,‘.,K)  =BLA(  II 
no  8C9  1=1, IFA 

909  0 (I,5,KI=EXX( I) 

00  8C5  I*l,I?PM 

805  0( I,fc,K) =XMM(I) 

00  8C8  1*1, JC 
OII,7,K)=X3m 
0( I,8,KI *X4(II 
nCI,9,K)=X5(I» 

0(I,10,K)*X6(It 

808  0(I,11,K)=AZZII,25I 
IF(XPSI  806,807,807 

806  00  125  1=1, JC 

125  0(I,12,X)*AZZII,?4> 

80  7 CONTINUE 

00  8C9  I=1,JC 

809  n(I,13,K)*B(I,7} 

810  IF(TANP»  812,812,911 

811  00  126  1=1, JC 

126  0(I,!4,'K)*Bil,8) 
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^■1 


i 


CONTINUe 

OEX(l(l»K)s4.0 

DEX(2,1,K)=11.0 

OEX{3»l«K)xi.O 

0EX(1,2,2»=-0EX{1,2,1) 

JTxll 

CONTINIJ*: 

IC*-2 

JM=0 

JQ=i) 

IK=1 

NF=0 

NTT=0 

IF(SHPI  90,91,90 

00  10073  NKsljNXVV 

CONTINUE 

NT=0 

NTT=NTT»l 
IF(NK.GT.l)  100*2 
00  10072  NIT*1,ID0 
NT=NT*1 
XVV=0 (11,1, IK) 

IV9SXVV 

XROM=0(11,1,IK) 

IF(SHP)  152,152,153 
IF(NTT.EQ.NXVtf4l)  GO  TO  151 
CONTINUE 
DO  020  1*1, IVV 
VEL(I)*0(I,1,IK) 

1 EHP(I)=0(I,2,IK) 

CONTINUE 

00  1 1*1, IVV 
CftV(I)*¥FL(I) 

CAE(I)*EHP(I) 

DIA=0 (10,2,IK) 
EHAKE=n(ll,2,IK) 
ETHRUS=0(1,3,IK) 
HEA0=0(Z,3,IK) 

RHOsO (3,3, IK) 

XZZ=0  (A, 3, IK) 

XEA=0(5,3,IK) 

XPS*0 (6, 3, IK) 

SHP=0(7,3,IK) 

TAN0I=0(8,3,IK) 

C0*0(9,3,IK) 

TAN8=O(10,3,IK) 
0EN=0(11^3,IK)-  . , 

HUR*0(1L,4,IK) 

00  821  1*1, IZZ 
8LA(I)*0(I,«*,IK) 

L 9LAS*BLA(l) 

00  822  1*1, lEA 
EXX(I)=0(I,5,IK) 

2 EXXSsFXX(l) 

no  823  I*1,IRPH 

3 Xf1N(I)*0(I,6,IK) 
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no  S2h  i=itjc 
x!(i) «nci,r,iK) 
xiidi  :oci»0,iio 
XS(T>  *0(1, 9,10 
X6(T»=OII,10,1*{) 

921.  ft2ZCI,?5)=n(l,ll  ,IK) 

IF(XPS)  b25,82fv,0?f« 

925  no  127  1=1, JU 

12?  AZZ(i,?4>  = n(l,12,Itf) 

8’6  CONTINUE 

927  00  129  1=1, JC 

129  d(l,7)=0(I,13,lKl 

929  IF(TAN9)  83l,d31,d?0 

8’5  90  129  1=1, JC 

129  BU,8)=0C1, 14,10 

931  CONTINUE 

'10  15  1=1, JC 

15  aZZd  ,23)  =v3m 
JC1=JC-1 

03  16  I=?,JC1 

16  AZZ(I,19)=X3(I.1) 

AZZ(l,19)=X3d) 
no  4 1 = 1, 11 
AZ‘'(x,36)=A2Z(I,-’5) 

4 AZZd  , 37)  =A?Z{  l,?3) 

00  10072  IP=1,IZZ 
3(9,2)=3LA  (IE) 

XSX  = XP«:*(360.0/P('J,2)) 
asi=xsx/<i,o-X3(i) ) 

A$’=XSX-ftSl 
00  10072  Ke=l,IFA 
no  10  072  Ii..P=l,IFPN 
PP<4=XP^(IRP) 

Ea-R  = FXX(t<F) 

00  100  1 = 1, JC 
=UT,3)  = XMI) 

B(T,4)  = X4  d) 

B(T,5)=X5(I) 

IF(XFS)  14,1-»,13 

13  CONTINUE 

AZZd  ,24)  =ASl*X3d)*9S2 
AZZd  ,24)  =0.0 
CO  TO  166 

14  AZZd,?4)=AZZd,39) 

166  CONTINUE 

AZZd  , 3fl)=AZZd,?4) 
ir?  CONTINUE 

DO  10051  LF=1,JC 

10051  9(LE,6)  = (9LAS»£AP*X6CU.n/(0(9,2)*eXXS) 
DO  10052  LF=l,JC 

10052  AZZ(LE,Z5)=AZZ(LE,36)*^BaE»9) 

00  30  1=1,9 

VEUI)  =CAVd) 

30  FHPd)=CAE(I) 

0(5,2)=XMHdRP)/6O,O 

IVVxXVV 

IVAsl 


€S 


iooiie 
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10001 

10002 
10003 


10012 

10017 


10018 

10019 


IF(SHPI  51,52,51 
DO  53  10*1, IVV 
VELKir.lsVFLIIG) 

6HF1(IG)*EHP(IGJ 

IV*1 

GO  TO  SH 
CONTINUE 
IV*1 

IFCNTT.EQ.NXWV+1)  GO  TO  155 

VELl(l)*VEL(NtO 

EHPm)aEHP(NK) 

GO  TO  5% 

V£L1(1)=VEL<1) 

EHPl(l)aEHPa) 

CONTINUE 

00  5 1=1,11 

AZZ<I,25)=AZZ<I,36) 

AZZ(I,?3)=ftZZ(I,3T) 

8(6,2>  = C325.86*£HPiaW))/(VELl  (IV>*ETHRUS) 

8(7,  2)  = 1.688*VEL1CIV1 

AJJ=1.0 

NN=0 

3J=1.0 

J8=8J 

FORMAT { 313) 

CONTINUE 

00j=1.0 

JE£=TAN31 

JOO=OOJ 

B6(361)=1.0 

J0=.666667*FL0AT(JC) 

3(l,l)=e.O 
3(2,1 )=CO 
8(T,1)=TAN8 
n(4,l)=0.6 
8(8,1)  = «.  0 
B(7,l)=g.O 
B(8,l)=0.0 
3(9,l)=0.0 
9(1,2)=1,0 
8(2,2)=1.0 
B(9,2)=RHO 
3(5,1)=EWAKE 
9(3,2l=OIA 
B(4,2)=H£A0 

RSL=B(7,2l/(3.mi59?65*B(5,2l*B(3,2l) 

IF  (8(1,5))  10017, 10017,, 10019 
CONTINUE 

IF(B(5,D)  800,800,801 
B(5,1)=6(JD,4) 

CONTINUE  - 
DO  100181=1, JC 

B(I,5.)=,CRSL*SQRI<(8(5  .,l)).TSQRTf  8 (1,4) ) > /(8(1, 3)  *8  (4,  D ) 
NN=NN*1 

HU9L=X3(1)*OIA/2,0 
A0JS=.(X3,(ll)^0IA/2.  0)  «HUBt. 

00  35  1*1,11 


‘A  — "I-  n'  '’‘,1  I 


35  RA<(I»-CCX3a)*0IA/2.  0>-HUbU  ♦RAKE/AO JS 
If»RsIPR41 

IF(JPP,NE.2)  go  TP  100A7 
50H=0IA*ELF 
WRITE (6,130) 

13C  FORHAT(lhl) 

IF(T<-1)  115,115,116 
115  HRITE(6,117) 

117  FOPMAT (ATX,*FO»HARL  PROP  OF  C OUTRAOOT AT ING  SET*) 

GO  TO  ll« 

HE  KRITE(6,119) 

119  F0R9AT (45X,*AFT  P»OP  OF  PONTPAkOTATING  StT*) 

11*  CONTISt'E 

WRITE  (6,830)  (CO** ( I , I K ) , I = 1, 1 2) 
i<;hp  = 2,*SFhr 

IF(ShF,lQ.O.)  WHTE(6,410) 

IF(SHP.NC.O.)  W-ITf <6,100?0)  TSHP 
41'’  FORMAT  (’X, 'THRUST  OPTION*) 

1002C  FCRMAT (7X,*PS (KK) =*, IX , IFF  12.4) 

WRITE (6, 1CC26)  XVV,OSN 

1032'^  FCPPAT  (2X,*N0  OF  V=*,F5.1,4X,  •DENSITY  DF  PROP 

1 (>(G/M3)  =*,Fn.2,4X) 

WRITE(6, 1C027)  (SVL( I) ,I=1,IW9) 

10327  FORMAT (?X,*V (M/lEC) =* , 1F9E12. 4) 

WRITE (6, 1C045)  (SPt ( I ) , 1= 1, I V°) 

1034E  FnRMAT(2X,*PE(KW)=*,2X,lPQF12.4) 

w:iTF  (6,  ^020)  SjH,EWAKE,ETH?.US,ShD,SRO,XPS,RAKE 
10025  FORMAT (2X,*D(M) =• F 3, 4 , 4X , *1- WTT =• , F 8. h ,4 X , * 1-THD=*F 8 . 4 ,4X , *H ( M)=* 

2 ,F9.4,4X,*RHO(<G/K3) =*F10.4/2X,*TETS  OPT=  • , Fb . 3 , 3X , *RAKE  OPT=  *, 
7 P4.3) 

WRI’-F  (P«,l  0C29)  TANGI,Cn,TAN9 

130  2=  format (2X,*TAN9I  OFT=  • , F4 , 2 , 4X , * DRAG  OPT=  • , F7. 4 , 4X , *T An9  0PT=  ♦, 
1F4.2) 

write (6,1 0 030)  (OLA (I) ,1  = 1, IZZ) 

100  30  poRMm-^  (2X,*?  = *,  10X,1P9:12.4> 

WRITE (^ , 1 03  31  ) (EXX  (1)  ,I  = 1,IEA) 

10031  FCRMAT(2X,*AE/AC=*6X,1P9E12.4) 

WRITE  (6,100  32)  (X'-M(  I)  ,1  = 1,  IRFM) 

1003?  PCRMAT (2X,*N(6EV/MIN) = *,1P9E12.4) 

WRIT£(6, 10042) 

WRIT F (6, 100  33)  (X  3(I),I  = 1,JC),  (X4(I) ,I  = 1,JC) , (X5(I) ,I=1,JC) , (X6(I) , 
11  = 1,  JO  , (AZZCI  ,2‘--)  ,1  = 1 ,JC) 

10  3 3 2 FORMAT (2X,*X=  * , 1P9E1 2, 4/ 1 OX , 1P2E 1 2. 4/2 X , *1- WX=  *,1P9E12.4/ 

110X,1P2E1?.4/2X,*TAN<3I=  * ,1?9F  12, 4/ 10 X , 1P2E  12.4/2X,  *C/0=  * , 1P9 

2E12.4/10X, 1P2E12.4/PX ,*T/C=  • , 1P9E1 2. 4/ 10 X , 1P2E12 . 4) 

IFCXPS)  ?2,?3,23 

22  WRITE  (f, 10036)  ( A Z7  (1 , 24)  , 1=1  , JO 

10036  FORMAT (2X,*TFTS=  * , iP9E 12. 4/1 OX, l^ZFie . 4) 

P3  CONTINIIF 

IF(CO)  24,24,25 

24  WRITE (6,10037) (R(I, 7) ,1=1, JC) 

10037  F0RHAT(?X,*C0=  *, 1P9E1 2, 4/10X» IP2E12. 4) 

25  CONTINUE 

10042  FORHATC4X,*0(M)  =*1PI;  10.4-,*  N(REV/MIN)  x*lPEl0. 4,*  Z=*,1P£10.4, 

2 * A£/A0=*iPE10,4  ) 

10044  F0RHATCJX,*ETA0=*, 

2 1PE10.4,2X,*PS(KW)**,1PE10. 4,3X,*1-THO=*,1P€10.4, 


Copy  avcdI(3blo  to  DDC  does  not 
pcrcii^  legiblo  reproduction 
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h' 


2 TX,*l-«TT**,lPF10.<it3X, 

2 •V(KNOTS»**lPfcl0.i*,3X,*  DESIGN  THtN) iP£10.%, 

3 ^,78X,*V(H/SeC)a*tlPE10*4,IX,*CALCUJ.ATEO  TH<M>**,1PE10.A» 
WRITE  (6tlQ0<i2) 

WRITE(6»iOO<i3)  S0H»RPH»e(<4«2)  »EAR 
10087  V<»8(7*2)/1.6878 

VMS=B<7t2»*ELF 
DO  10090  IsltJC 
AZZai25)=AZZ<I,?9l*8CI  ,6) 

10090  B(I,3Q)-6(I}5) 

1002A  F0RHAT(1P9E12.41 
IF(JN)  65f65»66 
65  CONTINUE 

XHHl*RPH/60.0 
IF(TANP)  i003<»»10034»i01 
1003A  DO  10035  I=1|JC 
10035  B<I,8»=RSL/0(I»3)‘»BCI,4I 

CM1=8C7,'8)"  ‘ ‘ 

101  CONTINUE 
GO  TO  67 
66  CONTINUE 

00  68  1=1, JC 

CONST (Il=PI*XHMl*DIA*X3(I)/8(7,2) 

6(I)=UA(I)fX4(I) 

H(IJ =CONST(II-UT (I) 

68  8(I,8)=G(I)/H(T» 

CH2sB(7,8) 

00  31  1=1,11 

31  B(I,8>=CW1/CH2*B(I,8> 

IF(TK-2I  33,32,33'' 

32  00  34  1=1,11 

34  B(I,8)=0(10,2,1)*6(I,S)/C(10,2,2) 

33  CONTINUE 
67  CONTINUE 

IF(SHP)  65,85,86 
65  JN=JN+1 
86  CONTINUE 

ALoHA=6,2831853772, 

00  12  1=1,36 
00  12  J=l,72 
AC=ALPWA*FLOAT((J-l»*I) 

A(I,J)=SIN(A0»  > . ' ' 

12  89(1, J)=COS(AC) 

AJJ=1.0 

20003  IF(JEEI2,2,3 

2 AJJ=1. 

GO  TO  47 

3 00  209  IOC=l,2 
no  201  1=1, JC 

201  3(I,30)=8(1,30>*AJJ 

914<41)«.975 
B14(42)=i«  000 
B14(43)«1.025 
00  215  tJ=41,43 
IJT=IJ+4 

IJP=IJ4-8  . , 

00  216  1=1, JC 


C8 


- /»  ^ *V- 


51«.(1JT)=35(  jgg, 

''i«^<IJP^=06(36^l 
?15  CONTINUE 
315(i»l  ) = i,  0 
315(i«<,)  = i,  0 
B13(47)  = i,o 

81S(42)=3i<,,4j, 

B15(4,5)=nu((,?, 

B13(«»0)  =nU(43) 

315(46) =P14(42) 

0 i 5 ( 4 9 ) = 0 1 4 { i,  3 ) ♦ ♦ p 
^ If(0(l,2n  6,d.q 

CC (1 , 1 ) =014(45) 

CC (?, 1 ) =014 (46) 

CC(T, 1) =914(47) 

, ^’0  ^0  IB  ’•  2 *’•i‘•15J^7♦3(7,^)••?,/4  , 

CC(1,  l)=Oj^<,J^gj 
CC(2, 1 ) =914(50) 
con,  1)=014(51) 

<=’*<J-i)  h 
C(J,I)=gi5( k»40) 

^10(1M)M 

.r,,  , SlHfC(3,C,C'') 

CONTINUE  ‘^‘'<i»l>-TTT)))/(2,.j 

50  86(361)=-i.o 

JE?  = 0 

"•8  00  49  1=1  , jc 

5 5(I,5)=e(i,3g,4fljj 

continue 

CALL  SU9 

PP11=015(161)*ETmpms 
PP12=EHP1 (IV) /PPii 
TH^=a(»;,2, 

33  Uv^r’ 

36  IV=IV*1 

IP(IVV*IV)  1005(1  91  1 4 
10050  CONTINUE  >2l»^l 

00  10019  IX=l,lv 
FHP(I  X)  =£HPi(lx) 

10049  VEUIX)=VEL1(IX, 

00  40050  1=1, JC 

fJJIJ»2®>=8(I,l4) 

27)88(1,5) 


o 


■’ 


AZZ(I,28)*Cir-8«T,6) 

AZZ<I,29)sCII*BCI,6) 

AZZ(I,30)=8<Iy6) 

AZZ(I,31)*B(I,4I 
AZZd, 32)  =0(1,12) 

AZZ(Z,33)=B(I,13} 

40050  AZZ(I,?4)=AZZ(I,25)*2.0 
00  36  1=1, JC 

36  AX(I) =A2Z(I,34) 

IF(AZZ(!1,25)«NE.0.)  GO  TO  5555 
SLP=(AZ2(9,34)-AZZ(6,34))/(AZZ{9,23)-AZ2(6,25)  ) 

YINT= AZZ( 9, 34) -(AZZ (9, 23) ♦SLP) 

AZZ(10,34)=(SLP*AZZ(10,23) )»YINT 
AZZ(ll,34)=(SLP*AZZ(li,23) )*YINT 
5555  CONTINUE 

DO  40052  <=26,34 
DO  40051  J=i,ll 
B11(J)=B(J,3) 

40051  83(J) =AZZ(J,K) 

DO  40052  1=1,11 
S1=AZ2  (1,23) 

call  OISCOT(S1,51,B11,03,B3,>12O,JC,O,S2) 

40052  AZZ(I,K)=S2 

40056  FORHAT(1PE10.3,1P10F11.3) 

00  5050  r=i,JC 
AZZ(I,2)=8(I,2) 

8(I,20)=B(I,15) 

R(T,21)=9(I,5) 
fUI=ATAK(0(I,21)  ) 

3(T,22)s8(I,38) 

8(I,23)=B(I,6) 

8(1,24) =AZZ(I,24) 

8{I,25)=AZZ(I,25) 

5050  B{I,24)  = {a(I,24)*3(I,3))/C5r.  295S»C03 (9II) ) 

8T(11)=8(11,5) 

00  5001  1=1,10 
5001  31(1) =AZZ(I,19) 

DO  50038  <=20,24 
00  50020  J=l,ll 
811(J)=8(J,3) 

50020  P3(J)=B(J,<) 

DC  50030  1=1,10 
5:1=91  (I) 

CALLOrSUOT(Sl,Sl,Bll,B3,B3,>120, JC,0,S2) 

50030  9(T,<)=S2 

DO  50040  <=20,24 
00  50040  1=1,10 
50040  AZZ(I,K).=  B(I,<) 

CALL  STRFSS(AZZ,AREA, XBAR, YBAR, AYEXO, AYEYO, EMXO,EHYO , EMTB, EHQB ,STR 
INAX) 

C THE  FOLLOHING  STATEMENTS  HAVE  BEEN  ADDED  IN  ORDER  TO  SEND  THE  R 

C VALUES  TO  SUBROUTINE.. HEIGHT,  GHORO» THICKNESS,  AND  CAMBER  ARE  I 

C PITCH  AND  SKEHR  ARE  IN  RADIANS, 

DO  999  1=1, JC 

VS=VK»1.6d78 

X(t)=3(I,3) 
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IF(X(I).NE.O.n  GO  TO  60G1 
SIGMAsP(I,l9) 

GO  TO  600C 

6801  IFCXCn.LE.O.eS)  GO  TO  6000 
IF(X(n.GE.O.  7S)  GO  TO  6030 
SIGMfl*R(I,19} 

6000  CONTINUE 

CH0R0CI)»BCI«6»*DI*M 

thicknsiii^azzcns^moum/?. 

PITCH(I)»ATAN(R(I,5)> 

9BJ« I) =ATAN(BI I.ftH 
IF(CLl(I»)  ?A,?7.2A 
?7  F)f?n  = C.O 
GO  TO  29 
21  CONTINUE 

PXU)::!.  0/(1.0  MC  .2A32*TAN(PITCHm-9BJ(I))/CU  (I  ))» 

29  CONTINUE 

PtTAI{H=PITCHJI)  •57.295R 
SKEHR( II =A77( I »38>/57.295« 

A2Z(I,34l=AX(H 

AVm=<9(I,4)»9(I,l?n**? 

PV  (11=  (9(1,41/9(1, 81 -o<I, nil**? 

VSUBOSQ(  D=VS**2*  IAV(I)»flV(I>) 

VSUBP  (I)  =SQRT  (VSUBPSOdn 

PLFT(I)  =.S*RH0*VSU9P9Q(I)*CL1 (I )•  » ( I , f I *0 1 A*l 4, 5939 

999  CAM9-R (I )=.0679*9 (I,18)*0IAH 
PPS=ePM/6C, 

00  996  1=1, j: 

PXT8I(I)=PI*9(I,3)*9(I,5) 

PXTB(I)=PI*R(I,3)*n(I,8» 

996  CONTINUE 

EA1=(EAR*3. 14159»0IA**2» /4.0 
AL=3.14159»0,7*B(7,9) 

AP=EAl*(  1.067-0. 229*AL  J 
VA=VS*R(7,4) 

VR=SQPT (VA**2* (0.7*3.14159*RPS*DIA)**2) 

TC=2.3’9(6,2)/ (RHO*AP*VR»*?J 
SIGMA7=(64.4»HEAD) /VR**2 
IF(JPR.NE.2I  GO  TO  40D 
10C42  FORMAT  (IH  I 

HRITE(6,77)  PPl , PP3, PP2 , PP5 

77  FORHAT(4X,»CPTI=*,1PE10.4,4X,*CPSI=*, IPE10.4,4X, *CTSI=*,IPE10,4, 

2 4X  ,»CTSI/CPSI=*,1PE10.4I 

WRITE (6, 10042) 

WRITE (6,78) 

78  F0PMAT(6X,»  X* , 8X , »TANBI*,7X, *T AN  B* , 9X , *G* , 9X, • UT/2V*,7X, 

2 *UA/2V*,7X,*0CTSI*,7X,»0CPSI*,6X,*VR(M/SEC)*,5X,*CAVV*) 

no  79  1*1,11 

VRSI  = VSU8R(I)  *Elf 

79  WRITE (6,6010)  B ( 1 ,3) ,B (1 , 5 1 ,B f 1 , 8 ) ,B ( I ,l M ,B (I ,l 3 ) ,B ( 1 ,121 ,B(I,15) 

1 ,B(I,17I ,VRSI,B(I,19» 

6010  FORMAT (1P10E12. 4) 

WRITF(6, 10042) 

WRITE(6,6011)  PP6,PP8,PP7,PP10 

6011  FORMAT  (5X»CPT**1PE10. 4, 5X*CPS*»1P£10.4*5)(*CTS**1PE10.4,6X, 

2 •CTS/CPS**1PE10.4) 


HRITE<6,100*i2) 

WRITE  C6, 76)' 

76  F0RM4TC5X,*  X*, 6X,*CL*,6X,*ALIC0EG)»5X,*FM/C*,7X**C0/CL*,7X, 

2 •F(X>*,5X,*LI{N/M)*,3X,*TETS(0EG)*,2X,*CC/RD)LE»» 

3 3X,*(C/R0)TF»,5X,*T/RD*  > 

DO  6002  Ial,ll 

600  2 WRITE (6,200^6) BCIt3),CC0NE(X) *CCTWO (Zl » CCTHR (I) tCCFOR(Z) ,FX (I ) ,PLF 
1T(I)  »AZZ(Z«38),AZZ(If28) tAZZ(I,29) »AZZCIf 34) 

20046  FORMAT(1PE10.3,1P10E11.3) 

WRITE (6„10042) 

SHR=THR*SWT 

STHR=ATHR*SWT 

WRITE  (6i  10044)  PPU,PP12,  ETHRUS,EWA<E,VK,SHR,VMS,STHR 
WRITE (6,6003) 

600  3 FORK  AT (iHl, 4X,*X*,7X, •AREACHZ )*,4X,*XBARCH)»,2X, ♦TSAR (M)*,5X, 

2 •IX0(M4)*,4X,*IY0(M4)*,4X,*KX0(N-M)*3X,*KY0(N-M)»,3X,*MT8(N-K)*, 

3 ,3X  ,*H08(N-M)»,3X,*MAXSTRESSIPA)») 

400  CONTINUE 

00  122  1=1,7 
IF(I-l)  120,121,120 

121  SX(I)=8(I,3) 

GO  TO  122 

120  SX(I)  =B(I4-1,3) 

122  CONTINUE 

IF(JPR.NE.2)  GO  TO  10091 
00  124  1=1,7 
APR=EL2*AREA(I) 

XPR=EII»XBAR(I) 

YPRsEUl*YBAR(I) 

X0I=cL4*AYEX0(I) 

Y0I=EL4*AYEY0(I) 

X0K=SIt1*EHX0(I) 

Y0N=SIH»EHY0m 

T9H=SIF*EMT8(I) 

QBM=SIM»EM08(I) 

STRK=SHX»STRMAX(I) 

124  WRITE (6,20046)  SX(I) ,APR, XPR, YPR,XOI, YOl,XOM, YOM,TBM,QBM,STRH 
WRITE (6,10042) 

WRITE(6, 10042) 

WRITE (6,6004) 

6004  FORMAT (5X*X*8X»RAKE*5X»P1  XTANBI*3X*PI  XTANB*) 

DO  6005  1=1,11 

6005  HRITE(6,6006)  8 (I ,3) , RAK( I) , PXT8I (I) ,PXTB(I ) 

6006  FOR‘1AT(1FE10.3,1P3E11.3) 

10091  CONTINUE 

CTS(IK)=PP7 

CP9(IK)=PPe 

CPT(IK)=PP6 

CALL  WEIGHT(JC,SIGMA7,IPO,IPR,HU30IM,IK,HUB) 
HU9SPAC=2.0*PI»X(l)*OIA*6.0*SIN(PITCH(l))/O(l,4,IK) 
3LASPAC*HU3SPAC-AX(1)«DZA*6.0 
FILSPAC=HU8SPAC-AX(1)»DIA*6.0*1..9 
IF(JPR.NE.2)  GO  TO  6902 
WRITE (6,998)  TC,SIGHA7 
SBL=Bl.ASPAC/DIA/12. 

SFI=FILSPAC/DIA/12* 

WRITE(6,994)  SBL,SFI 
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«8  FORHAT(/20X,*8URRILl.  THRUST  CQFFF  TC**»E10.  V/20X,*3URR 
ITATION  COEFF  SIGMA ( 0. Tl  i £ 1 0. 41 
994  FORMAT  {/20X, ♦CLEARANCE  AT  HUB  BETHFEH  BLADE 

2 ?0X, ♦clearance  at  hub  6FTKEFN  FILLETS/3** ,F 13, 8 » 

690?  CONTINUE 

IFtJPR.EQ.l)  JP5=2 
00  64  1=1,11 
AZ7(I ,23)=AZZ(I, 1?) 

ABCd,  1)  = 3(I,3> 

A3n(I,2»=0(i,14> 

A-3C(I,3)  = AZZ{I,27» 

A9C(I,4{  = aZZ(I,2‘'> 

A3C(I,«»=AZZ(i,29) 

A3C(I,6)=ArZ(l,31S 
AlCd  ,7»=3(I,12S 
64  ABCd  ,8>=aZZd,34) 

Ai:(2,c)=9{9,2) 

IFdOrK.NE.3)  GO  TO  6? 

IFlK.NE.l)  GO  TO  f? 

IFCNK.EO.  1)  GO  T«  63 
GO  TO  62 

b’  IF(NIT.EQ.l)  Go  TO  M 
GO  TO  62 
61  9Z(1)=1.0 
3Z(2) ="1A 
3zn)=wFAr) 

3Z(4)  =y|*M(l)/6.].3 
azc"?)  =f'.931 
az(6) =THR*2.1 
BZ(71 =VEL(l)*1.6a76 
Bzcn  ==Ho 
3Z(9)=EMA<F 
3Z(1C) =9LA (15 
BZdl ) =9Z  dO) 

0Z(15) =CO 
3Z(16) =C. C 
BZ(17) =0. 0 
9Z(18»=0.0 
no  300  1=1,12 

300  9Zd»99»=COMd,l) 

I*J1  = 6 

DO  301  IN=18,90,18 
IN1=IN1*1 

IFdM.FO.90)  IN1=13 
no  301  1=1,9 
IFd-l)  302,303,302 

303  J=I 

GO  TO  306 

302  IF(I-9)  305,304,335 

304  J=I»2 

GO  TO  306 

305  J=I4-1 

306  3ZdN  + n=0(J,  IN1,1» 

301  CONTINUE 

on  307  lM=27,81,ie 

IN9=IN*9 

DO  307  1=1,9 
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BZ(IN+l)=BZtIN9*I> 

CONTINUE 
BZ(10)sO(l,**»l) 

8Z(J1)*QC1»‘*»2> 

00  308  I=cl,9 
IN=I»72  . 

IF(I-1>  309»3l0f3Q9 
J=I 

GO  TO  311 

IF<I-9)  3i2»313t312 
J=I^2 
GO  TO  311 
J=I>1 

0Z(IN)  =O(J»10j2) 

CONTINUE 
00  5000  J=l»9 

if(J.gt.i)  go  to 

B2{J«-45)=GA11(  J) 

GO  TO  5100 
IF(J.EC.9)  GO  TO 
8Z(J»45J sGAll (J»l> 

GO  TO  5000 

5002  8Z{J4^45I=6A11(J>2) 

5000  CONTINUE 

CALL  OLO(BZ|IPOI 
CONTINUE 
n(l0»2f 2)»8Z( 16) 

FA9=0  C10»2>2)/0(10»2»1) 

AFRsO (10»2|1) 70(10 >2»2) 

00  320  I=1»JT 
OEXd  »3,l)sFAR*0(I  »7»2) 

OEX(I»3t2)=AFR»0(I »7» 1) 

CONTINUE 
IC=IC*1 

IF(IC)  17,17,20 

CALL^FlLo(A8C,UA,UT,UR,IK,OEX,COf1,I‘’0) 

IFCIC)  18,18,20 
00  19  I*l»ll 
UAl<I,IK)sUA(I) 

UTia,lX)-UT(I) 

UR1(I,IX)-UR(I) 
no  326  1=1,11 

IFC0EXa,3,U-0(l,7,l))  325,325,20 
UAKI, 11=0.0 
UTKI, 11  = 0.0 
UR1CI,1)  = 0»0 
CONTINUE 
CONTINUE 
UAKI  ,1)  = 0.0 
UTKI, 11  = 0.0 
URKI,1»*0.0 

IF(JPR.NE.3)  GO  TO  46 
IF(IC)46, 95,46 

m FOWT(lHl////50X,MNOUCeO  V|C0CmES*///33*,*tC.13XVU»FVS* 

1UR/VS*,11X,*UT/TS*/22X,*0N  AFT*) 

■* 
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00  107  1*1,11 

WRITE (6,1 Oo)  OtX (1,3, is ,UA1 (I, it ,UR1( 1,1) ,UT 1(1,1) 
iO«  F0RWAT(23X,4F16.«*» 

107  nOSTINUF 

WRITE (6,109) 

109  FORMAT (///2?X,*ON  FW)*) 

DO  110  1=1,11 

WRITE  (f>,108)  OFX  (I  ,(,?)  ,UA1(I  ,?)  ,UR1(I,2)  ,UT1  (1,2) 

110  COSTINDE 
9097  rONTiNUF 

46  IFCjPS.fQ.O)  JPi?*! 

C'O  111  1=1,11 
UA(T)  ='IA1  (1,10 
UR(I) =U«1(I,IK) 

111  UT(I) =UTi (I,I<) 

IF(«:hF)  .-^1,81,69 

31  IXA=IVA*1 

IF(IVV-IVA)  96,70,70 

70  GO  TO  10046 

69  IF(NTT.EO.NXVVfl ) GO  TO  96 
IF(nT  .‘.’c  IJU-1)  GO  TO  97 
9‘»  NF  = NF*1 

IF(IK-I)  93,93,94 
9’  FW9SHP  (NF)::PP17 
GO  ”0  96 
94  NF=NF-1 

AFTShPCNF) =PP12 
GO  TO  ^6 

97  IF(NT,"0.ir0)  GO  TP  98 
P6  JM=l 

IF(JO)  71,72,71 

79  Xi(r9 

GO  TO  73 

71  K=1 

’7  CONTINUE 

TF(JQ)  74,75,74 
79  J0=1 

GO  TO  <'0 
74  JQ=0 
8C  CONTINUE 
103^2  CONTINUE 

IF(SHP.F'3.C.0)  go  to  10»)73 
IF(NTT.NF.NXVV)  go  to  10073 
00  99  I=1,NXVV 

811  (1)  = (FWuSHP<I)7AFTShP(I).)/-2.,0 
VEL(I)=C(I,1,I<) 

99  B3(I>=VEl(I) 

S1=SHP 

call  OISCOT  <Sl,«'l,811,P3,e3,-l20,NXVV,0,S2) 

V6LK  1)=S2 
00  154  I=1,NXVV 
811(I)=WcL(I) 

EH*' a ) =0(1, 2,  IX) 

154  93(I)=EHP(I) 

S1=VELK1> 

CALL  OISCOT  (Si,5;i,311,83,t)3,-i20,MXVV,0,S2) 
EHP(1)=S2  < , ,v 
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VELdJsVELiCi) 

IVV=1  . . 

IK=l 

GO  TO  ifS 
10073  CONTINUE 

OR-(OC1I)2»2}/O(10y2>in«»2 
CT?T=CTSa)»0R*CtsC2l 
CPST*CPS(1)  ♦0R*CPS(2) 

CPTTaCPT(l) ♦0R*CPT(2» 

ETAPsCPTT/CPST 
ETA0*0(l,3,ll*CTST/CPST 
HRITE(6,992)  CTST,CPST,ETAO 

992  FORMAT (IHlflOX,* PERFORMANCE  OF  SET  OF  CONTRAROTATING  PROPELLERS*// 
1//2X,*CTS=  *,F«.5//2Xt*CPS=  *»F8.5//2X,*ETAD*  *»F».5» 

STOP 

ENO 


f 
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SU3R0tTINF  SUB 

OIHENSTON  CHORO(38) tTHICKNS(36) ,CAHB?R(38) ,P1TCH(38) »S<FHR(38» 
1,XI(S8) 

C0MH0N/CW£IGHT/Xl,CM0lf0,THlC<Ni,C*HBEP,PlICH,5KfHR,0IAM,22|DEM 
1 ,-?AKF,SI,PI,PPy,PPtt|PP9,PPil,»'MA((F,rfS,RPS»SIGMA,r8R,8T 
OI1FNSICN  9(38*38)  *dl(181)  »B?(161)  *B3(  16 1)  > d <4  ( 362)  * 85  ( 362) « 3<>  ( 362) 
1 *12(181)  ,88(181)  ,B':i  (181)  *910(181)  *911  (181)  *812(151)  * S13  ( 181)  , Bl6( 
2181) *815(161) *AZ (38*38) *GH (38*38) ,C (38*38) *CC (38, 38)  * INDEX (38 *3) 
3,5(38,72)*OB(38*2?) 

COM-^ON  8*91  *3?*03*n4*P5,66*P2  *«38,"B*310  *911*  91?»B13*Bl4*ei5*AZ*3H* 
1C*CC*INO:X*A*09 
rONHON  IiJ*J9*JC*  JO,JOn*JFF 
CO-ION  CLKll) 

CO*(MON  CCONF(  11)  *CCTWO(ll)  THR  { 1 1)  , CCFUR ( 1 1 ) 
rOMHON  P^l*F??,Pf-  3*o=4*FP5  *P5'fa,POio 
10052  m 10)25  N=1,JC 
DO  20021  1=1, JC 
A4'‘.=  1,/9{I,‘5) 

A4H=8 (M,»)/n(I,3)*AAC 
AftOrATANCOd.Sn 
IF (AAh-AAG) 10019,10018*10015 
1001?  9?(I)  =COS(AAO) 

■•3(1)  =<;iN(aAO) 

GO  TO  20021 
10)1=  S = l.*Aft»-*«r 
T = SOtiT(S) 

’/  = 1.  f A?.G**2 
W = 'iO^T  (V) 
ac=T-w 

')  = =XP  (Ai.) 

R=(  ( C-1.  ) / AAH*(Aflr./(W-l.)  ))*U)  ♦*3(9,2) 

flC=1.5 

AT  = .2  5 

y= (1. /(2.*9(9, 2) *AAG) ) ♦ ( (V/3) ••AO) 

Y = ( (9,  • AAG**2) +2. ) / (V**au  ♦U3.*AAH  •♦  2-2 , ) / ( S** AC) ) 

Z=l./  (24.*C'(9,  •»)  )*Y 
IF (AAH-AAG) 10021,10021 ,lu020 

no?:  AF=i.  ♦i./(r-i.) 

AA=y* (1./(F-1,)-Z*AL0G(AF) ) 

n?(I>  =2.*9(9,2)**?*AAG*AAH*( I.-AAG/AAH) *AA 

R3(I)  =3(9*2)*(1.-AAG/AAH)*(1.»2.*B(9,2)*AAG  *AA) 

GO  TO  20021 

10021  AG=l» fl./(l./R-l  .) 

Al=-X* (l./(l,/9-l. ) »7*AL0GCAGI) 

•J?(I)  =^(9*2)*AAG*  (1.-AAH/AAG)*(  1. -2**9 (9*2) ♦AAG*A3) 

93(1)  =2.*0(9,?)**2*AA6*(l.-AAG/AAH)*AB 

20021  COnFINUF 
2002A  F0’HAT(9F12.4) 

IN=0 

DO  1 1=1,181*5 
IN=INfl 

1 Bl(IN)  = .5*(l.»B(l*3))-.5*(l.-Ba,3))*C0S(<FL0AT(X-i)/57.2957)  ) 

DO  2 I=1*JC 

311(I)=B(Z*3) 

B12(I)=B2m 

2 8l3(I)=a3(I) 
no  3 1=1*37 
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si=Bia» 

call  OISCOT(SifSi»Bll,Bi3»BlS.»»120',JC,Ot$3) 
B15CI)=S3 
00  5 1=1,37 
S1=B1(T>- 

CALL  0ISCOT(Sl,Si,Bll,B12,bi2,-lZ0,JC,0,S2} 
5 Bi4(I)=S2 

00  k 1=1,37 
82(11=814(1) 

4 B3(I)=815(I) 

no  10022  L=l,35 
N1=37*L 
N2=37-L 
82(N1)=82(N2) 

10022  83(Nl)=a3(N2) 

C2=2./72. 

NP=T2 

NH=36 
XNP=NP 
S=0.0 
SL=0. 0 

00  20  T=1,NP 
S=S*B2(I) 

20  SL=SL»e3(I) 

84(1)=S/XNP 

85(1) sCL/XNP 

00  40  1=1, NH 

S = 3.0 

5L=0.0 

S1=0. 0 

SL 1=0.0 

00  30  J=1,NP 

S=S»02(J)*B8CI,J) 

SL=5l+03(J)*BB(I,JI 
Sl=Sl*a2(J)*A(I, J) 

30  SL1=SL1*B3(J)*A(I,J) 

11=5* !♦! 

34(I1)=S*C2  w t 

85(I1)=SL*C2  ..  . 

B4ai*l)  = Sl*C2 
40  85(Ilfl)=SLl*C2 
3(l,9)=a4(l) 

8(l,10)*85(i)  , , • , ■ 

JCM1=JC-1 

DO  10023  LK=1,35 

L=L<^1  , . 

K=5*LK 

R(L,9)=f4(K^il 

10023  8(L, 10)=05(Ktl) 

9(37, 9)*B4(181I 
8(37, i0)=85(18i) 

CPm*((l.f0(l,3n-2.*8<N,3))/(l«-8(l,3)) 
IF(ABS(CPHI)-l.)  20051, 20051, ,2»0050,  , , 

2D05C  CPHI»1,  V . • , 

20051  B(N,ll)»ACOS(CPHI) 

8(1, 11)=. 0 

B(JC, 11)*3.1415927 
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C0S1* 3. 1415927 
no  10025  1*1, JC 
SM«»=SINCPlO*T(I)  *3CN,  IIH 
CHP=C OS C FLOAT (I»*9(N,lll) 

-IF  JM-  1)  10027,  100  ?6, 10  027 
10027  IF(M- Jf*)  10028, 10r?9, 100  «e 
10026  AZ»1=.0 
n?N=. 0 
N?=If  I 

no  20  026  K=1,N2 

IF(<- JO  10  070, 10  07  0,2  0026 
100^''  AZN=AZK»CON3*FLOAT  (I)*Q<K,9) 
oZN=nZN*i;ON3*FLOAT  (I)  ♦9<K,  lO) 

20026  CONTIN'IE 
AZL=. 0 
BZL=. 0 

IF(N2-JC)10Q6  0,1':0'0,  10030 

10063  N1=N2»1 

nn  20  336  rt=Nl,J''. 

aZL  = A ZL»FL0AT (L I ♦ ^ (M, 91*COk3 
200  ’6  = 9 7L*FlOAT  (L)  *9  10)*CON3 

GO  TO  1 C3TC 
100  2S  AZ'^=.0 
0 7‘J=.  3 
M2  = I + I 

on  20  3'’9  K=1,N2 

C<P=COS  (FLOAT  (i<-l)*n(N,il)) 

IF(<- J«-l  10C71,  1007  1,20 '»2P 
10  07  1 AZN=AZK-CON3*CHP*FLOAT  9)  ♦OKP 

37N=nZN-CON3*CHP*FLOAT(Il*b(<,10) •CK^ 
20020  CONTInUF 
A7l  =,  0 
9Zu=. 0 

IF  (M2- JG) 10  061,100  30,  1 0030 
10061  Nl=N2fl 

on  20  3 39  f1  = Nl  , JC 
L = *(-l 

C<==COS  ( FLOAT  (L)  *9(M,  ID) 
AZL=AZL-C0N3*CMP*FLOAT (L) * 0(1 , 9 > *CKP 
20039  eZL  = 9ZL-C0N3*CMP*FL0AT (t ) *6 (M , 1 0) ♦CK® 
GO  TO  10030 
10026  AZN=, 0 
nzN=. 0 

CONI*  3. 1415927/SIN  (QCM,U)) 

N2*I+1 

00  20328  K=1,N2 

C<P=COS (FLOAT (K*1)*8(N, 11) ) 

IF(K- JC) 10072, 10072,20026 
10072  AZN*AZNfC0Nl*SMP*9 (K, 9) *LKP 

BZN=3ZN»CONl*SMP*0(K,lO)*CKP 

20028  CONTINUE 
AZL=.0 

9ZLS.0  ..  •’  ■ 

IF (N2-JC) 10 0&2, 10030,10030 

10062  Nl*N2*l 

00  20Q38  H*N1,JC 

79 


SKO=SIN(FL0AT(L)  •8(N,  ID) 

A2L=AZL*C0Nl*CMP*8CM,q)*SKP 
200:^8  BZL=BZL*CONl^CMP*aCM,10)*SKP 
10010  AZ(I, N)=AZN*AZL 
8H(29N)x8ZN+BZL 
10025  CONTINUE 

IF (9  { i,  1)  IZOOi^a,  20043, 10034 
20043  CONTINUE 

00  10031  1=1, JC 

CC(I,l)  = (l.-B(l,3n*(B(I,5)/8(I,8)-l.)*B(I,4) 
no  10031  j=i,jc 

10  0 31  C(l,J)=FL0AT(J)*(A2CJ,I)*Ba,5)*8H(J,I>) 

910(181)=JC 

CALL  SIHEQ(JC,C,CC) 

10034  00  10015  1=1, JC 
8(1,12)=. 0 
0(1, 13)=. 0 
8(1, 14)=. 0 

no  10035J=1,JC 

IF(9(1, 1) ) 10036, 10036,10037 

10016  Q(I,12)=3(I,12)*-FL0AT(J)*CC{J,1)*AZ(J,I)/8(I,4)*(1./C1.»B(1,3))> 
8(1, 13) =8 (I, 13) AFLOAT ( J) *CC ( J , 1) «3H ( J , 1) /B(I ,4) • ( 1./ (1 .-B ( 1 ,3) ) ) 
3(I,14)=Ct(J,l)*SINCFL0AT(J)*B(I,ll>)/B(I,4)*B(I,l4| 

GO  TO  10035 

10037  B(I,12)*B(I,12)«-FL0AT(J)*B6(J»*AZ(J,I)/B(I,4)*(1./(1.-BJ1,3>)  ) 
8(I,13)=iJ(I,13H>FL0AT(J)*B6(J)  *8H  ( J,  I) /8  ( 1 ,4)  ♦ (1 ./( 1 .-8 (1 , 3) ) ) 

8(1,  14)=B(i,14)4-B6(J)*SIN(FL0AT(J)*B(I,ll))/8(I,4) 

10035  CONTINUE 

. , B(JC,.14)  = .0 

20001  00  10036  1=1, JC 

3(1,  15)  = (8(I,14)»9(I,4)*(B(I,4)/8{I,3)-B(I,13)*8(I,4) ) )*4.»3(9, 

12) 

8(I,16)=3(I,15)*9(1,4I 

0(1,17)  = (B(  1,4) /6(  1,3)  »B (1,14)  *8  ( 1, 4)  * (B ( 1, 4’)*6  ( 1, 12)  • B (1 , 4) ) ) *4. 
1*0(9, 2) 

3TT=ATAN«B(I,8)) 

bti=atan(B(i,5))  : ’ 

0<I,13)=?.*3.1415927*B(I,i4»  ^COS  C DTI)  / ( 1 . /B  (1 , 8) -8  (1 , 13) ). 

3 (I,  19)  =64. 31»  (8  (4, 2) -8 a,  3)  <^8  (3,  2)  /2.)  * (SIN(8TT)/(8(I,4)  *8  (7,2) 

1*C0S(BTI-9TT))>**2 
IF(I-1)9,9,6 
6 IF(I-JC)10,9,9 

9 3(1, 20)=. 0 
0(I,21)».O 
GO  TO  11 

10  CONTINUE  ‘ ’ 

8(I,2  0)  = (l.“BlI,n*8(I,6)/BMil8)'*B(I,5))*B(I,15) 
B(I,21)=(1.*8(I,7)*B(I,6)/B(I,18)/B<I,5))*0(I,17) 

11  CONTINUE 
B(I,22)=B(I,20)*0(I,4) 

01(I)=B(I,3)  * 

B2(I)=3(I,15) 

B3(I)=n(I,16) 

B4(I)=B(I,17) 

85(I)=B(I,18) 

87(I)=0(I,19) 


Copy  avcdlcwls  io  DDG  does  noi 
pciEiil  fa3y  leyiblo  repioducUotf 


RSCn  s9(l,?0) 

8<}<i»  =aci,2n 
Bioa)>ecZ(22) 
lOO’lS  CONTINUr 

PPl*SIHPUN(Bl,83,JC) 

PO?=SIMPUN(Bl|B?iiJC» 

P»^=aIMPUN(31t8%tJC) 

PP4=PP1/PP^ 

PPS=PP2/PP3 

o®6=SI^PUN(Bl,aiO, JC» 

P07=SIHPUN(dl»99,JC) 

poa=SIHPUN(81,?9»jn 

PP9=PPf./PP6 

PPn  = PF7/PC9 

no  1C0?9  1=1, JC 

00  100«*8  L = 1,JC 
vn=3(L,3)-aa,3» 

IF(VO)  860,«fi0,«71 
5«>C  XO  = O.T 

GO  TO 

3T1  rONTINUt 

XO-9(L,3) -P(I,3) 
if)l  CO'^TpuiP 

9’(L)  =xO*-3(L,2*‘) 
ion:  ? J(L)  =xO/9(L,3)*P(L,21) 

I'^(JCI-2)  1C091, 10041,10059 
10041  1(1, 25)  = . 0 
3(1,25)  = . O' 

3(1,27)=.  0 

n(I,?0)=.0 

GO  TO  lQf’39 

l"05q  1(1,25) =bIMPUN(91,B?,jC  ) 
1)**?)/(16.*6(9,2)) 

9(1,25) =SIHPUN(B1 ,33, JC  1 

l(5,2)*n(q,2)) 

JTI=ATfiH(9(I,5)) 

b9T=*;iN(9TI) 

C3I=COb(9TI) 

9(I,27)=P(I,25)*CBIf3(I,?6)*S0I 
Bd,?*!)  =9(i,25)*Sl3I-3(I  »26)*Ci)I 
no  59  CONTINDc 

DO  ?Q5  I--1.JC 

9(1,  12)=‘’(1,4)*9(I,12) 

9(1,1  3)=u(I,4>*9(I,13) 

206  0(T,14)=B(I,4)*9(I,14l 

20090  IF(JEF)200fl2, 20082, 10081 
2008?  96(35l)=-1.0 
20081  rONTINUt 

IF(9&n61i)  10  08.0,  lOOM, ,10001 
10090  CONTINUF 

915(181>=PP10 
00  10049  I»1,JC 
IF(9(I,6))  702,702,203 

702  CC1=0.0 
GO  TO  704 

703  CCi=8{I,18)/B(I»6) 


tl 


• (9(9, 2) *9(3,2) ••3* 3. 1415927* 3(7,2 

* 5(8,2)*9(3,2)**2^B(7,2)**3/(16.*3 


'«  » 


• • 


70%  CONTINUE 

CL1CI>=CC1 
CC2*1.5%*CC1 
CCT=.0679*CC1 
IF(CC1}  701t700»701 

700  CC4:Q.O 
B(I|38)«CC4 
60  TO  52 

701  CC%*BCIf7)/CCl 
8(I,38)^CC% 

52  CCONE(II>CCi 
CCTMO(IMCC2 
CCTHR(l)sCC3 
10045  CCF0R(II=CC4 
10081  CONTINUE 

R5C362)=PP7 
B6(362I<PP6 
20041  RETURN 
ENO 


SUBROUTINE  S1HF0( JC»C|CC) 

OIRfNSION  BH(36»3S)  ( 38,3S} , BIS  (1ft)  ,C(39|34)  tCCdftfSft) 

BPlaJC 
HPC=HPl*l 
MPHsHPl 
00  SO  1=1, HPl 
60  C(!,»-PC)=-CCa,l) 

00  SI  1=1, HPl 
DO  SI  J=1,KPC 
91  3H(I,J)=0.0 

00  ’0  T=l,HPl 

21  9h(I,  l)=-C(l,Id) /cU,l> 

00  21  1 = 1, PPl 
J = U1 

21  BHd.JJ^l.O 

< = ? 

93  92, 9?, ^4 

9’  CONTINUE 

PO  9f)  i = l,HPP 
f>0  “6  l = l,HP<~ 

6b  A(’',J)=  C (K,  J)  *3HtI , J) 

jO  31  1=1, MPO 
AN=3. 0 

37  J=1,«PC 
d7  ftN=AN+4(I,J) 

59  019(I)=4U 

HOJsMpn.j 

JO  S9  1=1, MPO 
00  o9  J=1,'"PC 

o9  OCd,  J)=-915(I*1)/P1S(1  )*3N(l  , J)  fOH  (I»l,  J) 

no  90  !=1,“P0 
00  9d  J = 1 ,**PC 
90  9H(T, J)=CC<I,J) 

< = <♦1 
00  ""O  P3 

94  OONTIN'): 

00  99  J=1,“PC 
99  CJ(J,1)=CC(1,JI 

1 F0PM4T  (lPf.El^.4) 

2 FOPHAT (6F8.4» 

HETURN 

FN"> 


Copy  o^dlable  to  DDC  do<M  not 
pen^l  fully  legible  loproducUon 


J 


SUBROUTINE  OXSCOT  CXA,ZA,rABXtTA8Y>TABZ>NC>Ny9NZ>ANS» 
OIKENSIQN  TABX(l),TABTCl},TABZ(i) »NPX(37) ,NPV(37) ,YT(3Zi 
CALL  UNS  (NC»tA,lOX,IOZ»XHS) 

IF  CNZ-ll  5f5fl0 
5 CALL  OISSER  (XA»TABX»1»NY»ZDX»NN) 

NNN*IOX*l 

CALL  LACRAN  (XA, TABX(NN) >TABY(NN) >NNNf ANS) 

GOTO  70 
ID  ZA»G»ZA 
IP1X=I0X4-1 
IP1Z=I0Z*1 
IF  (lA)  15«25ylS 
15  IF  (2ARC-TA8ZINZ))  25,25,20 
20  ZARG3TABZ(NZ) 

25  CALL  OISSER  (ZARG,TA8Z,1,NZ,I0Z,NPZ) 

NX=NY/NZ 

NPZL*NPZ^IOZ 

1=1 

IF  (IHS)  30,30,40  t .».>  - 

30  CALL  OISSER  (XA,TABX,i,NX, lOXyNPX)  ^ - 

00  35  JJ=NPZ,NPZL 
NPY{I)s(JJ-l)*NX+NFX(l» 

NPXdJ -WPX(l) 

35  I=I*1 
GOTO  50 

40  00  45  JJ=NPZ,NPZL 
:S=(JJ-1J*NX*1 

CALL  OISSER  (XA,TAeX,IS,NX,IDX,NPX(I)> 

NPY{I)=NPX(I) 

45  I=I*1 

50  00  55  I=1,IP1Z 
NLCC=NFX(I) 

NLOCYxNPV(I) 

55  CALL  LAGRAN  (XA,TA BX(NLOC» ,TA BY tNLOCY) ,IP1X, YY (II » 

CALL  LAGRAN  (ZARG, TABZ(NPZ) ,YY,ZP1Z,ANS) 

70  RETURN 
ENO 


■ i-i-' 

:84': 


SU-l^OUriNE  UNS  (ICtlAfXSXf  X02»XHS) 
IF  tl'il  5,5fld 
•:  THS=1 
NC=-i'r 
GOTO  15 
10  iH$=a 
NC=IC 

15  IF  (NC-lOO)  ?0|25,?5 

?D  ia=o 

GOTO  30 
?*?  ia=i 

NC=SC-100 
TO  IOX=NC/lC 
ID7=HC-IOX*10 
RPTi|«?N 
£NT 


SUTROUTlHs  LAG6AH  (XA« X»T»N,ANS) 
OlfFNSICN  xm  ,r  (II 

PO  3 1=1, N 

o^OO  = Y(  I) 

"0  ? J=1,N 
A=X(I)-«(J) 

IF  (A)  1,2,1 

1 0=(XA-Y(jn/A 
P^00=P90L*2 
? fontimi.ic 
3 SUH=SUH*FROO 
AN3=SU« 

»FTIJRN 

'■NO 


Copy  'ccrc^ablo  {o  DDC  does  nol 
penait  fully  Icyiblo  reproduction 


SUa<tOUTXNE'  OISSCK  fXAyTiltf Zt NXylBfNN): 
JilMEMSION  TA8C1I  r . . • 

N^T*IO»l 
NPB*NPT/2 
MPU*MPT-HP0 
IF  (NX-NPT)  tOySylO 
5 NPX=I 


RETURN 

10  NL0H*l4>NPB 

NUPP«I4-HX-«HPU*1J 
. ^ DO  15..I.I*KL0Wi.HUPP.  . . , 

, NLnC*lJ  • . ' • ; 

IF  (TABfXI)-X*)  15t20  *20 
15  CONTINUE 

NPX»MUPP-HP0»1 

RETURN 

20  NL*NLOC-NPB 
NUaHL+IO 
DO  25  JJ>NUtNU 

NOIS*JJ  . .... 

IF  (TAB(JJ)-TA.B(JJ*in,  25,SB»25 
25  CONTINUE 


HPXxNL 

R«^TURN 

3C  IF  (TAPiCN01S)-XA> 
35  N^XsNOlS-IO- 
- peT,uR;N  \ 

AO  NPX*NOIS+l  .j 
RETURN 
ENO 


A8t  i9fZ5 


FUNCriOil  S1NPUN(X,Y|N| 

C FO^TRAK  iV  FUNCTION  FOR  SlH^SONS  RM.?  INTECRArZON 
C AO:)ITRAkY  no.  and  LENGTH  INTERVALS  K.  MEALS  NSROC  COOt  8<t2  10-5'’67 

OIMENSTCN  X(  ?),Y(  21 
IF{N-2)  7,  5,4 

5 S*(r<l)»Y(2l)»(X<2l-XliI»/?, 

60  TO  f. 

4 MaN-l 

d n,10,lt 

11  Nsm-2 
60  TO  * 

9 S=(X(2)-Xin)/6.*tY(l)*(S*-(Xt2)«‘X(l)»/(X(3l-X(l»l)«>r(2»*t3.»{XC2> 
l-X(l)  l/(X(3)-X(2)>  )*Yt3>*CC{Xt2J-X(in**2)/t  (X  (X)  «X(  I »)  * TXU)  -X{2> 
2)»») 

L=X 

63  TO  12 
1"  $=0. 

L='» 

12  M=N-1 

30  1 K=L,‘^^^ 

IF(43S(*(i<-l>-X(in.6E.A0S(X(K»rX(l»)>  GO  TO  3 

iF(ae6(x(<»i»-x(in.GT.AesU(K)*x{t»n  Go  to  t 

T write  (A, 21  K,  XC<) 

2 F0«’‘<AT{23mCN0N  MONOTONF  X SiHPUN  I4,1PE12*4> 

7 S-1, 

60  “"n  . 

1 S=SM  X(K^1>-X(K-1  n/6.*<y<K-l»*  (3.**  (X(<tl)«X  <K-inV;iX{K)4X  (<•!)))♦ 
1 <Y«I  ♦ (!.♦  (X{K*l}-y(K-l)l/{XlK)-X<K-l))  ♦<XU)*X<,K-l)-»  / <Xt<tl>-X<<) 
1 »»  MY  (Kfl»*(2.-(X(K)*X<K-1>»/(X  (<fl)“XlK)  ))  ) ) ) 

e ST^PUNsS 
RET'IPN 
£N0 


Copy  otyailable^lo  DDC 

,p^t  My  Icgilsla  xeproducUoa  ^ 


SUBROUTINE  STRESS(«ZZ* AREAfX8ftR»T9AR«AYEX0f AyEYOr€HXO»CHTO,eMTt*€H 
IQBySTRHAX) 

OTHENSXCN  CHOROtSOI >TKZCKNS(S8>f CANS£R(38> yPITCHiSS) »SKEMR(3S) 
1»XI(58} >BT(li} 

CONHON/CMEIGHT/XItCHORO»THlC<NS>CAHBER»PITCH,SKEWRyOIAH»ZZtOCH 
i ,RAKEfSI,PI,PP7,PP8,PPq,PRlljEWAKE,YS,RPSfSISHA,EfiR,BT 

OIKENSION  AZZf38)33) 

OINENSIOK  AREA(7»;X9AR(7)»YeAR(ri»AYEX0(7)>AVEY0{r»t£HX0(7),CHY0(7 
a feHTB{7)fE?lQ0C7»,STRHAXf7J 

IILL  ppograr  converting  has  done  by  bob  hccalley 

JOHN  METZ  f AT  GLENN  ENGINEERING  SERVICES*  INC. 

ROCKVILLE*  NARYLANO  20856  PHONE  427-38 JO 


VERSION  VA20  is  A HOOIFICATION  OF  HY-74 
««««**»«?»  WHICH  approximates  THE  EFFECT  OF  SKEW. 

HOCIFICATION  BY  D.T. VALENTINE 
OOOf  544  august*  1970. 


PROPELLER  STRESS  CALCULATION  PROGRAM  VA20 


simple  BFAH  APPROXIHATION  INCLUDING 
BENDING*  CENTRIFUGAL  AND  TORSIONAL.  FORCES. 


,«»»**,«««««««•  PROGRAM  TYPE,  x ************ 
»««««««  jHIS  PROGRAM  FEAOS  IN  THE  NEW  TYPE  DATA  CAROS.  * * * * 
THIS  AREA  RESERVES  COMPUTER  STORAGE  FOR  ALL  ARRAYS  USED  IN  THE  PROGRA 
OIHENSICN  XE12Q) 

OIMENSICN  HA(20)  ,HAl(20)*PHl(20)*PH12(20).»XUiC20>*ri(20}*Qi(20)*CP 
1HK20)  *SPHI(20)*X4(20)*AEC20}*BF(20)*PE(2.0.) 

DIMENSION  A(i3)*B<13)  *C(13)  *0(13)  *£(13)  *F(  13)  * G(13)  *HU3)  *0  (13) , 
XP(13)  *Q(13)*R(13)  *S(13)*Thj)  *U(13).*V(13)  *H(13).*X(13)  *Y(13)  »Z  (13) 


PE  (2.0) 


*UI1^I  *Kf  141  *:>  1141  *i  114>  *ut. 

DIMENSION  Rl(7.,t6),  $1(7,161 
OIMENSICN  CENTST  (7),  CENTMO  (7) 

DIMENSION  FHXC7)  ,TX(7)  *FHHX(13)  * YTX(.13)  , SKEW, (13)  ,XU(10) 

DIMENSION  VOL(7)  ,CENT4(7),A1(7),A2(7J  * X2BAR(  7.)  *CENTS2  ( T)  *82(13)  * 

X F5(T3)  *P2(13), 02(13)  ,A.A.(ld)'*BB.(  1.0).*CENT42(7)  ,CENTHS(7)  » 

X TSKEHl(7),TSKEH2(7>.jASKEHl(7).,ASKEH2(7.) 

DIMENSION  V2<13).,02Ct'3),E2(13) 

OIMENSICN  ALPHIA(7)  ' 

DIMENSION  XMT«10),XL(10)^XMH0),XT(10)*STH(10)»STLT(10) 

SET  'J<»  GChSTANTS  USED  IN  PROGRAM  COMPUTATIONS. 

Hi  1)»1-. 

H(  ?)=4. 

H(  3) =3. 

H(  4) =8. 

H(  5) =4. 

H{  6)*9. 

K(  7)=4. 

H(  8) s8. 

H(  9)s4.  ; ' ' • 

H(lb)sA. 

k(ii)=3.  . , ' . 

M(i'2)s4.'  ' ■ ‘ 

H(13)*1,  . ' ; ' 

~ ./.,RAV<  ~ ;32..1A  r.v  ...V  ^ 


, • l«t 

- - - - INPUT  DATA  OESCPIPTIOkS  - - - i - - L .• 

H.  * TOTAL  NUM3EP,  0F  f;ROP.,EL.LER^ky^^^^  ,NA.PE  IN  THI^S  BAT 

.'m  m 


O U O O CJ  O U O O U O O O CJ  O O O U O O O (J  U O O * J r.J  vJ  <„>•_)  ;j  ■_><_/  t ■ J j o o o o 


PUNCHfO  IN  COL.  1 AND  2 OF  INPUT  CARO. 

PR,  10  > PROPELLER  10  COOE.  ANT  ALPHNUHERIC  CHARACTERS. 

PUNCHED  IN  COL.  i THRU  12  OF  CARO. 

DA,  TE  = DATE  OF  RUN.  ANT  OROFR.  PUNCHED  IN  COL.  13  THRU 
ZZ«  NUH8ER  OF  PROPELLER  3LADES.  PUNCHED  IN  COL.  25  THRU  2 
YOU  MUST  PUNCH  A OECIHAL  POINT  WITH  THE  VALUE, 

VS  = SPEEC  ADVANCE  IN  FEET  PER  SECOND.  PUNCHED  IN  COL.  29 
36.  A OECIHAL  POINT  HOST  ALSO  BE  PUNCHED  HITH  VALUE 
OEN  = DENSITY  OF  PROP.  IN  LBS.  PER  CUBIC  FEET. 

PUNCHED  IN  COL.  37  THRU  Li  INCLUDING  OECIHAL  POINT. 
OIAH  :i  OIAHETER  OF  PROP.  IN  FEET.  PUNCHED  IN  COL.  L2 
INCLUDING  A OECIHAL  POINT. 

RAKE  = PAKE  OF  THE  PROP.  IN  DEGREES.  PUNCHED  IN  COL.  L6 
49  INCLUDING  A OECIHAL  POINT. 

VLL  = VELOCITY  OF  PROP.  IN  REVOLUTIONS  PER  MIN. 

PUNCHED  IN  COL.  50  THRU  56  INCLUDING  OECIHAL  POINT. 
0(1)  = LOCAL  NONVISCOUS  THRUST  COEFFICIENT  AT  XU(I). 

PUNCHED  IN  COL.  1 THRU  9 INCLUDING  OECIHAL  POINT. 
T(I)  = TANGENT  OF  THE  HYORODYNArtIC  PITCH  ANGLE  AT  XU(I). 

PUNCHED  IN  COL.  10  THRU  18  INCLUDING  OECIHAL  POINT 
Ed)  = CRAG-LIFT  PATIO  AT  XUCI).  PUNCHED  IN  COL.  19  THRU 
INCLUDING  A DECIMAL  POINT, 

P(I)  = PITCH  TO  OIAPETER  RATIO  AT  XU(I). 

PUNCHED  IN  COL.  25  THRU  36  INCLUDING  DECIMAL  POINT 
C(I)  = COkO  LENGTH  IN  INCHES.  PUNCHED  IN  COL,  37  THRU  45 
INCLUDING  A OECIHAL  POINT, 

SKEH(I)  = SKEW  VALUE  IN  INCHES  FROH  LEADING  EDGE  TO  THE 
REFERENCE  EDGE  . ALONG  THE  HELIX, 

PUNCHED  IN  COL,  46  THRU  52.  PO 
XU(I)  s VALUE  OF  -X-  AS  HE  USE  IT  IN  PROP  REFERtNCES. 

PUNCHED  IN  COL.  55  THRU  61  , INCLUDING  DECIMAL  PO 
FMX(I)  = MAX  CAMBER  AT  VALUE  OF  U(I). 

PUNCHED  IN  COL.  1 THRU  7 INCLUDING  DECIMAL  POINT 
TX(I)  » MAX  THICKNESS  AT  VALUE  OF  U(I), 

PUNCHED  IN  COL.  8 THRU  16  INCLUDING  OFCINAL  POINT 
FMMXd)  * camber  / MAX  CAMBER  VALUE  AT  UM). 

PUNCHED  IN  COL.  28  THRU  34  INCLUDING  OECIHAL  PO 
YTX(I)  = ONE  HALF  THICKNESS  / MAX  THICKNESS  AT  VALUE  OF  U 
PUNCHED  IN  COL.  37  THRU  43  INCLUDING  DECIMAL  POI 
U(I>  = VALUF  of  X SUB  L USED  IN  PROP  REFERENCES. 

PUNCHPO  IN  COL.  50  THRU  56  INCLUDING  DECIMAL  POINT 
program  NAMES  JEFINED  BELOH 
CENTST  = STRESS  DUE  TO  CENTRIFUGAL  FORCE, 

CENTMO  = MOMENT  DUE  TO  RAKE  , THAT  IS  ASKEH2  - ASKEHl. 

CENTMS  = MOMENT  DUE  TO  SKEW  THAT  IS  TSKEH2  - TSKEHl  . 

CtNT42  = COMPONENT  OF  CEN74. 

CENT4  = CENTRIFUGAL  FORCE. 

X?BAR  = CENTROID  IN  FEET, 

TSKFHl  = TRANSVERSE  SKEW  AT  POINT  OF  INTERESTS.  0.2, 0.3 

TSKEH2  = TRANSVERSE  SK£H  AT  CENTROID,  X28AR,  IN  INCHES. 

ASKEHl  = LONGITUOIAL  SKEH  AT  POINT  OF  INTERESTS.  (RAKE* 

A5KEHP  = LONGITUOIAL  SKFH  AT  CENTROID,  X28AR,  IN  INCHES 

( RAKE  ♦ SKEH  ). 

VOL  s VOLUMNE  IN  CUBIC  FEET. 

CURING  THE  TIME  OF  PROGRAM  CONVERSION  HE  USED  VARIABLE  NAMES  THAT  HE 


M, 


oooooooooooooo 


TOO  LONG  FOR  USE  ON  THE  2BH  7t9Q  FORTRAN  STSTEH*  THEREFORE  NC  KAO  TO 
SHORTEN  THEN  TO  NO  MORE  THAN  SXX  CHARACTERS  LONG.  THE  F0U0MXN6  IS 
CF  THSES  SHORTENED  NAMES 

CEMTFOR  - CUT  TO  - CENTA 
CENTST2  - CUT  TO  - CENTS2 
CENTF0R2  - CUT  TO  - CENTA2 
CENTNOS  - CUT  TO  - CENTHS 

LSXFWl  - TO  - ASKENi 

LSKEH2  - TO  - ASKEN2 


MM  aiMtC  «•••#•  XH^UT  ********  ***^**********  ****  *^**^*^  ******  ^******** 

ONE  CARO  IS  READ  TO  SFT  UP  MAIN  LOOP  AS  TO  THE  NUMBER  OF  PROP  STUDIES 
NN=l 

18  F0RHAT<flF9.6) 

ZZ=AZZ(9,2» 

VS=AZZ(7,?> 

0IAMsAZZtJ|2) 

OIAsOIAH 

V£L*60.t*AZZ(5,2)  . 

IS'C=0 

AZ?(10»20)=Q.0 
AZZ(iO|2U*Bmi) 

AZZ(iQf22)«0.0 
DO  1000  X>l»i0 

o(ii=Azz(i;2or 
t<Ti*AZ2*<lV21>‘  ‘ " 

F<I>sA2Z(l9'22) 

P(I)  = 0.0 

C{I)sAZZ(I,23)*OlA»12.0 
'5KEHa)=C(I)/2.0-AZZ(I,?A)*OIA*12.0/2.0 
lOO"  XU(I>=AZZ(I,19I 
00  50  J-2tT 

30  AZZ(J|25)~AZZ<J*1,?5» 

DO  1001  1=1,7 
FMX(I>:8,0 

1001  TX(I)=AZZ(I,?5J»DIA*12,0 

FHMX(1)*.0000 
FMMX(2)=.2712 
FMHX(3)*.A482 
F«1MX(A>  = .6993 

FMMX(5)*.8635  . ; 

FHMXC6»*.9615 

FMHX(7»=1.000 

FMHXTd»*,9786 

FMHXI9)*.8S92 

FHMX(10).*,7027 

FMHX(41»s.3586 

FMMX(12)=.i713 

.FHMXii3»  = .b000  , . 

YTX(i)s,0d60  • 

YTXT21=i:206'6.  . 

VTX(3)s.2907 

'nx(;4ji=>>ob^  / 

.YI-X,<5.),?-iA637' 

YT>C(6):*i.A952 


VTX(7} -.4962 

YTX{8)=.4653 

YTXt9)s.li035 

YTX(10)*.3110 

YTX(ll)s,1877 

YTX{12)  = .lli*3 

YTX(13)s.033S 

IN=0 

U(l)  = 0.1 
U(2)=0.Q9 
00  1002  T=3,ll 

1002  UCI)  = 0.1*FLOAT<1N» 

U(12>  =0.95 
U(131=i.C 

25  no  50  1=1,7 
DO  60  J=l,13 

= (FMHX(J)*FHX(I)-YTX(J»*TX{I)>  /C(I» 

S1(I,J)  = (FMHX(J)*F«X(IJ+YTX(J)*TX(I)>  /C(I) 

60  CONTIHUe 
50  CONTINUE 

C calculate  the  value  OF  FI  FROH  INPUT  VALUES. 

26  Fl=l. 9905* (OIAM/2. 0) **3*VS**2*PI*6.0/ZZ 
FF1=F1 

00  215  1=1,10 
5?  P(I)  = T(1)*P1*XU(I) 

215  CONTINUE- 

^ ••*«•••«•»»#«•••  MM  «MMMM  MA* MMMMf t M ««MM MMMM  M t«ff«MM 

C CALCULATIONS  FOR  CONSTANTS  USED  IN  OETERHINATION  OF  TORQUE  AJiD  THRUST 
C CALCULATIONS  OF  BENDING  MO  lENTS  FROM  THRUST  AND  TORQUE, 
no  360  15=1,2 
FlsFFl 

RA01=01AH»0. 5*12.0 
IF(I5-?)55,56,56 

55  DO  21''  1=1,10 
PEa>=PfI) 

A{II=D(I)*(1.-E«I)*T(1)  5 
AE(I).=  AII» 

B(I>  = Oa)*(E(IUTCl)> 

3Em  =BCI) 

PHI(I)=ATAN(T(I)I 

CPHIlI)*COS(PHICI)li 

SPHI(I)=SIN«PHKI)I 

HA(I)  = <SKEH(I))*<CPHia)/(KADl*XU<I))) 

XUlCI»=XUtI»*COS(HA(I) 9 
HAlCI)=C(I)/2. 

210  IFCSkewm.EQ.HAKDJ  XU1(I)  = XU<I) 

GO  TO  5« 

56  DO  57  1=1,10 
P(I)  = OF(I) 

Aa>*AE(l) 

57  BCDsBEd) 

58  Fl=Fi76«i0 
DO  69  I=:i^7 

ii*i 

IF(I5«2»62,63,63 
62  X6=XU1(I) 


GO  TO  m 
65  XO«XU(II 
6<i  X5«f 

00  66  I?*ZlflO 
X3*I3*1 

ZFaS-2165t66*66 

65  X6«Z3)«CXUia2)*X0) 

XE(I3}>8UfZ?) 

GO  TO  67 

66  X6(I3I«IXU(121*X0) 

XE<X3)sXU(I2) 

67  TltI3l»X4a3»**!I2l 

68  Ql(Z31sX4a3)*8CX2) 

T(XMSXNPUN(XC,TltI37 

Q«I>«SXttPUN(XE»QifX3) 

T(IiaTH»*FFl 

69  Q(I»*a(X}*FFl 
C 

C LCOP  MHICH  APPROXIttATES  STRESS  DUE  TO  TORSION  RESULTING 
C FROM  SKEM 

C XT(I»  * LIFT  FORCE  » XHTCI)  * MOMENT  DUE  TO  LIFT 
C 

00  111  1*1,7 
IF(I5-2)02O,83a,83O 
820  !(MT(I»a0.00 

XK*1. 9905* (0IAM72. 0) ♦*2.*VS**  2 .*PI/{2. 0*Z2) 

XA*C(I)/2.0 
X9*TX(II/2.0 
00  222  J*l,9 

XL< Jl  *A8S (SKEW  C J)-. 45*C  (J) ) 

XLCJ|sXL(J)-ABS<SKfNCI)  *.5*Ca>» 

XT ( JJ  a ( A C J)  *0 . 1*XKV  /CCOS  <PHI  ( J ) fE  ( ) 

XHCJJ aXT(J)*XLCJ» 

XMT(I)»XMT(I»^XM(J» 

222  CONTINUE 

STM<I>aXMT(I»*2.0/<PI*)JA*XB»*2.) 

STLTmaXMT(I»*2,0/fPi*XB»XA**2.0) 

GO  TO  860 
830  STM(I)*8«00 
STLTdiiO.Ob 
843  CONTINUE 
111  CONTINUE 

550  FORMATClHO,29X,lHX,i3X,4HTAUMV12y,5HtAULE,iOXp7HM  SUB  T> 
500  FO9MAT(lH,20Xy4F18.6) 

600  FORMAT (lHi,5dx,32HS»«ARING  STRFSSfeS  OUt  TO  TORSION) 

C • ■ ' ' 

c - ■ ■ • ■ ' , 

c LOOP  MHICH  CALCULATES  AREA  (A)  OF-  SECTIONS’. 

DO  230  1*1*7 
AM)*0iT  ' 

00  230  J«lfl3 
'R-’ri)-*Rrn*j»  ■ ■ . ^ 

SU»*Sl(li3^-'  ' ■ ' " - - ' ' ''  ’ 

230  A«)*AM)  tC<I»**2’^<SC  J).;  -R<  J)‘ ' )*H(U|  < / 60. 

C LOOP  MHICH  calculates  Voti  50F-  SECTIONS.'-  . ; * 
VOLTOT*0.d 
DO  241  1*1,6 


W 


■ - ' ’ >--4. J "i''-  - 


IFCI.EO.IJ  U«I*3)*XU(I) 

238  VOUI>=A(I)*(U{I»-<»>-U(I«-3»)*3IAM/288,0 

241  VOLTOT=¥OLTOT*-VOLtI» 

242  tfOLCn  =AC7)»«1.0-UCIO))*OIAH/  576.0 

243  VOLTOT=VOLTOT4-VOL(7I 

C LCOP  WHICH  calculates  CENTRIFUGAL  FORCE  AND  STRESS. 
IF«I5-2»244,246,246 

244  00  245  I=l|6 

245  A1(II=XU1(I)*€{XU1(I+1)-XU1<I))/2.C> 

A1(7|=XU1(7)*({XU1(10)-XU1  (7>)/2.0) 

GO  TO  ?68 

246  DO  247  1=1,6 

247  AlCI)  =XU(I)+(  (XU  (I«-i»-XU(I)  1/2.0) 

Al(7»  =XO(7)4((XU(iO)-XU(7))/2.0» 

C LCOP  TO  TRANSFER  CONSTANTS  FOR  OETFRHINING,  X2BAR. 

24*  00  236  1=1,7 
X25AR(I)  = 0.0 
236  A2CI)  = Aid)  * VOL(I) 

C LOOP  TO  CALCULATE  RADIAL  CENTROID  ( X2BAR  ). 

00  251  1=1,7 

X2BAR(I)  =(  (A2(l)*A2(2)+A2C3)^A2(4)fA2t5)«-A2(6)4.A2(7))  / VOLTOT  ) 
X * (OIAP/2.01 

A2(II  = 0,0 

C UNCORRECTEO  FORCE  AND  STRESS  FOR  OUTPUT  OF  ANSWERS  WITHOUT  THE  EFFECT 
C rake  AND  SK:-W  TAKEN  INTO  CONSIDERATION, 

26A  CcNT4(l)  = 0£N*4. 0*PI**2*VEL**2*VOLTOT*X2BAR(I)/ (3600, 0*GRAV» 
CENTST(I)  = CENf4(I>  / 1(1) 

251  VOLTOT  = VOLTOT  - VOL (I) 

C LOOKING  AT  THt  EFFECTS  OF  RAKE  AND  SKEW  IN  THE  PROPEtLER* 

00  263  1=1,10  ■'  ■ 

IF(I.EQ.l)  U(I»3)=XU(I) 

AA(I)  = PI»U(I+3) 

263  93(11  = SQRT(AA(l)<‘*?!^P(i)'**^) 

DO  267  1=1,7 

TSKEWKI)  = (C(I)/2.0  - SKEHd))  ♦ AA(I)/9B(I) 

K<  = 1 

146  IF(X23AR(I)-U(KK»3)*DIAM/2.0)  149,149,151 
151  KK=  KK+l 

IF(KK-lfl) 146,149,149 

149  TSKEH2(I)  = (C(KKI/2.0  - SKEH(KK))  • AA (KK) /BB(KK) 
ALPHIA(1)=ATAN(TSKEW2(I  )/(X29AR(I)  *12,0)  ) 

CENT42(I)  = CENT4(I)*C0SULPHIA(I») 

CENTMS  ( I)  .=  CENT42  (I).*  (TS.KEH2  (I).  -■  XSKEHl  ( I)  > 

CENTS 2 (i)  = CENT42(I)  / A (I) 

ASKEWl  (I).a(TSKEHl;(I).fP(l)/AAII)  );  f (UCI>3)*CIIAN*6.0* 

X fAN(RAkE*Pt/i8b.O)) 

ASKEW2(I)  = (TSKEH2(I)*P(I)/AA(I> ) ♦ (U  (KK»3)*OIAM*6,0* 

X TAN(RAKE^Pr/lSU*0>)  . 

267  CENTHO(I)  = CENT42M)  * ( ASKEH2(I)  - ASKEWi(I)  ) 

C LCOP  TO  CALCULATE  RESULT  AND  HOHENTS  FOR  BOTH  ANSWER,  PAGES. 

DO  281  1=1,7 

0(1)  s |(T(I)>CENTNO(I))*AA(I)MQ(I)-CENTHS(I);)*P(I))/BB(I) 

£ < I)  a ( ( T ( I ) *0  ENT  MO  (1)  ) *P  (I  )>  C 0(1 ) -CENT  MS  (I)  ) ♦ A A ( t ) ) / B 8 ( I ) 
02(I)=(T(I)*AA(I)-PQ(I')*P(IH  / BB(i). 

281  E2(I)s(T(I)*PCI)-Q(I)*AA*W  /'  BBdi? 

C PROGRAM  CONTINUES. 

DO  350  1=1,7 


m: 


X(I»x(j,o 

YtDxQ.Q 

00  2«»0  J*l,13 
R(J)  X fiidfj} 

S(J)  X 

X(nxX(I)»C(II**3*CS(J>  -R(J)  / 60./ACX) 

240  Y(I»=Yai*C(n**3*  CS(J)**2-R(J)**2>*HCJ)  /120./A<Z} 

GtDsO.Q 
H(I)x0.o 
no  250  J=iil3 

G(I»xGCi)*C(n**4*<SC,l)**?-ff  (J»**3>*H(J»/18  0, 

259  H(I)xH(I)+C{I)**4*CS(J>  -R(J)  S*H( J)*U(J>**2/  60. 

G{T)xG(I)-A0S<A(I))*YfII**2 
H(nxH(i)-ABS(A(in*xai**p 

F’xS (11 
no  270  J*2|13 
IF  (F2~S(JI)  260,  270,  270 

260  FpxSCJ) 

270  Xxo 

no  290  J=i,13 
IF  (F2-SCJI>  330,  280,  290 
280  K=K*1 

7(K) xU(J) 

290  KxK 


00  300  Lx1,k 

Ba)s<<C(I)*Z(Ll  -X(II)*Eail7KCI)-{CCCI)*F2-YCin*0CI))7G(I) 
X+CENTS2CI) 

e2(L)s(  (C(I)*Zai-X(I))*F2(in/HCI)-((Ca)*F2-YCI)l*02(i)>  /G(I)  + 
X CFNTSTdl 

V2(U-A8S(82(L)) 

300  V(UxABS(B(L1) 

F3xV(l» 


F4  X V2(1I 
F(I)x9(l) 
F5(II  X 82(1) 
DO  3f0  L=1,K 
Fd)xV(L) 

320  F5d)xV2(U 
GO  TO  340 
330  FdlxO.O 
F5(I»  X 0.0 


349  P(I)  = -Xd)*E(I)7H(I)r  (C(I)*S(l)-Yd>)*0d)/G(Jl)>CENTS:2  AD 
P2(I)  *-X(I)*t2d)/H;il-(C(I)»5(l)-Yd))»02dl/C(jE')^CENTSja) 
02d)  = (C(I)-X(I))*E2(I)/H(l)-<-Td))*02d)/G(;I)>CENTST^dl 

350  0(I)s(Cd)-Xd))*E(I)/H(I)-(-Y(I)>*0d)/G<J)>CENTS2  (I) 

00  100  1*1,7  ‘ . 

AREA(I)xAd) 

X8AR(II*Xd) 

YBAR(I)*Y(I)  : > 

AYFXOdlxGd) 

AYEYO„(D*H(i) 

eMxoM)xn(i) 

STRHAX ( DxAZZ.I tl>. 

ERTB(  iy*T  ( tr  ■ ■"  ' 

115  0 EMQB(i).*Q(i> 


r.-'y- 


IF (15-2)  351,352, 


352 


f 


o o o o 


351  CONTINUE 

CALL  PFKSTR  (F,P,0,STM,STLT,AZZ) 
CO  TO  T60 
35?  DUHMY=nUNMY 
NN  = NNfl 
360  NN=NNfl 
RETURN 
ENn 


SUBROUTINE  PRNSTR  ( XX, TY,ZZ,S1,S2, AZZ) 

OINENSTCN  AZZ<38,38) 

OIMENSICN  XX(IO)  ,YY(10) ,ZZ (10) ,S1 (10) ,S2( 10) 

calculation  OF  PRINCIPLE  STRESSES 
nUF  TO  TORSION  AND  BENDING. 

OIHENSICN  XI2(10) ,XI3(10) 
no  333  K=1,7 
XI2(K)=-S1(K)*SUK) 

XI3(K)=-S2(K)*S2(K) 

333  CONTINUE 
XXXsO.l 
00  L=itZ 
XXXsXXX^O.l 
00  555  H=1,3 
IF(H-2) 72,22,33 
72  XT1=XX(L) 

00=(AeS(XIi))*»2.0 

X0=00-4.*XI2(L)  ^ 

CC=<ABS (X0))**.5 
GO  TO  kk 
22  XI1=YY(L) 

GO  TO  66 
33  XT1=ZZ(L) 

66  00=(ABS(Xr'l)j»*2^.0 
Xa=b6-4.*'XI3(L!) 

CC=<ABS(X0'))**.5 
44  SIGMA  1=  (Xli’VCC) / 2V0 
SIGMA2=(Xll-CC)/2.0 
AZZ(L,M)=XXX 

azz(l,hho)5:sighai 
555  AZZ(L,M>20)=SIGHA2 
444  CONTINUE 

700  FORMAT ( IH, 20X,F12. 2, 6X,2F20.6) 

350  FORMAT  {1H0,33)(,.1HX,12X,6HSIGHA1,1(JX,6HSIGHA2) 

360  FORMAT  (lH0,2)(,99HStR£S'S,ES  AT  EACH  X ilATION.ARE  GIVEN  IN  THE  FOLLO 
XHING  OPCER^  MIDCHORO,  LEADING  EDGE,  TRAILING  EDGE.  I 
RpTURN 
■ENb 


SUBROUTINE  HEIGHT (JCtSIGHA7» IPO, IPRtHU301H,IKfHU8) 

HEIGHT  COHPUTFS  THE  HEIGHT  AND  CENTER  OF  GRAVITY,  THE  VALUES  F 
CHORD,  THICKNS,  CAHBER,  PITCH  AND  SKEHR  COHE  FROM  CHAIN.  OIAH, 
DEN,  RAKE  AND  PI  ARE  SET  IN  STRESS,  OTHER  VALUFS  ARE  COHPUTEO 
MAKING  CERTAIN  ASSUMPTIONS. 


COMMON  /HRT/  JPR 

CQHM0N/CHE1GKT/X,CH0P0,THICKNS, CAMBER, PITCH, SKEHR,OIAK,ZZ,OEN, RAKE 
1 ,SI 

1,‘>I,CT3,CPS,EP,PC,HAKE,VS,RPS,SIGHA,EAR 
OTHENSION  CH0R0(?8),THICKNS(3B) ,CAH8ER(3a),PITCHC38)  ,SKEHR(38) 
l,xn8) 

OIHENSICN  0ISTHF{38» ,A(28) 

DIMENSION  Rt9),PHT(q> 

OIMENSICN  HUaOIH(6,2) 

DATA  CNSTNTl,CNSTMT2,CNSTNt3/. 3581, .8071, .0238/ 


•♦♦♦VALUES  COMPUTED  AND  DATA  OUTPUT^*^* 

THE  HUa  OIAMEIER  IS'ASSUKiO  TO  SE  THE  OIAMETtR  TO  THE  FIRST  RAO 

RATIC  TO  BE  CONSIDERED  AND  THE  HU3  ASSUMED  TO  BE  CYLINDRICAL. 
HU^niAt'=X(l)»OIAM 

THE  HI  3 length  IS  ASSUMED  TO  EQUAL  THE  HUB  DIAMETER  AND  THE  OIS 

THE  RC-ERENCE  LiNt.  FROM  THE  HUB  FACE  IS  TAKEN  AS  HALF  THE  HUB  L 
0I<;REFL=HUB0IAMif2. 

INPUT  DATA.  ANO  ASSUMED'  DATA  WRITTEN  OUT. 

HUSLEN=HUBOiAM 
IF(HU3.EQ.O.)  GO  TO  50 
FH00IAM=HU901M(1,IK) 

AFIO1AM=HU0OIM<2,IK) 

HU3LENsHU80IM(3,IK) 

F330RE=HUB0IMC4,IK) 

AOeORE=HU0OIM<5,IK) 

OISRPFL=HU9DlM(6,I»f| 

FHORAO=FHOOIAM/2.0 
AFTRAO=AFTOIAM/2.0 
HU30I AM=X(1I^DIAH 
HU3RAO=HU6OIAM/2.0 
FR80RE=F090R£/2. 

AR90Rc=A080RE/2, 

GO  TO  270 

50  CENGRVH=HU8LEN/2c 
27 C CONTINUE 


C ♦♦♦♦WFIGH7  CALCULATION**^* 

00  10  1=1, JC 

10  A(I»=CHOkO(I)*THICKNS(I) 

C HEIGHT  OF  THE  BLADES 

0SA1=SIMPUN(X,A,JC) 

HEIGHTB=CNSTNT1*DIAH*0EN*ZZ*BSA1 
C McIGHT  OF  THE  HUB 

IF(HU8,EQ,0.)  GO  TO  200 
IFCFWODIAH-AFTOIAN»>  201,2(10,201  ‘ 

201  MEIGHHT=DEN»PI*HU8LFN/4,*(CFHDRAO«>AFTRAO)**2»(tFHORAO-AFTRA3»»*2/3 

1,0))  ' • ; 

WEI6HBR=DEN*PI*HU8LEN/4.*((FRBORE*AR3ORE)**2  + tLFR'0dREKA'RBORE)**2/3 

1.0)) 

WEI6HTH=WEIGHHT>HF1GHBR^  ‘ “ 

96  ‘ 


1 


CEN6RVHSHU3LEN-  { ( ?HUBLEN»  <FHDRA0*»2f  2.*FM0RA0»AFTH*0+3.*AFTIIA0**2 
l)/(4.*(FHORAO**2*FHOR<kO*ftFmoaFTMO**2>n*HElGHHT-CHUBLEM*(FR0OR 
2E*»2f  2.<»FRB0RE*ARB0RE^3.*A»B0RE**2>/C4.*(FR80RE**2*FRa0RE*ARB0RE+A 
?R3CRE**2)>»*HEIGHBR)/HEIGHTHI 
GO  TO  292 

20  9 «EIGHTHs:PI*HU80IAM**2*HU8LFH»DEN/4. 

202  CONTINUE 

HEIGHT  OF  THE  PROPELLER 
HEIGHT  F*HEIGHTB»HEIGHTH 

••♦♦CENTER  OF  GRAVITY  CALCULATION**** 

00  20  1=1, JC 

?0  OISTHF(I)=CNSTNT2*CAHBER(I»*SOStPITCH(l>>*CNSTNT3*CHOROfI) 
1*SIN(PITCH{I) >*DISREFL 

THE  EFFECT  OF  RARp  AND  SKEW  ARE  AODEO  TO  THE  DISTANCE  OF  THE  CE 
GRAVITY  FRO^  THE  HUB  FACE  FOR  EACH  SECTION. 

DO  30  1=1, JC 

OISTHFn>=OISTHF<l)-SKEHR{U*XtI>/2.*OIAH*tAN(PITGHm)-TAN(RAKE* 

lPI/180.1*OIAH/2.*(X(ll)-X(li) 

39  A(I|s CHORD {I)*THlCKNSa)*DlSTHFCl) 

BSA2=SIHPUN(X,A,JC) 

CEMGRV3=9SA2/9SA1 

CENGRV8=02SREFL-C£NGRV3 

center  of  gravity  considering  rake  and  skew 

CENGRV1= (HEI&HT8*CENGRV3*HEI6HTH*CEN6RVH) /HEIGHTP 
CPNGRVF=DISREFL-CENGRV1 

IF(JPR.NE.2)  go  to  53 
SHT=4. 44822 
U9=HEIGHTB*SHT 
UP=HEIGHTP*SHT 
CFL=CENGRVF/OIAf* 

C9L=CFNGRVB/OIAM 

9A=HUeLEN/0IAM 

33=FH0DIAM/0IAH 

3C=AFT0IAM/0IAH 

30=DISREFL/DIAH 

8E=HUB01AM/DIAH 

3F=FO3OKE/0IAM 

9G=ADBCR£/0IAH 

••♦♦RESULTS  OUTPUT*^»* 

IFtHUB.FQ.O.)  GO  TC  55 
PRINT  197,  U3,UP,CFL,C0L 
PRINT  108,  8A,08,BC,B058E»BF,BG 

GO  TO  53 
55  CONTINUE 

PRINT  104,  U0,UP,CFLiC9L 
PRINT  110,  BE, BA, BO 
52  CONTINUE 

NININUH  EXPANOFO  APFA  RATIO  CALCULATIONS; 

AJS=VS/-(RRS*DIAHV, 

AJA=HAKE*AJS 
AKT.=PI*CTS*AJS**.2/8  . 

AK0=CPS*AJS»*3/16. 

EARMIN=62.6*0.6*ZZ)*AKT./<SIGH|k7*UJA**24:(.7*PI>A*2))*.,t5 


IF(JPR.NE.2)  RETURN 
PRINT  lQS»EARHlNtAJS*AJA,AKT,AKQ, 


PC 


104  FORHATt  //20X, •HEIGHT  OF  BUOFS(N»**,F15.4//20X, •HEIGHT  OF  PR0P(8 

ILADEStCYLINORICAL  HUBXN)**  ,F15.4//20X, •CENTER  OF  GRA 

2VITY  OF  PROP  REFERFMCFO  FROM  HIDCHORO  OF  ROOT  SECTION  i-  FHO,  ♦ AF 
3T)/0=*  .F9.6//20X, •CENTER  OF  GRAVITY  OF  BLADES  REFERENCED  FROH 

4 NIDCHORD  OF  ROOT  SECTION  <-  FHD,  ♦ AFTI/D**  ,F9.6J 

105  FORMAT (/20X,^KELLFRS  HINIHUH  EAR*^,E10.4 

1//20X, •SPEED  COEFF  V/(ND)  JS*^, £10. 4//20X, •ADVANCE  COEFF  V( 
21-MTTI/(N0)  JAs«,E10.4//20X, •DESIGN  THRUST  COEFF  KT**,E 

310.4//20X,*TOROUE  COEFF  kQs*,Ei0.4// 

4 20X,*PROPULSIVE  EFFICIENCY  FTAO«^,£iO. 4) 

110  FORMAT(/20X,^HUB  OIHENSlONS/0*  11X,^HUB  OIAM  s*F9.6/47X,  •HUB  L 
lENGTH  =^,F9.4/47X, •HIDCHORO  OF  ROOT  SECTION  TO  AFT  END  OF  HUB  =*,F 
2<s.4) 

100  FORMAT (6F8. 41 

107  FORMAT ( //20X, ♦HEIGHT  OF  BLADES (N)*^,F15.4//20X, ♦HEIGHT  OF  PR0P(8 

1LA0ES«-TAPERE0  HUBI(H>=*  , F15.4//20X, •CENTER  OF  GRA 

2VITY  OF  PROP  REFERFMCFO  FROM  HIOCHORD  OF  ROOT  SECTION  (-  FHD,  * AF 
3T)/0=*  ,F9.6//20X, •CENTER  OF  GRAVITY  OF  BLADES  REFERENCED  FROM 

4 MIOCHORO  OF  ROOT  SECTION  {»  FHD,  ♦ AFTI/0«*  ,F9.6) 
lOe  FORMAT(/20X,*HUB  OIMENSIONS/D*  , IIX ,*LEN6TH*^,F9. 4/47X, •FHC  OIA 
lHs*,F9.4/47X,*AFT  OIAM**, F9i4/47X, •MIOCHORO  OF  ROOT  SECTION  TO  AFT 
2 ENO  OF  HU8**9F9.4/47X,*HUB  OIAM  AT  MIOCHORO  OF  ROOT  SECTI0N**,F9. 
34747X,*FH0  OIAM  OF  80RE=*,F9.4/47X,*AFT  OIAM  OF  B0REs*»F9.4) 

ENO 


SUBROUTINE  FIELD (ABCfUA,UT)UR«IK«OEXyCOH»ZRO) 

HAIN  FPV-7  FIELD  POIHT  VELOCITIES  AUG  28,1969 

DIMENSICN  Ai(792),K3(20,<»l},SlC792),Ul(792),XR(ll),X6(ii)  ,XTt(li}, 
IXSLCll) ,XST{ll),XVXCil),XUACll>,XTZ«lll,2<36>,P{36» 

DIMENSION  ABC$L1,9) 

COMMON  7HRT/  JPR 
DIMENSION  DEX(il,S,21 ,C0H(12,2) 

OIHENSICN  UA(ii) ,UT(11) yKVtiil ,UR(11) 

COMMON  A{<»2,42>,A3(24,il,3l,Sl24,ll,JJ  ,U(2li,il»5),XINPUTm,l6J 
,SINKNI20,2A)  ,XSTARCli 

IJ  yCOSI  (%2>  ,C0EXt5l  ,C0SKN(20,24>  ,6BC20,V2)  ,CLR<28).,GLRZC28)  ,C<28>  ,6 
2HAU0  0)  ,GLT  t20)  ,GTC20,42J,GTLC20),MHOB(ir>  ,NLE(28}  yNT£(28)  ,NUHf<»l} 
3,PHI  (421  ,REHARK(18J  ,R(20»  ,RVLAM(20,11)  ,RV(11)  ,RVSQf  ID  ,RZUH(20)  ,R 
3ZRV 

42(20,11*  ,RZRVI20,1  j),RZ(20),RZSQ(20)  , SB  (20, <»2), SINK  42) 
5,SL(20),SMA(42),SQLAM(20),STARI20),ST(20)  , SPACE(42).,T(24, 2)  , WEIGHT 
6(9},X  (42)  ,XCL(11)  ,XHAP(11|  ,AY,  AH, AA jAB, AC , AO, AE, AF,.AG, AH, AI 
7,AL,BUG,BLAOO,BL&OE,BBL,BB,BOG,BUB,CHOR,COSlKN,COSy,CCA,CCL,C,OECR 
8EE,OELT,0ET,O2jO516,075,O8iDELM,0,bELTA,bU,0V,0W,EX,E,EXGMU,GHUi,G 
9NU, GL MAX, GLMIN,GMU2,H,IMAX,JT,KT, LINE, HAX, HOUSE, MT, MIN 

1, NT,NX2,NSTOP,NTHICK,NX,NIN,NHAX,NHIN,NLEH,NLL,HTEM,MVV,NCOSE,MOGO 

2, PP,aQ,Q,RH,RBASE,RMAP,SLH,STH,SINIKN,SINY,SSL,TTKlCK,TP,V,XL,XP,A 
3NGLE(33),P,Z 

OIHENSICN  CN(4,3) 

EQUIVALENCE  ( A,  KB) , ( A3,AD  ,,iSl  ,S> , (Ui,U)  , (XR,XINPUT)  , 

1 (XG,XINPUT(12)),  (XiTB,XlNPUT<23)),  (XSL,XINPUT(34)), 

2 (X3T, XINFUT ( 45) ) , ( XVX, XINPUT (56* ) , (XUA, XINPUT (67) ) , 
3(XTZ,XINPUT(78)) 

NOGQsiO 
MHUB( 1)=0 
MHU3(2)*4 
MMUB(3)*12 
MHUB(4)=30 
HH'JB(5)=60 
MHUB(6)=120 
MHUB(7)=240 
MHU9(8)*360 
MHUB(9)=720 
MHU3( 10)*4320 
DO  51  N=l,NOGO 

J=NfNOGO-l 
M=NOGO-NH 

ANGLE ( J)=FLOAT ( HHUB (N) )♦. 17453293E-01 
51  ANGLE (M)=-ANGLE(J) 

NOGO=2*WOGO-1 

IX=0 

NSTOP=D€X(l,l,IK) 

NcX=0£X(3,l,IK) 

1 CONTINUE 
100  FORMAT (14) 

IF(NSTOP)  2,90,3 
95  RETURN 
109  FORMAT (18A4) 

3 CONTINUE 
NX=11 

KT=ABC(2,9) 

LINES 1 


HTslZ 

NT*12 

NTHZCK«23 

AT«0*8 

BU3<6.Q 

TTHICK*.! 

on  87 

00  87  J«l»6 

87  XIHfUTCX»J»*ABCCIf Jl 
XIHPUTIlf9)*0.0 
XIHPUT«2,9)*8.5 
XINWJT(3,9)*i.25 
Xlf<#UTOi»9l*5«0 
XlMPUT(9f9)*10*3 
XIHPUT(8»9)*15*8 
XlHPUTf7,9»*20e0 
XIMPUT (•t9)«25.0 
XIHPUT(9,9)*30.8 
XIHPUT<1,13)*35.0 
XINPUTiZf 10l«40.0 
XIHPUT(3f 101*45.0 
XIMPUT (4,101*50.0 
XIMPUT(5,10)*55.0 
XINPUKS, 101*60.0 
XINPUT(7, 101*85.0 
XIHPUU8, 101*70.0 
XIMPUT(9, 101*75.0 
XINPUTd, 111*80.5 
XIHPUTC2, 111*85.0 
XINPUT(3, 111*90.0 
XINPUT(4, 111*95.0 

XINPUT(5, 111=100.0 
00  93  1=6,9 
93  XINPUTd, 111*0.0 
00  94  1*1,9 
Q4  XINPUTd, 121*0.0 
XInPUTU, 131*0.0 
XINPUT(2,131».9120 
XINPUTd, 131*1. 004 
XINPUT (4, 131*1.0811 
XINPUT(5, 131=1. 0818 
XINPUT(6, 131*1. 0927 
XINPUT  (7, 131  *1.1 00,4 
XINPUT  (8d3l*l. 1080 
XINPUT (9, 131*1.1095 
XINPUTd  ,14)  * 1. 1120 
XINPUKZ, 141*1. 1142 
XINPUT  (3,14);*  1.1 153 
X INPU  T (4 ,141*1 .1 139 
XINPUT (5, 14) *1. 1102 
XINPUT  (6 ,14)  41 . i0>.6 
XINPUT.IT',  l41,*1.0970 
XINPU,t(8  rl,4.1*  1.0  88.8; 
XINPUTiij  i4)»i.0T4^^^^ 
XINPUT  (lfi5)*lV0.5^ 
XINPUT  (2,  iSipl.Q  3836; 


XINPUT{3,15)*1.0102C 
XINPUT(^,l5»*l,l«*f 
XINPUT(5,15)*0.0 
00  SB  1x6,9 

88  XlNPUT(I,15lx0,0 
DO  89  1x1,9 

89  XIHPgT(l,l6)xO.B 
00  15  1=1,9 
Z(I)  = XIHPUTtI,9) 

7<I>9)=XINPUT(I,10) 

Z(I+18)*XINPUTtl,lll 
ZCIt27|xXlNPUT(I,12l 
Pa»  = XIHFUTa,iJ) 

P(H-9)=XINPUT£I,m) 

P{I»18)=XINP0T{I,IS) 

15  P{I«-27»=XINPUTa,16> 

IFtNTJ  52,99,53 

52  NT=-NT 
NHELIX=-1 
GO  TO  54 

53  NHELIX=1 

54  RHxXRCl) 

0£LT=.17453293£-Qi*8UB 
IF{XT8<2n  99,32,33 

32  EUGxXTBd) 

00  N=1,HX 

34  XT3(N»=flU6/XR(N) 

33  IF(XVX(in  99,35,36  . . ' 

35  no  37  N=1,MX 

37  XVXtN>=1.0 

36  IF£XUA(2)»  99,38,39 

38  aUG=XUAll) 

DO  40  N=1,NX 

40  XUa<Nt  = (XTB(N)*XR(NJ/ByG-XVX(N))/’(1.3+XT3{H)**2) 

39  IF{XT2(2))  99,41,42 

41  BUG=XTZfl) 

00  43  N=1,NX 

XTZ(N)=aUG*(1.0-XR(N) ) - 

42  00  4 N*1,NX  , ■ 

XMaP(N)=COMAP(XRlN),RH) 

XGL(N)=XR(N)*XTB(N)  • 

9UG=XGL{N)*{XSt(N)-XSI.CNH  • 

IF(3UG.GT.0,000i)  GO  TO  55  . 

XSTAR(N)=2.0*XSTAR(N-l»-XSTAR{N-2>  . — 

GO  TO  4 - L.  . ' 

55  XSTAR(N»=XTZCN»*0El.T*(X(CCNJ**2*XR(N)^*21iMXVXCN)/»:XUA',CN))/ 
1 (XGL(N)*(XST(N)-XSL£Nn)  - ; - 

4 CONTINUE  ’ ' ■ ; * i , 

NX2=NX-2  ^ 

DO  5 Ixl,,HX2  ' ./ 

0(1,1  )=XG(I+i)  ; . f 

00  5 j=i,NX^  ' 

5 A(I,vl)=SIN(FtOAt(J)*XHAP<im>  , 

IFf-XG.t'l-}')  99,,6;-jT  ^ ^ 

6-  OE-T=i.b  ■ ‘ •■■■<  r 

call  H,ATiNV(A,NX2iB,l,0ET,H0USEV 
GO  TO  £8  ^9#, HOUSE 


Mdi,. 


i COHTINUF 

»^(NX2«1, 11*0.0 
8(NX,1)*0.0 
GO  TO  11 
7 XG(1I=0.C 

XG(NX)=f.Q 
00  9 1*1 t NX? 

XG(I»li*0.0 

no  9 J*1,MX? 

9 XG(I+l)sXG{I*l)**atJ>*§U,l) 

11  IF«JPR.NE.3)  GO  TO  60 
300  FORHATdHlJ 
WRITE  (6t3Q0t 

WRITE(6tl02)  (CQt1(N,XK),N«l,12)  »NX,aKE>HT,8US>lCT»AY 
WRITE  (6»lC;»}NT,NTKICK*<XR<N}  «XT8fN)*XGL{N»y  XSi.£N)(XST<N}»XC(N),  8 
1 (Stl) |XVX(N} f XUA«N) fXT?(Nl sNxl«NX) 

WRITE  (6»103}TTHICW,(Z(N}|PfN)  »N*i»NTHICiO 
102  FORNAT C////38X,32H» ••♦••••••  KXT-FPV-T  SSHPROPiLLE 

IR  FIELD  POINT  VE10CITIES/25X, 12A6/11X, ZlHHUMiE^  OF  INPUT  RADII, tX, 
2I2,?flX,18HTyPE  OF  CHORD  L0A3,5X,Il/ltX, 19HLATTICE  ARRANGEHEMT ,32X, 
31TH(1=NACA  A-SERieS>/13X,26HN0.0F  FUU  RADIAL  SPACES  I2,2ix,17H(2 
L=EILIPTICAL  »/i?X,?SHAN6ULAR  SPACING  -DEGREES-  F4, 1,19X,17H (3*FL 
SAT  PLATE  )/llX,16HKUMBER  OF  BLADES, 12X, 12, ZIX, 16HA-0ESIGMATI0H  0 
6F/11X,18HN0.0F  FIELD  POINTS, 33X,17HTfPE  1 CHORD  LOAO,6X,F4.2) 

104  FORMAT (llX,14H9FTMEfc«  BLADES, 13X, 13, 21X,25HH0. OF  POINTS  IN  THICXNE 
1SS/11X,34H(NEGATIVE  FOR  HELICAL  COOROIHATSS) ,17X,16HF0RK  INPUT  TAB 
2LE,6X,I2///13X,1HR,6X,22HTAH  BI  LANBOA  I SL,7X,2HST,rX,tHG, 8X, 
T’5HC(N)  VA/VS  UA»/¥S  TO// (7X,6F9. 4,Fll,6 ,Ffl.3,F6.4,F9.4M 

101  F0?MAT(///24X,S4HVEL0CITY  DISTRIBUTION  OF  2-0  THICKNESS  FORK  WITH 
1T1/L=F6.4//29X,51HFERCENT  CHORD  VELOCITY  PERCENT  CHORD  VELOC 
?ITY//(33X,F6.2,6X,F5.3,llX,F6.2,6X,F5.3)) 

6?  CONTINUE 

DO  44  N*1,NTHICK 

44  Pi[N)=<c(N)-l.O)/TTHICK 
KAX=MT+3 
“IN=HT*2 

OELMs ( 1. 0-RH) /FLOAT «MT) 

0«=nEL«*0.l25 

0?=OELH*0.5 

0516=OFLM*0,3129 

07S=OEL«*0.7S 

RZ(1)*RH*08  ' ' 

RZIMAX)sl.0-O8  ’ ' t . 

RZ(2t=‘fh*02 
PZ(NIN)  = 1,0-02 

»ZC3) =RH+D£LH  ^ * 

R(1)»RH*D516 
R1HIN>  = 1.0-05,16 

P(2»-RH*D75  ' ■ 

R(nfl)'*l.0-D75 

R3ASE»RH-02  *'  ■ - ■ . : 

DO  12  M=3,MT  '■  ' 

Art=H-l 

RZ(Hf  l)=RH+AM*OELM-  : * ' ” ' v 

12  R('H>=ReAk«-AM*dELN  " ■ i - ■ - . 

no  13  'H*1/MIN-  ' ,u  , . 

RMAPiCdHAP.CRi'NlYRH).  ~ ,v.;  - r ■ - 


SL  sFILtINCRMAPf  mP#KSL  »M« 

STtHI *FILLINCRHAP, XHAPfXST  fHXI 

GLRfM}=fILLlHIRCH)  tXRt*SI.>MX> 

STAR( H) *FILLIM(R(HI »XR,XSTAR*  HX) 
GLRZ( M>  *FILLIH  CRZi HJ  t XR*XCL#HX5 
G{M)*O.I 

00  13  JsltHXZ 

13  G(rt)*C*{Hf+BCJ,l>*SIH(FLOATtJ»*RHAP> 
6LRZ?M*X)sFILLlN(RZ{HXX>*XR»XSl.»HXI 

00  17 

00  17  Jfsl»23 

G9(rt«NI»0«0 
GT(H»N)s0*0 
17  SB(H»H)sO«0 
CALL  FAN 

00  1^1  N»l»41 

KUHCH)xN 

DO  t<» 

m k3(H» NJ*i OOO6B»O*G0tN»NJ 
00  45  K«l,41 

00  45  N*!,!? 

45  K8«H,N»*1(10060.9*6TCH*M) 

00  46  N«lf41 

00  46 

46  K9(rtf N)*1000d0.0*SB{M,N» 

2 JTsO£X«2,l»IK) 

N€XsN£X-l 

1XSIX4-1 

£X=0tX(lX,2,IK) 

00  105  J=1»J7 
105  R^(J)*0EX(J,3|IK) 

9UGS0.0 

IFINHELIXJ  47,48,46 
47  BUG=-£X/XGLfNX) 

48  0lAOD=6«2331853/FLOAT CKT> 

8U3'BLAC0/FL0AT<NT) 


50 

49 


18 


19 


0LAOE=9UG 

00  49  N=1,NT 

T(N, 1)--9UG/«174S3293S-01 

00  50  K=i,XT 

C0S<N(K,N)=C0S{BUG) 

SINKN(K,NJ=SIH(BUG) 

eUG'BUG»eLAOb 

BUGsBLAOE'BUB 

ELAOE=BUG 

00  18  Hsl,MAX 

RZSQCH)*RZ(H»^*2 

RZLAM(M)=RZiM)*GLRZ<«) 

DO  IS  J*1,JT 

RZRV{M,J)»PZ(H>*RV<JI 
RZRV2 (H,J)=2.0*R7RV(Hf J) 

00  19  Jsl,JT 

RVSOi J»*RV(J)**2 

00  19  HaliMAX 

RVLAHCMf  J»=RV(JJ*GtR2'*.H) 
BUGsEX/<«l7453293£-0i^XCL‘CNX)J 

00  31  «-l,NT 


31 


TCH»2>*TCNfli-SUr. 

D»«XQL(NX&/m,e»XGLfHX)*»2l»XUA{NX}«1.0» 

02*9.  SG96Q«)%*D8/XVX  fKX) 

8L*!)€*PLn*T(KT> 

HST0P*I*9S<HST0P» 

60  TO  f20tait2it29>tNST(P 

20  CALL  LTPAXL 
00  23  J*l,JT 

23  CALL  SUIIOUT(EX,ltV(J}yT»StSfNTt JfltCNAi<X*ly02f0S»8LADEfCNl 

21  GOTO  (1*2<»,2S,2U  yNSTOP 
2k  CALL  VORCFS 

00  26 

26  CALL  SUX0UHEX,RV(J}yT»UyStNryJyR£KAtK»2y02»0et8LA0£»Crn 
00  22  J*1,JT 

22  CALL  SUn0UT<FX,RV(J}«TyA3,S;NTyJf^EHAAX,SyQ2»0B)BLA0EyCN) 
25  GO  TO  (itii27»27),KiArOP 
27  CALL  97RAIL 
00  28  J*1?JT 

?!  CALL  SUn0Ur(FXyev{J)»TyU»StNTtJ»REHARXy4>02*08«8LAQE»CN} 
GO  TQ(1,1,1,29),NST0P 
29  PQ  30  J*i,JT 

call  SUnOUT(FX,PV(J> tT,S»S»NT,JyREHARK}5t02»C8yBLA0FfCH) 
UA(J)  sCNdyi) 

U9(J) sCN(l,2) 

33  UT(JJar»l(l,3» 

IF(NEX)  20Q|20l}200 

200  N5TnpaU 
GO  TO  1 

201  NSTOPal] 

GO  TO  1 
FN*1 


c 


FAN  ♦♦♦  FPV-7  FItLD  POINT  VcLOClTlES  AUG  20,1969 
SUBROUTINE  FAN 

OIHENSICN  AK792)  ,KB(20,A1)  ,S1  (792)  ,UK792>,  XRCii)  yXG(li}  ,XT8(li>  , 
IXSL(ll)  ,XST(il)  ,XVX(li},XUAUl)  ,XTZ(li)  ,2(36)  ,P(36I 
COMMON  A(<»2,(>2)yA3(24,ll,3),S(2<f,ll,3),U(2A,ll,3l,XZNPUT(li,l6},8( 
142,2) ,SINKN{2Q,24) ,XSTAR(11 

1) ,C0SI(42) ,COEXC5} ,COSKN{20,24) ,6BC20,42) ,GLR(20) ,GLRZ(20I ,6(20),G 
2MA(10Q) ,GLT (20) , GT (20 , 42) ,GTL(20} ,HHUB(17> ,NLCf20l yNTECZC) ,NUN(41) 
3,PHI(42),R£MARK(18),R(20),RVLAK(20yll),RV(il),RVSQfii),RZLAH(20),R 
3ZRV 

42(20,  11)  ,RZRV(20, 11), RZ(20)yRZSQ(2S), 88(20,42), SINK42) 
5,SL(20)9SHA(42},SOLAH(20),SrAR(20),ST(20}  ,SPACE(42) ,T (24,2) , HEIGHT 
6(5),X (42) ,XGL(ll),XKAP(ii) ,Ar, AN,AA,A6,AC,A0,AE,AF,AG,AK,AI 
7,AL,3UG,9LAOQ,6LAOE,BeL,BB,BOG,BUB,CHOR,C.OSIKN,COSr,CCA,CCt,C,OEGR 
eEE,QELT,0ET,02,0918,076,08,0ELR,0,DELTA,DU,Dy,0H,£X,£,EXCNU,GHUi,G 
9NU,G4 MAX, GLHIN,GHU2,H,IHAX,JT,KT, LINE, HAX,HOUSE,HT, KIN 

1, Nr,NX2,NSrOP,NTHICK,NX,NIN,NMAX,NNIH,NLEH,)M.L,NTEN,NVV,NCOSE,NOGO 

2, PP,00,Q,RH,R8ASE,RMAP,SLH,STN,SIN|KN,SINy,SSL,TYHICK,TP,y,XL,XP,A 
3NGLP(33),P,Z 

EQUIVALtNCE  (A, KB) , ( A3,A1) , ( S1,S) , (Ul, U) , (XRyXlNPUT) , 
1(XG,XINPUT(12) ), (XTB,XINPUT(23)), (XSL,XINPUT(34)), 

2 ( XST , X INPUT (45) ) , ( XVX ,XINPUTC56) ) , ( XU4 , XIMPUT (67) ) , 

3(XTZ,  XINPUT(7«)J 
00  30  N=l,20 

BUG=FLCAT(N)*OFLT 
PHT(N+21)=aUG 
M0USF=21-N 
PHl(MOUSe.)s-BUG 
SINI(N«21)>SIN(BUG) 

SINI (MCUSE)=-SINT(Nf21) 

COSI(Nt21)sCOS(BUG) 

30  C0SI(M0USE)=C0SI(N+21) 

PHI(21)=0.0 
SINK  2l)  = 0.0 
COSI(21)=1.0 
00  1 M=l,MIN 

P.OOT=  SORT  (R(H)**2f  GLR  (H)  »»2) 

0=0ELT*R00T 

SLH=SL (M) /ROOT-0. 25*0ELT 
STM=ST (M) /ROOT+0 .25»0ELT 
00  2 N=l,41 
NLFH=N 

IF  (PHI(N)-SLM)  2,2,4 

2 CONTINUE 

4 DO  5 N=NLEM,41 
NTEM=N-1 

IF  (STM-PHI(N))3,3,5 

5 CONTINUE 

3 NyV=NTE«-NLEH*l 
XL=-SUM)  +R00T+PHI  (NLEM) 

CHOR=STCM)-SL(H) 

IF(NVV-l)  99,33,34 

33  GHA(1)=1.0 
SMA{1)-0.0 
GO  TO  35 

34  CALL  CHORD 

35  NLE(M)=HLEM 
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NTE(n)=hTEH 
< = 1 

00  I N=NLtH,f4TEM 
GTCM,N)=Gth)«GM4CKI 

1 X=<*1 

=»UG=0,f‘ 
f'O  17  S=if23 
CT(1,K)=-GB(1,N>*BUG 
17  QUG=GT(1,N> 

0UG=BU6*G{1) 

GTUm=-G(l> 

NTFH=MaX0(21,NTF(l)-l>  . 

no  16  N=?1,NTEM 
GT(i,N) x*G6(l,N) «8UG 
IS  8UG=GT(1,m) 

no  21  «=2,«IM 

R'JG=C.O 

no  '’0  N=1,2C 

GT(‘^,N)=G3tM-l,N)«G3(H,M)tr.l}& 

?0  3UG=oT(H,N) 

&TL(«)  =G(H“i)-G(*«) 

3UG=3UG-GTL(H) 

fJT.1=rtA)(C{?l,NT£  (H-ll-i,NTf  <«)-l> 
00  21  H=21,N7FH 

GT(M,N)=G8(M-1,N)-G3{«,N)*1UG 
?l  3U0=GT{M,N} 

OUG=C.T 

00  2?  N=l,20 

GT (MflX,N)=G3(MIN,H) fOUG 

22  0UG=oT ( KfaXjNl 
0UG=3'JG-G(MIN) 

GTUHaxisG(MIN) 

NTrMryflx-HeijNTFCrtT*))-!) 

00  23  N=21,NT£M 

GT(MAX,M  = &q{MIN,N)*OUG 

23  0UG=GT{MAX,N> 

RHTIJPN 

09  STOP 


k 
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C CHORD  FPV-7  FIELD  POINT  VELOCITIES  AUG  28,1969 

SUBROUTINE  CHORD 

OIHENSION  AK792)  , K3t  20,*»1)  ,511792)  ,U1(792)  ,XR(11)  ,XGC11}  ,XTS(11) , 
IXSL(ll) ,XST(11),XVX(11),XUAC11),XTZC11) ,ZC36) ,P(56) 

COHiON  A{i»2,i»2),A3(2<»,ll,3),S«2%,ll,3)  ,0124,11,5)  ,XINPUT (11, 16)  ,B( 
142,2)  ,SINKN(20,24) ,XSTAR(11 

1) ,COSI (42) ,C0EX(5) ,COSKN(20,2%) ,GB(20,42) ,GLR(20) ,GLRZ(23) ,G(20) ,G 
2Hft(10  0),GLT(20),GT(20,i,2),6TLC20)  ,HHUBtl7)  ,NLE(20)  ,NTE(20)  ,NUHC41) 
3,PHI{4?),RFKARK(18),R(20),RVLAM(20,11) ,RV(11) ,RVSQ(11) ,RZLAM(20),R 
3ZRV 

42(20,11) ,RZRtf(2t,ll),RZ(2fl),RZSO(20),SB(20,42),SINI(42) 

5,SL(2  0) ,SMA(42) ,SOLAH(20) , STAR (20) ,ST(20)  ,SPACE(42) ,T (24,2) , HEIGHT 
6(5),Xi42),XGL(ll),XMAP(ll) ,AY, AK,AA,A8,AC,A0,AE, AF,AG,AH, AI 
7, AL, BUG, BLAOO, BLADE, BBL,ee, BOG, SUB, CHOR,COSIKN, COSY,  CCA, CCL,C,DE6R 
8Ee,0ELT,DET,Q2,O516,075,08,0ELH,O,DELTA,DU,0V,0H,tX,E,EX6NU,GHUl,G 
9NU,GLHAX,GLMIN,GMU2,H,IMAX,JT,KT, LINE, MAX, MOUSE, HT,HIN 

1, NT,NX2,NST0P,NTHIC(C,NX,NIN,NHAX,NHIN,NLEM,NLL,NTEH,NVV,NC0SE,N0G0 

2, pp,QQ,Q,RH,RBASE,PHAP,SLM,5TM,SINIKN,SINY,SSL,TTHICK,TP,V,XL,XP,A 
3N6LP(33),P,Z 

EQUIVALENCE  (A, KB) , ( A3,A1) , ( Sl,S) , (01, U) , (Xfi,XlNPUT ) , 

1 (XG,XINPUT(12) ), <XTB,XINPUT(23)) , ( XSL, XINPUT (34) ) , 
2(XST,XINPUT(45)) ,( XVX ,XINPUT (56) ), (XUA , XINPUT  (67) ) , 
3(XTZ,XIHPUT«78)) 

B(l,l)=1.0 

3(l,2)=Q.O 

on  1 N=1,NVV 

A(1,N)=1,0 

DO  1 Hs2,nVV 

1 A{M,N)=l.0/(FLOAT(N-H)»0.5) 

NMIN=NVV-1 

BUBslOO.O/CHOR 

00  2 M^IjNMIN 

X(M)  = XL-KFLOAT(H)-0.5)»0 
AG=X(H) ♦BUB 

2 9(HU,2)::FILLIN<AG,Z,P,NTHICK) 

GO  TO  (3,5,6) ,LINE 

3 IF(AY-0.gg  ) 4,4,17 

C CASP  1 CONSTANT  LOAD 
17  DO  7 M=2,NVV 

7 B(M,1)=0* (ALOG(1.0-X(H-1)/CHOR)-ALOG{X(M»1)/CHOR) )/CHOR 
GO  TO  11 

C CASE  2 A SERIES  MEAN  LINE 

4 £=1.0-AY 

00  8 M=2,NVV 

V=1.0-X(M-l)/CHOR 
Q=AY-X(M-1)/CH0R 
IF(ABS(Q)~0.0001)  14,14,15 

14  QQ=0*0 
GO  TO  16 

15  0Q=Q*AlOG(ABS(0)) 

16  PP=V*AL0G(V) 

8 8(H,1  )=2,  0*D»  ( (PP-QQ)/E-AL0GCX  (M-D/CHOR)  -1.  0)/(  (AYU.  0)*CHOR) 

GO  TO  11 

C CASE  3 ELLIPTICAL  LOADING 

5 no  9 M*2,NVV 

9 B(M,l)=4.0*O*(l,0-2,0»X({1-l)/CHOR)/CHOR 
GO  TO  11 
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t 


LTRAIL  FPV-7  FIELD  POINT  VELOCITIES  AUG  20,1969 

subroutine  LTRAIL 

OIHENSICN  A1(792),KB(2Q,41I,S1(792} ,U1 (792) , XR(li)  XG ( lit , XTB  ( 11) , 
iXSL(ll),XST(ll),XVX(ll> ,XUA(11) fXTZCll) ,Z(36},P(36) 

COMMON  A(42,(>2)  , A3 ( 24, 11,3)  ,S  ( 24, 11,  T)  ,U(24,  11,3)  ,XINPUr(ll«16)  ,3( 
142,2) ,SIN<N(20,24) ,XSTAR(11 

1) ,COSI(4?) ,COEX(5),COSKN(20,24) ,GB(20,42) ,GLR(20) ,GLRZ(20) ,G(20) ,G 
2MA(i00),GLT(20),GT(2a,42),GTL(20) ,MHUB(17) ,NLE(20) ,NTE<20) ,NUM(41) 
3, PHI (42) ,REHARK(18) ,R( 20 , RVLAH (20, 11) ,RV(li) ,RVSQ(il)  ,RZLAH(20) ,R 
3ZRV 

42(20, 11) ,RZRV(20,il) ,RZ(20),RZSQ(20) ,SB(20 ,42) ,S1NI(42) 

9,SL(20) ,SHA(42) ,SQLAH(20) ,STAR (23) ,ST (20) , SPACE (42) , T ( 24,2) , HEIGHT 
6(5),X (42) ,XGL(il) ,XHAP{11) ,AY, AH,AA,A8, AC,A0,A£, AF,AG,AK,AI 
7,AL,BUG,BLAOO,8LAOE,GSL,88,BOG,BU3,CHOP,COSIKN,COSY,CCA,CCL,C,OEGR 
8EE,0EL7,0FT,O2,O916,O75,O8,OELH,O,0ELTA,0U,OV,CH,EX,E,EXGNU,GMUi,G 
9NU, GLHA V,GLHIN,GHU2,H,IKAX,JT,KT, LINE, MAX, MOUSE, HT,HIN 

1 , NT,NX2,NST0P,NTHICK,NX,NIN,NMAX,NHIN,NL£H,NLL,NTEM,NVV,NC0SE,N0G0 

2, pp,QQ,Q,RH,RBASE,RMAP,SLH^STM,.SINIKN,SINY,SSL,TTHICK,TP,  V,XL,XP,A 
3N6LP(33) ,P,Z 

EOiJIVALENCE  ( A,  KB)  , ( A3,  Al)  , (S 1 , S)  , (U1 , U)  , (XR,XINPUT ) , 

1 ( XG, X INPUT ( 12 )), (XTB, XINPUT(23) ) , ( XSL , XINPUT (34) ) , 

2 (XST , XINPUT (45 ) ) , ( XVX, XINPUT ( 56) ) , ( XUA , XI NPUT ( 67) ) , 

3 (XTZ, XINPUT (781) 

PQ  5 N= 1,792 
Sl(N) =0.C 
BUG=40. Q 

00  1 J=1,JT 

00  1 H=1,MAX 

90G=AB«;(rV(J)-R2(H)) 

IF(30G-8UG)  2,1,1 

aUG=B0G 

TP=EX/GLPZ(M) 

CONTINUE 


i 


r> 


IF(TP-ANGLE(in  3,3,4 
4 no  7 M=1,N0G0 

NL0H=M 

IF{AN6L£(M)*TP)  7,7,8 

7 CONTINUE 

8 IF(NLOW-l)  10,10,11 
10  GMA(l) =TP*ANGLF(1) 

^ GO  TO  12 

n GMA(1)=0.0 

12  N=? 

DO  9 M=NLOH,NOGO 

GMA(N)=TPfANGLE(M) 

9 N=Nfl 

NCOOE=NOGO-NLOM+1 
no  6 I=1,NC00F 

OELTA=GMA(iH)-GMA(I) 
C0EX(1)=. 046910*OELTA 
COEX(2)=.230765*OFLTA 
C0EX(3)=.5*0FLTA 
C0EX(4)=.769235»0ELTA 
COEX(5)=.953090*OELTA 
HEIGHT (1)=.059232*OFLTA 
WEIGHT(5)=WfIGHT(l) 
HEI6HT(2)=.119657*0ELTA 
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HEIGHT  (<i)=HEIGHT  (2) 

HEIGHT <3)a.l42222*0ELTA 
DO  6 L*l,5 

GMU>GHA(I)«COEX(U 
COSY»COS(GNU) 

SINYsSiHtGNU) 

DO  6 K=1,KT 

DO  6 N=i,NT 

COSlKN=COSY*COSKN (K,N) -SINY»SINKNCK,M» 
SINlKN=SINY»COSKHCKfM|»COSY*SINKN(KtN) 

DO  6 J=1,JT 

DO  6 Hsi»HAX 

BOG=SQRT  ( { (EX-GLP7  CMI  •GNU)  ♦•2«>RVSQ(  J)  ♦RZSQCH)  -RZRV2C H,  J) * COSIKNI* 
t»3) 

B0G=HEICHT<L)*GTUH)/80G 

FXGNU=(£X-GLRZ(H)*GMUI*RZ«H) 

S(N,J,l)  = £.N,J,i)»(R2SQCH)-RZRV(H,,J)*C0SIKN)*B06 

SCN,Ji2)=S(N«J)2)«-(RZLAH(H)»SZNlKN»EXGNU».C0SZKM^B0G 

S(N,J,T)-S(NtJ93)f(RVLAH(HfJ)«RZLAH(i4)«COSIKN«-EXGNU»SINlKN)«BQG 

CONTINUE 

RETURN 

FNT 


c 


3 


4 


6 

7 
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VOPCES  FPV-7  FIELD  POINT  VELOCITIES  AUG  20,1969 

SUBROUTINE  VORCES 

DIMENSION  Al(r92)  ,KB<20,M»,S1  (792)  ,U1  (792) , XR  (111  , X6  ( 111  ,XTB(11), 
IXSL(ll) ,XST(11),XVX(11I ,XUA(11) ,XTZ{11) ,Z(36) ,F(36) 

COMMON  A(42,42},A3(24,lly3),S(24,ll,3) ,U(24,11,3),XINPUT(11,16) ,B( 
142,2)  ,SINKN(20,24>  ,XSTAfi(H 

1) ,COSI(%2) ,C0EX(9) ,COSKN(20,24) ,GB(20,42) ,GLR (20 ) ,GLRZ ( 20) ,G(20) ,G 
?Mft (100),GLT(20) ,GT (20,42), GTL( 20) ,HHU3(17) ,NLE( 20) ,NTE { 20) ,NUM(41) 
3,PHI(42) ,RFMARK(18),R(20) ,RVLAH(20,11) ,RV (11) ,RVSO(ll) ,RZLAM(20) ,R 
3ZOV 

42(20, 11) ,RZRV(20,11),RZ(20),RZSO(20),SB(20,42) ,SINI(42) 

5,SL(2  0) ,SMA(42),SQLAM(20) , STAR (20) , ST (201  , SPACE (42) , T ( 24, 2) , HEIGHT 
6 (5),X (42) ,XGL(11) ,XHAP(11) , AY , AH, AA , AB, AC , AQ , AE, AF. AG , AH, AZ 
7,AL,BUG,3LA00,BLA0E,86L,66 ,B0G , BUB,CH0R ,COSlKN,COSY, CCA ,CCL ,C , OcGR 
8EE,0ELT,0ET,O2,0516,O75,D8,DtLM,D,0FLTA,DU,0V,OH,EX,E,EXGHU,GMUl,G 
9NU, GL HA X, GLMIN, GHU2,H,IHAX,JT,KT, LINE, MAX, MOUSE, MT, MIN 

1, NT,NX?,NSTOP,NTHICK,NX,NIn,NHAX,NMIN,NLEM,NLL,NTEM,NVV,NCOSE,NOGO 

2, PP,00,Q,RH,RBASE,P.MAP,SLM,STM,SINIKM,SIMT,SSL,TTHICK,TF,V,XL,XP,A 
3NGLE(33) ,P,Z 

EQUIVALENCE  ( A,  K8I , ( A3, A1 ) , (SI ,S) , (U1 , U)  , (XR,XINPUT) , 

1 (XG,X INPUT (12) ) , (XTB,XINPUT(23)) , (XSL,XINPUT (34) ) , 

2 (XST,  XINPUT (45) ) , ( XVX, XINPUT ( 56) ) , (XUA, XINPUT  (67) ) , 
3(XTZ,XINPUT(78)) 

00  3 N=l,792 
UKN)  =0.0 

at(N) =0.0 
no  5 1=1,41 

no  1 <=1,KT 

00  1 N=1,NT 

C0SIKN=C0SI(I)*C0SKN<K,N)-SINI (I) »SINKN(K ,N) 
SINIXN=S1NI(I)*C0SKN(K,N)«-C0SI(I)  *SINKN(K,N) 

00  1 J=1,JT 
93=-2.0»KV(J)*COSlKN 
00  1 MC=1,HIN 
IF(G3(MC,I) .EQ.0.0)  GO  TO  1 
AA=6X-GLR(MC)*PHI( I) 

A9=AA**2+RVSQ( J) 

AC=4. 0*  (AA**2  + RV^Q(J)*SINIKN»*2) 

MO=MC«'l 
BUG=-1.0 
AG=0.0 
AM=0. 0 
AI=0. 0 

00  8 M=MC,M0 

1F(ABS  (AO-Q. 00001)  4,4,6 

AE=0.5*BO+RZ(M) 

AO=“0.25/AE**2 

AF=-0 . 5/AE*O.125*0B/AF**2 

GO  TO  7 

AE=AC*SOKT(RZSQ(M) fRZ (M)»BB+A9) 

AO=(2.n*RZ(M)+9B)/AE 

AF=.(3B*RZ(M)f2.0*AB)/AE 

AG=AG+nUG*AO*GB(HC,I) 

AH=AHfRUG*AD*SB(MC,I) 

AI=AI+BUG*AF»SB(MC,I) 

BUG=-0UG 

U(N,J,  1)=U(N,J,1)-kV(.))»SINIKN*AG 
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Try: 


U(N,J,^>xUtN,J,2)♦AA♦SINIKN•A6 
U(M,J,3)xU(N»J)3)-«A«C0SIKN»AG 
AJ(N,  J,l)=A3(N,Jtl>4>AA*AH 
A3(N,J,2)=A3(N,Jt2f)^(RV(J)  ♦AH-C0SIKN*AI) 
A3<N, J,3J*A3(N,J,3>-SINIKN*AI 
1 CONTINUE 

5 CONTINUE 

RETURN 
ENO 


^ pi 


C 8TRA1L  FPV-7  FIELD  POINT  VELOCITIES  AUG  20,1969 

SUWOUTINE  9TRAIL 

OINENSICN  A1  (792»,K8(  20,61),  SI  f 792)  ,Ul(792)  ,XRU1)  ,XG(11)  ,XTB{11), 
iXSL(ll) ,XST(11),XVX(11),XUA(11) yXrZCli) ,Z(36) ,Pf36) 

COHHON  A (42,62) , A3 (26, 11, 3) ,S ( 26, 11,3) ,U(26,11,3) ,XINPUT(11, 16) ,8( 
162,2) ,SINKN(2B,26) ,XSTAR(11 

It ,C0SI(62) ,C0EX(5) ,COSKN(20, 26) , G3( 20 , 62) , GLR (20) ,GLRZ(20) ,G(20) ,G 
2HA(100) ,GLT (20) , GT (20 , 62) , GTL ( 20) ,NHUB(17) ,NLE(20) ,NTE(20) ,NUK(61) 
3, PHI (62) ,R£MARK(iej ,R(2B) ,RVLAH (20, 11) ,RV(11) ,RVSQ<11> ,RZLAN(20) ,R 
3ZRV 

62(20, 11) ,RZRV(20,11),RZ(20),RZSQ(20),S3(20,62) ,SINI(62) 
5,SL(20),SHA(62),SQLAH(20),STAR(20),ST(20) ,SPACE(62) , T ( 26,2) , HEIGHT 
6(5),X(62) ,XGL(11) ,XHAP(11) ,AY, AH,AA,A3,AC,AD,AE,AF,AG,AH,AI 
7,AL,BUG,eLAOO,BL A0E,3BL,SB,B0G , BUB,CH0R,C0SIKN,C0SY, CCA ,CCL,C ,0EGR 
8EE,0ELT,0ET,02,0516,075,0S,0ELH,0,0ELTA,0U,0V,0H,EX,E,EXGNU,GHU1,G 
9NU,GL MAX, GLHlN,GhU2,H,IHAX,JT, XT, LINE, HAX, HOUSE, HT,HIN 

1, Nr,NX?,NSTOP,NTHICK,NX,NIK,NHAX,NHIN,NLEH,NLL,NTEH,NVV,NCOSE,N060 

2, PP,QQ,Q,RH,R8AS£,PMAP,SLH,STH,SINIKN,SXNT,SSL,TTHICX,TP,V,XL,XP,A 
3HGLE(3-<)  ,P,Z 

EQUIVALENCE  (A, KB) , ( A3,A1 ) , (S 1 ,S) , (U1,U) ,(XF,XINPUT) , 
i (XG,XINPUT(12)) , (XT8,XINPUTI23)) , (XSL,XINPUT(36I ) , 

2 (XST,  XINPUT(65)) ,(XVX,XINPUT(56) ) , ( XUA , XINPUT ( 67) ) , 
3(XTZ,XIHPUT(79)) 

00  3 N=l,792 

3 Ul(N)s0,0 

00  1 L=l,69 

GHUlsPHKL) 

OELTA=FHI (L»1)-GHU1 
COEX< l)aO,211325*OELTA 
COtX(2)=0.783675*OELTA 
HEIGHT (1)=0.25*OFLTA 
HFIGHT(2)=HEIGHT(1) 

00  4 1=1,2 

GNU=GHI)1«-C0EX(I) 

COSY=COb (GNU) 

STNY=STN(GNU) 

00  4 K=l,KT 
00  4 N=1,NT 

C0SIKN=CUSY*C0SKN(K,M)-SINY*SINKN(K,N) 

SlNIKN=SINY»COSKN(K,N)+COSY*SINKN(K,N) 

00  4 J=1,JT 
00  4 H=1,HAX 

IF(GT (H,L) .EQ.0.0)  GO  TO  4 

BUG=SQRT( ( (EX-GLKZ  (H) •GNU)»*2+RVSQ( J) ♦RZSQ ( H) -RZRV2 ( H, J) • C3SIKN) 
1»*3) 

BUG=HEIGHT(I)»GT  (M, D/BUG 
EXGNU=(eX-GLRZ  (H)*GNU)*RZ(H) 

U(N, J,l)=U(N,J,l)*(RZSQ{M)-RZRV(H,J)*CnSIKN)*BUG 

U(N,J,2)=U(N,J,2)+(PZLAH(H)»SINIKN+EX6NU*C0S1KN)*8UG 

U(N,J,3)=U(N,J,3)f(RVLAH(M,J)-RZLAM(H)*C0SIKN»£XGNU*SINIKN)‘^BUG 

4 CONTINUE 

1 CONTINUE 

RETURN 

PNO 
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. -'i*  *V-'-  ^ 


! I 


i 


i 
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SUHQUT  FPV-7  FIELD  POINT  VELOCITIES  AUG  21 » 1949 

SU3ROUTINE  SUHOUT(X,R,T*U#SfMT, J»^EHAtKrKODE9A,B»C»CN» 
OIM6NSICN  T {24t2)  • 0{2A»11»  3)»  S (24»  11>  31  »NEHARK  tlAI  *SN{4*3) 
OIHENSION  CN(<9>3) 

GO  TO  a,2,3j4tS)»K0nF 

1 GO  TO  21 

2 GO  TO  21 

3 GO  TO  20 
«♦  GO  TO  20 
« CONTINUE 

20  00  12  I«l»3 
L=J*I-1 

12  CALL  HAHMAN(UC1» J»I)»T,CN(1»I) »SN(1,I>»NT,L) 

SPsA*  CB*SN( 1, 1)»R*SNC1,3> ) 

COsA* (B*CN<2fl)fR*CNC2tS)) 

PKsSQRT  (SP*»2>CP*»2> 

IF(ABS(SP)-fl.OOOOi)  i3,la,Hi 

13  PT=90.0*(1.0-SIGN(1.0»CPn /C 
GO  TO  15 

14  IF(ABStCP) -O.aOOOl)  16jl6*17 

16  PT=90.C*t2.0-SI6N{l,0,SPJ)/C 
GO  TO  15 

17  PTs(90. 9-ATAN (CP/SF) / • 17453293E-01) /C 
IF(SP)  18|18,15 

13  “TsPT+iaO.O/C 
19  CONTINUE 

GO  TO  (9f  8f 8,8,9) ,KOOE 

8 00  7 N*1,NT 

00  7 K=l,3 

S(N,J,K)=S(N,J,K)«-U(N,J,K) 

7 CONTINUE 

9 RETURN 
ENT 


; ***  COMAP  ♦♦♦  FPV-7  FIELD  POINT  VELOCITIES  AUG  20,1969 

FU^CTION  COHAP(TEMP,RH) 

IF(TEMP-.999)  1,1,2 

2 C0HAP=3. 1415926 

GO  TO  19 

1 CN=C1»0*RH-2*0»TEMP)/<1»C“P.H) 

IF(ABS<CN)-.OOOOl)  17,17,18 
17  COMAP-1. 5707963 

GO  TO  19 

l«.  CTN=SQRT(1.0-CN**2)/CN 
CO‘<APsATAN(CTN) 

IF(CTW)  20,19,19 
20  C0MAP=CCMAP+3. 1415926 
19  RETURN 

C END  OF  ***  COHAP 

FNO 
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FILLIN  ♦**  FPV-7  FIELD  POIHT  VELOCITIES  AUG  20,1969 
FIND'S  YCXJ  FROF  TABLE  OF 
A9(N)  AND  OR(N)  CONTAINING  NO  POINTS. 

FUNCTION  FILLINfX,AB,QQ88FL,NO) 

QIHENSICN  a8{3),QOQOFL{3) 

C OIHEHSIQN  AS(3i,OR(3t 

ANTRA  (O901FL,Q00  2FL,a003FL,380<iFL,Q005FL,QO06FL,Qa07FL)xQ00SFL*(QQ 
lO«»FL-Q082FL>*(Qa04FL-Q003FU  / C (Q001FL-Q002FD*  (Q001FL-Q003FLI  )♦  QO 
206FL*(Q90iiFL-Q00lFU*(O00%FL-O003FL)/(  (00  02FL-QOOIFLI*(0002FL-0003 
3FU)»Q087FL*{QflB‘»FL-QOOlFU*CQ004FL-Q002FLJ/(  (Q803FL-Q001FLJ*  CQOO 
<*3FL-Q002FLn 
IFCX-ABfin  1,3,2 
3 Y=n800FL(l) 

GO  TO  99 

1 Y=ANTRA(AB(1) ,A3(2} ,A6(3t ,X,300DFL(1) ,QOO OFL (2> , QOOOFL C 3) } 

GO  TO  99 

2 IF(X-A9C2))1,6,5 

6 Y=OOOOFL(2) 

GO  TO  99 

5 00  7 1=3, NO 

H=I 

IF{X-A9(It>8,9,7 
9 VsOOOOFL'I) 

GO  TO  99 

7 CONTINUE 

Y=ANTRA«AB(M-?),AB(H-l),Aa(«)  ,X,QOOOFU«-2),QOaOFHH-l)  ,QCOOFU(M)) 
FILLIN=Y 
RETURN 
END 


HARMAN  FOV-7  FIELD  POINT  VELOCITIES  AUG  20,1969 

SUBROUTINE  HARMAN<  X ,T  , A.,B,NI,  JUMP) 

OIHENSICN  C(2A,3),S(24,3),T(2<,)  ,A«,>,B(A)  ,X(24} 

IF(JUMP|  2,1,2 
0=NT 

00  3 N=1,NT 
ANGL=FLCAT (N-l)*360.0/0 
DO  3 K=l,3 

E=FLOATCK)*ANGL*.17453293E-Ol 

C(N,K)=2.0*COS(E)/D 

S(N,K)=2.0»SIN(E)/0 

00  4 K=l,4 

A(K)  = 0«0 

B(K)=0.0 

00  9 N=1,NT 

Atl)-A(1)*X(N) 

00  5 K=l,3 

A(K«^l)=A(Kfrl)«X(N)*C(N,K) 

B(K}=B(K)«X(N)*StN,K) 

A(l)sA(l)/0 

RETURN 

END 
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V- 


C 

C 

C 

C 

c 

c 

c 


c 

c 

c 


matrix  inversion  with 

N0VEH3ER  1692  S GOOD 


ACCOMPANYING  SOLUTION 
DAVID  TAYLOR  MODEL 


OF  linear  equations 
basin  ah  MATl 


SUBROUTINE  HATINV(A,Nif e»Ml*0E7ERMf Ml 


general  FORM  OF 
OIHEHSIOH  A*. 


dimension  STATEMENT 
, ),B(  , l.IHOEXt 


»3> 


OIMENSXCN  A{A2,A2I »IM0EX(42!^» 
equivalence  ClROMjJROHl*  (ICOLUM, JCOLUH) , 


{AKAX,  T»  SWAP) 


INITIALIZATION 

H=M1 

N=N1 


I 


I’ 


I 


1 


C 

C 

c 


10  OETERH=l.O 

19  00  20  J=1»N 

20  INDEX  U»3)  * 0 
30  DO  550  I=i»N 

SEARCH  FOR  PIVCT  ELEMENT 


40  AMAX=0*0 
45  DO  105  J=1»N 

IF(TNDEX(J,3)-1>  60,  105, 
60  no  100  Ksl,N 

IFdNOLX  (K,3)-l)  90,  100, 
80  IF  ( AMAX  -ABS  (A<U,K) 
85  IROHsJ 
90  ICOLUMsK 

AMAX  = ABS  (A<J,K)1 


60 

715 

M 85,  100,  108 


100  CONTINUE 

105  CONTINUE  ^ , 

INDEX (ICOLUM, 3) =INDtX (ICOLUM, 5) ♦! 


260  INOEX(T,ll=IROH 
270  INDEX (l*2l=ICOLUM 


C INTERCHANGE  RONS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

^ 130  IF  (IRCH-ICOLUM)  140,  310,  140 
140  OETERM=-OETERM 
150  DO  200  Lsl,N 
160  SNAPS  A ( IRON, L) 

170  A(IROH,U=A(ICOLUH,L) 

200  A(ICOLUM,L)=SWAP 

IF(M)  310,  310,  210 
210  DO  250  L=l,  M 
220  SMAP=B<IR0W,U 
230  B(IROH,L)=B(ICOLUM,L> 

250  B(lCOLUM,L)sSHAP 

C DIVIDE  PIVOT  ROM  BY  PIVOT  ELEMENT 

^310  PIVOT  sAdCOLUM, ICOLUM) 

0£T£RMsOETERM»PIVOT 
330  AdCOLUM, ICOLUM)  si. 0 
340  00  350  L=1,N 


I 


t 


o o o o o o 


A(XCOLUKtL,l«AtICOLUK»L$  «P1V0T 
35;  ir(K)  380»  380)  380 
360  DO  370  L:=itH 

370  B{ICOLUn,LlxBaCpLUKyU/PXVDT 

REDUCE  SON-PIVOT  RONS 

380  00  550  Li=i)N 

390  IFCLl-ICOLUH)  i*0Qf  550,  «*00 

<,00  T=A(L1,IC0LUK) 

1,20  A(Ll,lCOLUH)*0«e 
430  00  450  L^lyN 

450  A(Ll,LtxA(Ll,LJ-A(ICOLUH,LJ*T 
455  IFCM)  550,  550,  460 
460  DO  500  L=i,H 

50  0 BCLl,L)=8(Ll,L>-aaC0LUM,L)*T 
550  CONTINUE 

INTERCHANGE  COLUMNS 

600  00  710  1=1, N 
610  LxN+l-I 

620  IF  (IN0EX(L,l)-INDEX(t,2)J  630,  710,  630 
630  JROH=INOEX(L,l) 

640  JCOLUH*lNOEXtL,2) 

650  00  705  X=1,N 
660  SHAP=A(K,JROH) 

670  A(K,JROM)sA(K,JCOLUH} 

700  A(K, JCOLUH)=SHAP 
705  CONTINUE 
710  CONTINUE 

00  730  K = 1,N 

IF(IN0EX(K-,3>  -1)  715,720,715 
720  CONTINUE 
730  CONTINUE 
10  =1 

740  PETURN 
715  ID  =2 

GO  TO  740 
FNO 
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SUBROUTINE  OLOCBZylPO) 

DIMENSION  ftZZ(38,35J 
OlfEHSICN  BZdlll 

OIHENSTON  CZ(108),TLti2),B(99l  t*^C9»9)  y AQC9,9»  ,AU«9,9J  ,AJ«C9,9)  ,AXC 
19,9> ,AYC9,9),0AC9,9),FT«9),CGC9,i) 

OI»€NSICN  EA(9t8),e9(9,8),EC(9f8>,ED(9f8J jPA{8»jPE<3,3)yPGC3,3) ,PH 
1C3,1) ,PI(3,1) jPF«3) »XZ(9y9)fMA(9»9»»XA(9»9) 

OINFNSICN  FLO(9» ,AL0<9> ,FOCT( 9) » AOCTCtJ ,F0CP(9I , A0CP<9I ,F0PTC9) , AO 
1PTC9) ,F0HT(36) ,AGHT<56) yFOPaCSe) ,A0Ma(S6l  ,FHT8(9) ,AMT3t9) ,FHQ8C9J , 
2ANQ8(9) ,FHXO(93 ,AMX0(9J  y FHT0(9) #AHYO{9J 
COMMON  /HRT/  JPR 
OIMENSICN  INDFX(3y3i 


COMMON 

CZ 

, TL 

, B 

t AP 

f aq 

, AU 

COMMON 

AH 

, AX 

, AY 

« QA 

, FT 

y CG 

COMMON 

30 

, by 

y KO 

t M3 

y CMP 

yWAyXA 

00  30  I=lyl2 
30  TL(I}=BZ(I*99J 
KO=l 

1 KOUNT=n 
KT=0 

CMP  = 0,0 
CA  = 0.0 
DO  5 I=ly9 
3(I«^80)-BZ(I481) 

5 9(I4-903*BZCI+90) 

CO  TO  S 

6 00  7 1*1,9 

9(I>80)»9Z(i5) 

7 9(I4-901*BZ(15) 

e 00  9 1=1,14 

9  B(I)=9Z(I) 

00  10  t=l,45 

10  3(I+14)*fcZ(Ifl9l 
B(60»  =eZ(15) 

00  11  1=1,9 

3(It60l=BZ(I+63l 

11  9(I«-70)=8Z(I+72> 

O000FL=8(10) 

IP(JPP«NE.2)  GO  TO  100 
HRITE  (6,12) 

12  FORMAT (IH1,44X,32HC0NTRA-R0TATING  PROPELLER  DESIGN//) 
HRITE  (6,13)K0 

13  FORMAT  (110X,5HPAGF  12) 

WRITE  (6,14) 

14  FORMAT (SSXylOHINPUT  OATA///) 

HRITE  (6,  15) (BZ(I) ,1=1,91) 

15  FORMAT  (6X,9F12,4) 

IF  (8Z(15))17,16,17 

16  HRITE  (6,  15) (BZ(I),I=82,99) 

100  KO=KO*l 

17  ECC=B(10) 

B(12»=.95 
B(13) =1.0 
B(14)=1.05 
PM=9(12) 

KK=0 

BI=3.0 
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1 


I 


I 


BJ-0.0 

9J1=0.IJ 

SJZsQcO 

Jl=1 

RSL=Bt7)/t3.Hil592r*a(«il*B{2)J 

NZlsBdfl) 

HZ2*BCil) 

no  3B0  ICl»i,2 

IF  t<r-2l 201,202*202 

201  KJ=.T 

GO  TO  203 

202  KJ*l 

203  CONTINUE 

no  20 t KCsl,KJ 
ECCsNZl 
BdOlaNZl 
CALL  SUil 

19  DO  20  H0al,9 

HA(MO,l)**NCNO,ll 
XACHO,U  = AX<MO,i} 

20  QA(HO,l)=CG(HO,l) 

21  K2=MB 

GO  TO  (25,22>,K2 

22  N=9 
JP=9 

23  FORHATdlJ 
KlsJP-N 
IF(KlU,2<»,i 

2<»  GO  TO  *2 

25  00  26J=1»9 

EA(J,KK>=AU<J,1» 

F3(J, KKI=AH5J,1> 

FC(J, KK)=AX(J,1) 

26  EO(J, KK)=AY(J,1) 

8(10)=NZ? 

CALL  SU81 
XK'='ii 

00  226  J=l,9 
EA(J,KKl=AU(J,ll 
E9(J,KKI=AHCJ,1J 
EC(J, KK)=AX(J,1» 

226  EO(J, KK|sAY(J,l) 

57  call  SU82(eCC,E9,EC,ED,HR,IP0> 
BZ(16)=HR 

PAT=PA(%>>C1*0-FT(9))**2*PA(8I 
00  58  1*1,9 
IFd-l)  400,‘»00,401 

401  IF(I-9)  402,400,402 

400  FLO(I)=.0 

ALO(I)=tO 
GO  TO  403 

402  CONTINUE 
FLO(I)=CZ(I*9)/9(I^60) 

FLO(I»  = Ba*80J/FLD(I) 
ALOCI>*CZ«I*63>7Ba^7  0» 
ALOa>=a(I+90)/ALD<l> 
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«»I3  COHTINUE 

FOCTa»*C1.0-FLDCn*C2<I>S*CZ£I+18) 

AOCT { I ) » r 1 . Q-ALD ( X ) *02 (l*54n *C2 CI+72J 
FOCP(I)a(l.C«>FLOai/CZCI>J  •CZCI*27> 

ADCPtl)  *C  1.  04-AL0  {!)  /C2  (I*54J  I *CZ  i I»ei) 
F0PT(1)*«CI+2S)*F0CT(II 
58  A0PTCI)sR{I^32)*A0CT(lJ 
FCT-SIHPUMtB(i5) iF0CT,9J 
ACTsSXMpUN(CZ(46) }A0CT,9) 

FCP*SI«PUKIB115» ,F0CP,9) 
ACP‘SlHPUN(CZC<«6)tA0CP)9} 
FCPT=SIPPUHC8(15)»F0PT#9t 
ACPT=SI*IPUN«CZC46»  ,A0PT»9* 

FFE=FtPT/FCP 

AFE=ACPT/ACP 

FCOCsFCT/FCP 

ACOC=ACT/ACP 

TFCT=FCT*a.  0-FT  (9)  )**2*ACT 

JtsJl+1 

K=Ji+12 

BU2I -P(K) 

3Y=TFCT 

200  Pfi{KC)=3Y 

IF{<I-2I204,67,67 
204  CONTINUE 

no  31Jsl»3 
P£«Jell*liO 

31  PG(J|1»*1.0 
DO  T2J=lj3 
PH(J,l)sPA(J|**2 
PE(J,3)aPN(J,l» 

PG(J,3}sPE(J|3» 

PECJ,2)=PA(J» 

32  PG(J,2)*PE(Jf2) 

PI(lfl)=PM 

PI(2,l)=9(13) 

PI(3,l)=8<14) 

CALL  HATINS(PE,3,3,Pl9i,l,0ETEP«,I0,IN0EX> 
IF(I0-n33,35,33 

33  pi^INT  34 
GO  TO  82 

34  FORMAT aeH  CTSl*  IS  SINGULAR* 

35  PK=8Q 
PLsPK^^Z 

PR=PI  Clfl)+PK*PI(2,i)»PL»‘»I<3,lJ 
IF(CA)37,36,37 

36  A4=PR 
CA=1.0 

3?  BU2)=PR 

60  CMPsl.O 

DO  61  1=1,9 

9(I»23»=BZ(I»27) 

61  B(I«>60)»BZ(I»63) 

62  00  63  1=1,9 

63  B<I+50)*A4*B(I^50> 

64  CMPsl.O 
8(12)  =1. 


llS 


-T  ' I ■— V-T- , 

'fp-f.  V.tS’ 


't 


I 

I 

i 

I 


300  CONTINUE 

67  TFCf»  *FCP+(l.Q-FTC9n**2*AC? 

TCTP  = FCPT  ♦ (l.3-FTC9})**2**CPT 
TFFF  s TCTP/TFCP 
TCOC  s TFCT/TFCP 

FHHT=8(a)*8<2)**3*3*H»159265*8(7)**2/Ci6.8*B(10»l 

ANHTsFNKT*8(10»/3(lll 

FNHQ=  (R («)*B(2I*»2*6C7)**3» 7(16. 0*8 ) 

ANMQ=FNRO*8(10)/S(liJ 

H=0 

DO  63  J3l»4 
00  68  1=1,9 
M=m-i 
K=J*2 

F0HT<HI  = (S(I*lii»-8(K«^13n*FDCT(I) 

A0HT(M)  = {CZ(i+45)-CZ(K+%‘*n»ADCT{I) 

F0R0(H)  = (B(I>1<*)-9(J(*13))/8(I*14»*FDCP(I) 

63  A0MQ(HI»(CZ(I»‘i5)-C2CK*44>)/C2a»45>*A0CP(I) 

DO  69  1=1,9 
FNT9{II=0.0 
AMT9CI)=0.0 
FHQ8(I)=0.0 

69  ANQB(I)=0.0 

FNT3(l)=SlMPUN(B(i9>,F0HT(l>,9) 

AiTB(  1)*SIMPUN(CZ(<.6I  ,A0HT(1)  ,9) 

F<<Oa  ( 1 ) =SIHFUN  (8  (15),  FDMQ ( 1) , 9 ) 

AMC9(1)=SIHPUN(CZ(46)  ,A0HQ(1),9) 

FHT3(3)=SIMPUN(8(17),FDMT(i2)  ,7) 

A1T3(3)=SIPPUN{CZ(*i8l  ,ADHT(12)  ,7) 

FMOR ( 3 ) =SIMPUN  t B < 1 7) , FOMO ( 121 , 7) 

AH09(3)=SIMPUN(CZC48) ,A0M0(12) ,7) 

F'iT8(5)=SlHPUN(3(l?),F0HT(23)  ,5) 

AMT9(5)=SIHPUN(CZ(50) ,A0HT(23) ,5) 

FM08(5)=SIMPUN(B(19>,F0MQ{23) ,5) 

AHC8(5)=SIMPUN(CZ(50) ,ACMQ(23» ,5) 

FMT8(7)=SIMPUNC8<21),F0MT(34) ,3) 

AMT8(7)=SIMPUN(CZ(52)  ,A  (PTC34)  ,3) 

F**nS  ( 7 ) =SIMPUN  (B  ( 21)  , FOPQ  ( 34)  , 3) 

AHC8(7)=SIHPUN(CZ(52) ,A0H0(34) ,3) 

00  70  1=1,9 

FHT3(I)=FNHT*FMT9(I) 

AMT9(I)=ANMT»ANT6(I) 

FN09(I)=FNNQ*FMQ8(I) 

70  AM09(I)=ANMQ*AHQe(I) 
no  71  1=1,9 

FMXO ( I ) =FMTB (I) *COS  (C Z ( I*36) ) ♦FMQBd) *SIN  (CZ ( I+36)  ) 
A1)fO(I)=ANTB(I)»COS(CZ(I»90))»AMQB(I)*SINCCZa^90)) 
F‘(YO(I)=FMT8(I)*SIN(CZ(I»36)>-FHQ3(I)*C0S{CZa*36)) 

71  AHYO<I)=AMTB(I)*SlN{CZ<If90))-AMQB(I)*COS<CZ(I^90)) 

IF(JPR«NE.2)  60  TO  91 
Hi?ITF  (6,12) 

WRITE  (6,13)K0 
WRITE  (6,72) 

72  F0RMAT(52X,14HF0PE  PROPELLER//) 

WRITE  (6,73  )FCT,FCP,FCPT,FFE,FCOC 

73  F0RMAT(9H  CTS= 1PE1C.4,9H  CPS»1PE10. 4,9H  CTP«1PE10.4»9H 

1 tF=lP£10.4,l3H  CTS/CPS*1PE10.4) 


HRITE 

7k  fORHATUH  ) 

«RITE  C6,74) 

MRITE  (€,75) 

75  FORHAT  (17X,^HXI,8X,7HEPSILON,7X,3HHTB,9X,3KHQB,9X,SHMXO,9X,3HHYO, 
ilQX,lKH,9X,4HG(H)) 

WRITE  (6,74) 

91  KO=KO^l 

00  92  Isl,9 

92  QA(2, 1)sQA(I,1)/OOOOFL 

IF(JPR.Nt.2)  GO  TO  94 
DO  76  1=1,9 

76  WRITE  (6,77  ) B(l4l4) ,FLC(I) ,FHTB(I) ,FNQBC I) ,F«XO(I) , FHYOt I) ,I ,QA(I 

1,1) 

77  F0RMAT(12X,1P6E12.4,6X,I1,5X,E12.4) 

WRITE  (6,74) 

WRITE  (6,74) 

WRITE  (6,74) 

WRITE  (6,78) 

78  FORMAT (53X,13HAFT  PROPELLER//) 

WRITE  (6,73  )ACT,ACP,ACPT,AFE,ACOC 
WRITE  (6,74) 

WRITE  (6,75) 

WRITE  (6,74) 

94  00  93  1=1,9 

93  CGa,l)=CG(I,l)/B(ll) 

IF(JPR.N£,2)  GO  TO  22 
00  79  1*1,9 

7C  WRITE  (6,77  )CZ(If45) ,ALO(Z) ,AHTS(I) ,AHQB(Z) ,AHXO(Z) ,ANYO(I) , 1,CG( 

11,1) 

WRITE  (6,80) 

80  FORMAT (IH  //////) 

WRITE  (6,81) 

81  FOPMAT (57X,5HT0TAL///) 

WRITE  (6,73  )TFCT,TFCP,TCTP,rFFE,TCOC 

82  CONTINUE 
RETURN 

PNO 


SUBROUTINE  SU31 

DIMENSION  CZ(i08),TL(12)«Bf39) f Af (9*91  * AQ(9*9) *AU(9*9) *AW(9*9) rAX( 
19*9) * AY (9 >9) *QA(9*9)*FT(9) *CG(9*1) 

OIHENSICN  P(A1) »AH(10,9) 

DIMENSION  AJ(9*9I*AL(9)»A0(9*9)»AR(9*9) *AS(9*9) *AT(9»9) *AV(9*9) >AZ 
1 (9*9) *C(9) tE(9)tP(9)*G(9)»H(8i) *8A(9*9) *BB(9*9) *8C(9*9) *60(9*9) *BE 
2(9*9) }BF(9*9) *BG(9*9) * 3H(9 *9) * HA( 9*9) *XA(9*9) 

DIMENSION  CF(9*1) * CH( 9* 9) * Cl ( 9* 1) 

OIHENSICN  6K(9)*BS(9)  *C0(  1 0* 9)  * Q (9)  *0(Si) 

OIHENSICN  IN0EX(9,3) 


COMMON 

C2 

, TL 

* B 

* A«» 

* AO 

» AU 

COMMON 

AH 

* AX 

* ay 

* QA 

f FT 

* CG 

COMMON 

BQ 

, 8Y 

* KO 

* 

, CMP 

*MA*XA 

CA=4.0 
00  A 1=1*9 

C(I)  = {( 1.0*B(15)  )-2.0*DCl4ll»)  )/(1.0-9(15)  ) 


AH(2,I)=C(I) 

C«^=C(  I)**2 
A3G=( IcO-CE) 

IF(ABC) 1,2,3 

1 A9C=“A8C 
GO  TO  3 

2 A9c=.  ocei 

3 Cn(?,I)=SQRT(Aer) 

AH{1,I)=1.0 

<4  CD(t,I)  = C.O 
AH(2* 1)=1.0 
AHC»,9)=-1.  0 
C0(2* 1)=0.0 
C0(2,9)=0.0 
DO  8 J=3,'10 
L = J-1 

00  9 K=l,9 

AH(J,K)  = (AH(L,'<)*AH(2,K))-(CD(L*K)*C9(2*K)) 

5 CO(J, K)=(CO<L,K) »AH(2*K))*(AH(L*K)*C0(2*K)) 

DO  6 1=2, « 

6 E(I)=3.1415927’/CD(2*I) 
RSL=B(7)/(3.1415927:*9(4)*B(2)  ) 

IF  {9(91)-0.0)  9*7*9 

7 DO  9 J=l,9 

8 B(J*5  0)=RSL/(8(J  + 14)*B(5)) *SQRT(3(9)*3(J^23) ) 

9 00  10  J=l,9 

1C  F(J)=B(12)*9(J»50) 

11  00  12  J=l,9 

12  6(J)=1.0/F(J) 

DO  13  1=1,9 
JN0=I-1 

00  13  J=l*9 
X=9*JN0»J 

13  H(K)=B(I>H*)/B(J*14>f‘G(J) 

DO  14  J=l,9 

14  Q(J)  = A-^AN(F(J)) 

00  25  N=l,9 
KN0=N-1 
LN0=N-1 
HN0=N-1 
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DO  25  I»l,9 
J=9*HN0*I 

16,15,16 

15  Ks9»KN0*I 
L-9*LN0»I 
0(K)=C0S(Q(II} 

P(L)=SINCQtin 
GO  TO  25 

16  S=1.0+HCJI**2 
T=SQRT(S» 

V=l. 0+G(I)*»2 
W=SQRT(m 
A^=T-H 
U=«=-XP(AE) 

R-=C(CT-1.0>/HCJ)*(G(l)/<M-1.0)))*U)**Ba0) 

AC=1.5 

A0=,25 

X=(l.  0/(2.0*8(10»*G(I> ))♦( (V/S»**AO) 

Y={(9.0*G(II**2»*2.0)  /<V**AC»+M3.0*H(J)**2-2.0)/CS**AC>) 

Z=1.0/(2‘*.0*B(10))*Y 

IF(H(  J}-G(D)  21,21,17 

17  AF=1.0f l.O/CR-1.0) 

IF(AF)ia,19,20 

18  AF=-AF 


GO  TO  25 

19  AFs.OOOl 

20  AA=X*  (l*0/(R-l,0»-Z*ALOG<AFn 
Ks9*<N0*I 

L=9»tN0*I 

O{K)  = 2.0*8<1O5**2*G(I)*H(J)*U,O-G(I>/HU))*AA  ! | 

P<U  = 8(10)*(l.f'-G(I)/H(Jn^(1.0  + 2.0*3(10>*6(I)*AA) 

GO  TO  25 

21  AG=1.0*1.0/(i .O/P-1.0) 

IF(A6)  2?,23,2<» 

22  AG=-AG 
GO  TO  Zh 

23  AG=,00ni 

24  A3=-X*(l,0/(1.0/R-1.0)+Z*ALOG(AGn 
K=9*KN0*I 

L=9*LNC*I 

0(0  = 8(10)*G<I)*  (l.O"MCJ>/GCI))»{l,0-2.0*B(10>»Ga»*ABI 
P(L)=2.0^8(10)»»2*G(I>*(1*0-GCI) /H(J)»*A9 
2’5  CONTINUE 
00  26  1=1,9 
IM0=I-1 
00  26  L=l,9 
X=9»IH0*L 
26  AP(L,I)=0(K1 
00  27  1=1,9 
J10=I-1 
DO  27  L=l,9 
K=9*JH0»L 


? 27  AQ(L,I)=P(K1 
i 00  28  1=1,9 
; DO  28  1=1,9 


DO  29  L*l,9 
29  AO(L, I)=AHCI,L) 

CALL  MATINS(AJ,9,9,Af»j9,9,0£TER«,I0,IMDEXJ 
IF(IO-lJ30f32,30 
39  PRINT  31 
GO  TO  110 

31  FORMAT  (26H  1(A)  IS  SINGULAR  FOP  Z(F)) 

32  CALL  MATINS(AO,9,9,AQ,9,9,0£TERM,IO,INOEX) 
IF(I0>1)33,35,33 

33  PRINT  3li 
GO  TO  110 

3«*  FORMAT  (25H  1(T)  IS  SINGULAR  OF  Z(F)) 

35  AS(1, 1)=AP(1,1)*AP(2,1) 

DO  38  J=2,8 

36  AS(J, 1)=AS(J-1,1)*AP( J+1,1) 

AS(9,1)=AS(8,1) 

00  57  L=l,9 

37  AM(L, l)=FLOAT(LI*AS(L,l) 

00  38  L=2,9 

K=L-1 

38  AU(K,  1)=AP(L,1)*FL0AT(K) 

AU(9,  1)=0.0 

AV(9, 1)=0.0 
00  39  L=2,9 
JslO-U 

39  AV(J» l)=AV(J>l»l)^AU(J^lf 1) 

00  40  LSI, 9 

40  AZ(L»l)=(AW(Lfl) +AV(L,1))* 5. 1415927 
00  41  1=2,8 

00  41  L=2,9 
K=L-1 

41  AR(K, I)=AP(L,I)*AH(L,I) 

00  43  1=2,8 

AS(1, I)=AP(l,I)*AH(l,IlfA«(l,I) 

00  42  J=2,o 
K = J-1 

42  AS(J,  I)«AS(K,I)«-ARCJ,I» 

43  AS(9,  I»=AS(6,I) 

00  44  1=2,6 

00  44  L=l,9 
J=L*1 

44  AW(L,  I)=CO(J,I)*AS(L,I) 

00  46  1=2,8 

00  45  L=2,9 
<=L-1 

45  AU(K,I)=AP<L,I)*rO(L,I) 

' 46  AU19,I)=0.0 

00  47  1=2,8 
AV(9,I)=0.0 
00  47  L=2,9 
J=10-L 

47  AV(J,I)=AV(Jfl,I)*AU(J+l,I) 

00  48  1=2,6 

00  48  L=l,9 
J=Lfl 

48  AX(L,I)=AH(J,I)*AV<L,I) 

DO  49  1=2,6 
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I 

r 

V 


t 


¥ 


00  «»9  L*l,9 

49  AZ(L,I)'(AHa>I)>AX(Lt^>)*Fm 
00  50  L*l,4 

J*2*L 

50  Af*(J,9)»-l«0*AP{J,9) 
AR(lp9)*AP(l,9>+AP(2,9l 
00  51  L*2,8 

51  ARCL,9J*AR(L-l*9)»APa»i,9) 
AR(9j9)xAR(8,9) 

00  52  L=l,9 

52  AS(L,9>=FL0AT(L)*AR(L»9I 

00  53  L=2,9 

<=L-i 

53  AU(K,9)=FL0ATiK)*AP(L,9) 

AU(9,9)=0.0 

AV(g,9)=Q.O 

00  54  L=2,9 

J=10-L 

54  A\f{J,9)=AV(J^i,9)  ♦AU(J*1,9) 

00  55  L=lt9 

55  A2(L,  9)  = ( AS (L|9)»AV(Lf9n*3. 1415927 

00  56  L=l,4 

J=2*L 

56  AZ(J,9)=-1.0*AZ(J,9) 

9A(1, 1|=AQ(1,1)+A0C2»1) 

00  5?  Js2f8 

57  aA(J,l)seA(J-l,l»4-AQ(J+l,l) 
9A(9,1)=8A(8,1) 

00  58  L=l,9 

58  99(L|  l)sFLOATa)*0A<L»ll 

00  59  l.s2,9 

K=L-1 

59  eC(K,l>=AQ(L,l)*FLOAT<K) 

3C(9)  D^O.O 
e0(9jl)=0.0 

00  60  L=2,9 

J=10-L 

60  eO<J, l>=BO(J+l,l)tBC(J+l,il 

00  61  L=l|9 

61  BH(L,  l)  = (8B(L,l)4-80a,l)}*3.1415927 
00  62  1=2,8 

DO  62  L=2,9 
K=L-1 

62  BA(K,I)=AQ(L,I)*AH(L,I» 

00  64  1=2,8 

B8(l,I)=AQ(l,I)*AHa,I)+BAa,I) 

00  63  J=2,9 

K=J-1 

63  B9U,I)=BB(K,I)«‘eAU,I) 

64  BB(9,I)=BB(8,I) 

00  65  1=2,8 

00  65  L=l,9 

J=U1 

65  BC(L,I)=CO(J,I)*BBa,I> 

00  67  1=2,8 

00  66  L=2,9 

K»t-1 
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66  SO{X,I)sAO(LyX}»COCL,II 
sr  B0{9fi)«a.o 
DO  66  x«2t6 

8E<9,n*0.0 
DO  68  L'2|9 

J*10-L 

66  8E{J,I)sB£(J>iiXI^BO(J4]fX) 

DO  69  I=2f8 

DO  69  '..=  1,9 

J=L*i 

69  3FtL,  I)a6H{J,I)»BEa,I) 

DO  70  1=2,6 

DO  70  L=l,9 

7 0 BHCL,n*(8C(L,II^8FCL,I))*ECl) 

00  71  L=l,4 

J=7*L 

71  AOU,9)*-1.0*AQCJ,9) 

0A(1,9)*AQ(1,9J+A0<2,9) 

00  72  L=2,8 

72  3A{L,9)=eA(L-l,9l*AQan,9) 

6A(9,9)=BA(8,9) 

no  73  L=l,9 

73  B9(L,9)=FL0AT(L» ♦8A(L,9» 

00  74  l.=2,9 

K=L-1 

74  BC(K,9)=FL0AT(K)  *A0<L,9> 

BC(9,9)=0.0 

90<9,9)=0.0 

00  75  L=2,9 

J=10-L 

75  90(J,9)sB0(J+l,9UBC(J*l,9) 

00  76  L=l,9 

76  BH<L,9l  = (e8(L,9l+B0(L,9))*3.  1415927 

no  77  L=l,4 

J=2»U 

77  BH(J,9)=-1.0*9H(J,9» 

T = 3.0 

no  78  J=l,9 

7fl  CP(J,l)  = (8(J«-23)  - <eCJ«-14)/RSL*F<jn)*(1.0-B(15/> 
on  79  J=l,9 

79  HCJ>=CF<J,1)»RSL/CC1.0-BC15>)*BU«>14))-F(J) 

00  80  M«l,9 
00  30  1*1,9 

8',  AT(I,rt)  = <H(I»*9H(H,I)-AZCM,I>»*  (FLOAT(H»7B(10>) 

81  00  82  J=l,9 

82  0(J)  = (G(J)fH(J))/B(Jfl4>'^((1.0-B(15»>/2.0) 

00  83  J=l,9 

83  CG(J,i»*CF<J,l> 

00  84  «*1,9 
L=Mfl 

00  84  1=1,9 

84  AFa,H)=C0a,I)»0(I) 

00  35  1*1,9 

no  85  J*l,9 

85  AQ(J,I)*AT(J,X)-AP(J,X) 

00  36  J=l,9 

00  86  K=l,9 


PG 


86  CHCKy J)sAQ(K,J) 

CALL  HATtNS(AQ«9t9«CGyi,l,0£rERH,X0tIN0EX) 
IF(I0-li87,89f87 

87  PRINT  68 
GO  TO  110 

86  FORMAT  (?7K  G(NX)  IS  SINGULAR  FOR  Z(F)} 

89  DO  91  J»l,9 
SUMUAxG.O 
SUNAU^cO.D 
00  90  L*l,9 

SUMUA=SUKUA»FLOA  T C L >»CG (L, 1) * A2  CL » J> 

90  SUNAU=SUHAU»FLOAT(L)«CGtLyil*BH(L»J) 

91  AU(J, l)aSUHUA/SUHAU 

ss.nos 
00  94  J=l,9 
CI(J,1)=AU(J|1)-G{J1 
IF^CI  CJ,inS2,93,97 
9?  CKJ,  l)=-i.O*CICJ,l> 

9T  8H=S-CICJ,1) 

IF(8H)95,94,94 

94  CONTINUE 
GO  TO  101 

95  T=T-1,0 
1FCT>98,96,96 

96  00  97  J=l,9 

97  6(J)=AUCJ,1» 

GO  TO  «1 

98  CONTINUE 
200  FORMAT (9F12. 4) 

99  FORMAT  (20H  TOO  MANY  ITERATIONS) 

100  FORMAT  C12A6) 

101  00  103  J=l,9 
SUMUTsO.O 

00  102  K=l,9 

102  SUMUT=SUMUT^FLOAT(K)*CG(K,1)»0HCX,J> 

103  AM(J,  1)=SUMUT/{8{10)^<1.0-BU5))) 

00  105  J=lf9 

SUMUAs0«0 
00  104  Xslf9 

1 0 4 SUMUA  =SUMUA*FLOAT ( K) *CGCK, 1) » AZ CK, J> 

10  5 AX(J,  1)=SUHUA/(B{10)*C1.0-.8C15))) 

00  107  J=l,9 
SUMG=o.fl 

00  106  K=l,g 

L=K>1 

106  SUMGsSUMG4CG(K,l)*C0CLfJ) 

107  AYU,l)aSUMG 
AYClf 1)=0.0 
AY(9,1)*0.0 
00  108  1=1,9 

80=AY  (I,  1) **2/  C 2.O*0Cin4) » 

8V=8(in4)/RSL 

loa  9K(II  = CCRV-AH<I,in*AY(I,l)*3O)*8*0 
8L*SIMP0N(BC15)»BKC1) ,9) 

Bt=9L 

6M=B(2)**2 

0P=9(7)**2 


W 


109  MS=1 

GO  TO  111 
113  «3=2 
111  RETURN 
HMO 


SUBROUTINE  SUB2I CC* FByFCfFOyHRtIPOl 

DIMENSION  CZ(10S)»TL(i2)te(99)»AP(9,9}fAQ(9»9})AU(9r9})AH(9,9)»AX( 
19,9>  ,AY(9,9I  ,0A(9,9)tFH9)  yC6l9»l) 

OIHENSICN  TABX(9)tTASV(  ,F8(9»  8)  ,FC(9»5)  »F0(9»9)  ,FUt9>  ,FVt9)  *FH(9 
1),HA(9,9}»XA(9,9) 

DIMENSION  FEC9J,FF(9)  ,FK9J,FP(9),FQ(9)tFRC9»  ,FS«9) 

OIHENSICN  FX(9i*FY(9l >GA(9)»GB(9) »GC(9t  >G0(9> >GF(9) » GG C9) »6H(9) ;HV 
1 (9)»GI(9))GX(9)tGL  (9)  ,GH(9)fSN  (9)  f HE(9)  »HKC9)  »Hi.t9)  >HP  (9)  »HQ(9>  «KS 
2C9> 

COMMON  /HRT/  JPR 
OIMFNSION  HEK9J  ,GC1(9) 


COMMON 

CZ 

1 tl 

. B 

. AP 

> aq 

, AU 

COMMON 

ah 

1 AX 

♦ ay 

» QA 

» FT 

f CG 

COMMON 

BO 

f BY 

, MO 

t H9 

» CMP 

»HA>XA 

1 RSL=8C7)/»3.14i5927»8<4)*a<2)l 

2 GU=3.0 

00  4 J=l,9 

4 FE(JI=FO(J,4)/C2.0*8(J+  14)*FB(J,4H 
00  6 J=l,9 

6 FF(J»=FO(J,8)/(2.0’B< J*14>  «F8(Jf 8H 
FI=3.0 

FJ=0.0 
F<=1. 0 

7 00  12  Jsl,9 

« FG=FC (J,4)*FE(J> 

9 FH=FC(J,8)*FF(J> 

1C  FL=<FH*(1.0-BIJ+41)))*F«' 

12  FM(J)=B(J  + 23)tFG+-FL 
00  16  J*lf9 

13  FN=FK*FCiJ,8)*FF(J) 

14  FO=FC (J,4)*FE<JI*< 1.0fB<J+41)l 
16  FPCJ)=e<Jt32)*FN*F0 

00  18  J=2,9 

18  FQ(J) =1.0-(FH(Ji/FP<J)) 

FOCll *0.0 

19  G0=(8(16)-8(15)»/2.0 
FR(1> =0,0 

00  ?1  J=2,9 

GO*{8(J»14)-8<Jfl3)J/2.0 
21  FR(J)  = (F0(J»1H-F0(  J»)*BCJ»14»*G0 
FSdJ  =0.0 
00  23  J=2f9 
23  FSCJ»=FS(J-1)»FR<J) 

FTCl) =0.0 
00  25  J=2,9 

25  FT(J)=FSCJ)/(B(J^14)**2I 
6V=0. 0 
00  27  J=2,9 

27  GV=GV»FT(J) 

FI-GV/o.O 

28  FK=(1.0fFI)*(1.0+FJ) 

00  30  J=l,9 

30  FU(J>»FX*FCU,8)*FF(JI*(1.0-B<J*41.)) 
00  32  J=l,9 

32  FY(J)  *B{J*23)«-FCCJ,4) +FU<JI 
00  34  J»l,9 
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GQ=F<*FC(J,8J 

GR^FCCJ,4I*FECJ)*(1.0»B(J»«»1> » 

J%  FH!J}*9«J*32S  ♦GQ^GR 
00  37  J»l,9 
35  GS=Ftt(JJ/F«»(JJ 
37  FX(Ji-CFV(J»/CFJS(J»*GSI>-l.O 
GT=0.0 
00  39  J»l,9 

39  6T=GTfFK(J) 

FJ=GT/«».0 

6U=GU-1.0 

40  IF  (GUI  42f4i»41 

41  FKs(1.0*FI)*Cl.l+FJ» 

60  TO  7 

42  00  46  J=l,9 

43  FY(J)=9tJ*l4)/RSL 

44  F7=FV(JI-FC(J,4> 

46  Gft(J»=F2/FY(J> 

00  48  Jsli9 

48  G3(J) =FV(J)/(FY(J)-F8{J,4> ) 

00  50  Jsl,9 

49  GC(J1=ATAN(GB<J)) 

50  GC1(J|5ATAN(GB(J)»»57. 2957795 
DO  53  J*l,9 

51  H0=2. 0*3.1415927/CC 

53  G0(J>=H0*CF0{J»4)/FVCJH*SIHtCC(J)) 

00  55  J*l,9 

55  GF(J|sF0(J,4l/CC 
00  60  Jsl,0 

56  GH=(FV(Jl/SIN<GC(J)n**2 

57  GY=2.  0*CC/3.1415<^27 

58  GG(J)sGY*GO(J)*GH»COSCGC{J)» 

6 0 GH(JI  = GY/RSL*GOCJ»*GH*SIM(GC( J))»B{J«>14> 
00  64  J=l,9 

61  HiJ=9(3)-(a(J*14>»8l2l/2,0) 

62  HH=(SIN(GCCJn/C0(7)*FV(J)  >)»'»2 
64  HY(J)-=64.31*HU*HW 

00  66  J=l|9 

66  GIU)  = ni.0-FT(J))/(1.0-FT(9)n*3(Jfl4) 
00  68  J=l,9 

68  GK(J)=GI(J)*a.0-FT(9))/RSL 

00  72  J=l,9 

69  GL(J) =2.0*F3(J,4)*FE(J)*C1,0^FT<J)) 

70  GZ=3(Ji-32)fFC( j,4)*FECJ)«(l.lJ*B(J*41») 

72  G*1(J)=GZ/(GK(J)»6L(jn 

00  74  J=l|9 

74  GN( J) *FM ( J)/ (GKC J> -FK’FBC J,6) ♦GL C J) » 

00  76  Jal,9 

75  HE(JI=ATANtGN(JH  , 

76  ^fEi<J»  = ATAN(6N(j»'l*$7. 2957795  ’* 

DO  81  J*l,9 

77  HF*2.0*3.1415927/B<11) 

78  HI*FD<Jf6l/FHCJ> 

79  HJ-(l,0fFX(J)»/<1.0-FT(9)) 

81  HK(J»*HF»HI*HJ*SlN<HEf j»l 

00  84  Jxl,9  ■ 

82  HMx9(11)*B(J<-32>*C1.0-FT(9I) 


SV  HLCJJ =(FDCJ,8)*<1.0»FXCJI> )/HH*8(J^32) 

00  85  J*l,9 

HN=2.0*8(ll>/3.i<*15927 
H0=CFHCJ»/SINtH!:«J)))'<>*2 
HF(JI*HK*HK(J) ♦H0»C0S  CHE? Jt> 

85  HQ?JI=HH/RSL*HK{J)*HO*SINCHECjn*GI(J)*(l,0-FT(9>) 

66  HR=C1.0-FT(9))*3(2) 

00  90  J*lt9 

87  HTa9<3»-GI?J)*HR/2.0 

88  HHa(SIN(HECJH/?8(7)*FHCJni-**2 
CTSIF=SIHPUNC8CJ5I  ,G6,9l 
CPSIF=S1KPUNCB{15J ,GH,9) 

CTSIAaSIMPUN(6I,HP,9> 

CPSIA=SIHPUNCGI»HQ,9l 
9C80»  =CTSIF4.(1  .0-FTC9H**2*aTSIA 

90  HStJ)=6%.31*HT*HH 

IFCJPR.NE.2)  GO  TO  120 
IFCCMP) 107,107,91 

91  WRITE  (6,92) 

92  FORMAT (IHl) 

WRITE  (6,  93)TL 

93  FORMAT  (12A6) 

WRITE  (6,  94IK0 

94  FORMAT  <110X,5HPAGE  12) 

WRITE  (6,  95»HR 

95  FORMAT (60H 

1 1PE12.4) 

WRITE  (6,  96)  FI 

96  FORMAT (60H 
IBAR  1PE12.4) 

WRITE  (6,  97)FT(9) 

97  FORMAT (60H 
lOF  11PE12.4) 

WRITE  (6,  93) FJ 

98  FORMAT(60H 
13AR  1PE12.4) 

WRITE  (8,  99) 

WRITE  (6,  99) 

99  FORMAT  (5H0  ) 

WRITE  (6,  100) 

100  FORMAT  (53X,13HFWD  FR0PEUEK7/) 

WRITE  (6,  101) 8C70J ,CTSIF,CPSIF 

101  F0RMAT(8H  KalPEl0.4,8H  CTSIkIPEIO . 4, 8H  CPSI=1PE10 . 4) 

WRITE  (6,  99) 

WRITE  (6,102) 

102  FORMAT  (5X,4HX  (F)  ,5X,10HTAN  8ETA-I,3X,8HTAN  BETA,5X»5HCUl/.0,8X,4HG( 
IS)  ,6X,6HUT/2tfS,6X,6HUA/2Vs,tX,5H0CTSt,7X,5H(jCPSr,6)(,THSI6MA  X//> 

DO  103J=1,9 

103  WRITE  (6,104)B(J«'14),GB(J)  ,G&  C J)-,GO(J)  »GF  (J)  ^FB  C3>4) , FC  ( J,4)  ,6G(J) 
1,GH(J) ,HV(J) 

104  FORMAT! 10C1PE12. 4)) 

WRITE  (6,  99) 

WRITE  (6,  99)  , 

WRITE  C6,  105) 

WRITE  (6,  101)CZ(i80),CTSlA,&PSIA 
WRITE  (6,  99) 

105  FORMAT (53X,13HAFT  PROPELLER//) 


O-A 
DELTA 
DELTA  X 
ZETA 


WRITE 

00  106  J>1,9 

106  MRITt  (6,10ii)GI(J)«CNf  J)»r>K(J)  >HK(J)  «FBCJ»S)  yFCCJ»S)  ,HP(J1  ,K 

iO(J) 

120  X0=K0fl 
107  DO  1081=1,9 
C2(IM69(I) 

CZ(I*9J»G0(II 

CZa^l8)=GG(n 

CZ(I»27)=GHCn 

CZ<I^36)=GC1(I) 

CZ(I»45)=GI{I) 

CZ(I^54J=GN(1) 

CZ(I*631=HKC1> 

CZ(H-r2)=HP(I) 

CZ(r*81»=HO(II 
103  CZCI»90)=HEiCI) 
no  1091=1,9 
T43X(I»=9<I*14) 

109  TA3Ya)=8(I*70l 

O'!  1101  = 1,9 

X=GI(I) 

CALL  OISCOT  (X,X,TAPX,TAB‘f ,TABY,30,9,0,Y) 

110  9(I»70)=Y 
GO  TO  111 

111  RETURN 

END 


if 


ooo  ooo  ooo  ooo  ooo 


SUBROUTINE  HATINS CA|NRtNi,S,NCyHl»OSTERH,IO»ZNOEX} 

EQUIVALENCF  (IROWyJROK),  (ICOLUHy JCOLUH) » (ANftXy  T»  SWAP) 

DIMENSION  AiNR»NR)y  B(NR>NC)»  INDEX (NRySl 


INITIALIZATION 

N=Nl 

M=Mi 

OcTERN  = 1.0 
OQ  20  J=iyN 
20  IN0EX(J*:«)  = 0 
00  5S0  I=i,N 


SEARCH  FOR  PIVOT  ELEMENT 


AHAX  = Q.G 
00  105  J=1,N 


IF(IHOFX(J,3)-l) 

60, 

105, 

60 

60  00  100  K=1,N 

IF(INOEX(K,3)-l) 

80, 

ICO, 

715 

80  IF  ( AHAX  -ABS  (A(J,K))I  85,  100,  100 
85  IROM=J 
ICOLUM  =K 

AHAX  = ABS  (A(J,Kn 
IQO  CONTINUE 
105  CONTINUE 

IN0EX(IC0LUH,3)  = INDEX aC0LUH,3>  ♦! 

IN0EXa,l}  = IR0H 
INDEX (T,?)=IC0LUN 

INTERCHAnGF  rows  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 


IF  (IROV-ICOLUH)  140,  310,  140 
140  DETER M=-0ETERH 
00  200  L=1,N 
SWAPsA(IROH,L) 
A{IROH,L)=AaCOLUM,LI 
200  A(ICOLUH,L)=SHAP 

IF(H)  310,  310,  210 
210  DO  250  L=l,  M 
SHAP=8aROH,LJ 
9(IR0H,L»=B(IC0LUM,L) 

250  9aCOLUM,L>=SHAP 

DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 


310  PIVOT  =A (ICOLUM, ICOLUM) 
QETERM=0ET£RH*PIVOT 
330  A(ICOLUM,ICOLUM)=1.0 
00  350  L=1,N 

350  A (ICOLUM,  U=A(  ICOLUM,  U /PIVOT 
IF(M)  380,  380,  360 
360  00  370  L-1,M 

37  0 BCICOLUM,L)=BaCOLUM,L)/PIVOT 
REDUCE  MON-PIVOT  ROWS 
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380  00  550  L1-1)N 

IF(L1-IC0LUH}  4B0»  550,  400 
400  T^4{L1,IC0LUH> 

AfLl,ICQUUH)«8.0 
DO  450  L=1,N 

45  0 A(Ll,L»*AtLl,L»-A(ICOLU!«,U»T 
IF(H|  550,  550,  460 
460  no  500  L=1,H 

500  8(U,L)=Q(L1,L»-B(IC0LUH,L»*T 
550  CONTINUE 
C 

c interchange  COLUHMS 

C 

00  710  1=1, N 
L=N«*1-I 

IF  (IN0EXCL,1)-IN0EX(L,?))  630,  710,  630 
630  JROH=INOEX(L,l» 

JC0LUH=IN0EX(L,2) 

00  705  X=1,N 
SNAP=A(K,JR0H) 

ACK,JRON)=:A<K,JCOLUK) 

A(K,JC0LUH)=SXAP 
705  CONTINUE 
71C  CONTINUE 

00  730  K s 1,N 

IF(IN0EX(K,3)  -1)  715,720,715 
720  CONTINUE 
750  CONTINUE 

n = 1 

010  RETURN 
715  in  = 2 

GO  TO  810 
END 


.fe; 


V' 

f ' 


000000>“»00000  0 000000000<-)00  0 0000000000r>0000000000  0 00  0 00  0 00 


SUBROUTINE  GMHaSINP,Nh,Y,R» 


GHH&S,  k FORTRAN  IV  VERSION  OF  SHARE  SUBROUTINE  AH  GHHAS, 
PROGRA»*HED  FOR  IBH  360  CONVERSION. 

PROGRAFHER  - M.  GOLDEN,  CODE  892,  8-20-68. 

ALGORITMH -7 

FOR  A SET  OF  Y(I)  POINTS  f O.LE.I.LE.K-IJ , WHICH  CORRESPOND  TO  A 
SFT  OF  equally  SPACED  X(I»  POINTS,  CONFUTE  THE  ACJI,  3(JI  , C(J), 
ANT  PHI(J)  TERMS  OF  THF  FOLLOHING  SERIES,  WHERE  H IS  THF  NUMBER 
OF  HARMCNICS  OESIRFO. 

H 

Y=A(0 ) fSUH{A(J)*COS (J*X) frB (J> *SIN { J*X> ) 

J=1 


OR 


H 

Y=A(0) ♦SUM(C( J)*SIN(J*X+PHIC J) ) ) 
J=1 


HHPoe  C(J)  IS  THF  AMOLITUOE  AND  PHI(J)  IS  THE  PHASE  ANGLE  OF  THE 
JTH,  HARMONIC,  THE  FUNCTION  IS  ASSUMED  PERIODIC  WITH  Y(0»sY(K). 

ARGUMENT  DEFINITION — — 

1.  NP  - -^HE  NUM8*^P  OF  INPUT  POINTS  ( K ABOVE  ). 

2.  NH  - THE  NUMBER  OF  HARMONICS  DESIRED  ( H ABOVE  ). 

3.  Y - THF  SET  OF  INPUT  POINTS  DESCRIBED  ABOVE. 

i,.  R - THF  OUTPUT  ARRAY,  WHICH  CONTAINS  A ( 0)  , 0 (0 ) ,C  ( 0)  ,PHI  ( 0)  , 
C(0)  /C(MAX),  A(l) ,A(NH),3<NH),CtNH) ,PHI(NH),C<NH)/C(MAX) . 

p.ESTPICTIONS  

1.  0 MC)5T  3E  DIMENSIONED  AT  LEAST  5*(NH»l)«-2  IN  THE  CALLING 
PROGRAM. 

2.  Y MUST  BE  DIMFNSIONFO  AT  LEAST  NP  IN  THE  CALLING  PROGRAM. 

NOTES  — — — - 

1.  R(5*NHf6»  CONTAINS  THE  CHECKSUM 
K/2 

Y(0)=A(0)  ♦SUMCACJ))  FOR  H EVEN 
J=f 

OR 

(K-D/2 

Y(0)=A(0»»SUM(A(J)»  FOR  H POO. 

J=l 


2.  R(5»NHf7>  CONTAINS  THE  CHECKSUM 
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C K/2-1 

C ¥(l)-'f{K-lJ  = 2*SUHfO(J>*SINCJ**LPHACO)) ) FOR  H EVEN 

C J*1 

c 

C OR 

C 

C (K-ll/2 

C Ym-Y(K-1)  = 2*SUM(8(J)*SINIJ*ALPHA(0)))  FOR  K 000 
C J=1 

C 

C 3.  ALPHA  IS  DEFINED  AS  2*PI/< 

C 4.  ALL  ARGUMFhTS  ARE  SIN6LF-PRECISI0N,  AND  ALL  COMPUTATIONS  ARE 
C DOUBLE-PRECISION. 

C 5.  the  checksums  rill  not  be  computed  if  a FULL-POINT  ANALYSIS 
C IS  NOT  MAOF. 

C the  phi*s  are  in  degrees*  not  radians. 

C 7,  THE  SIGN  OF  PHlCO»  (AND  PHI (K/2>  IF  COMPUTED)  IS  MADE  TO  AGREE 
C WITH  THE  SIGN  OF  A(OI  (AND  A(K/2)>. 

C 8.  PHI(I)  IS  NOT  COMPUTED  FOR  A(I)»B(I)*0*  FOR  1 .LE. I.LE.NH. 

C 

C— - 

c 

DIMENSION  Y(l)  ,R(D 

DOUBLE  PRECISION  C *Cl  ,C2,S ,S 1, XNP, ANGLE , ALPHA 
EQUIVALENCE  (S,C) 

C 

C 

C ZP'^O  the  OUTPUT  ARRAY. 

C 

NHl=5*NH+7 
00  10  1=1, NHl 
n R(I)  = 0. 

C 

L COHOUTE  At'.  CONSTANTS 

C 

Cl=360./F. 2831653 

XNP=NP 

C2=?./XNF 

ALPHA»f),2831853/XNP 

C 

C COMPUTE  ACC),  3(0)=0,  C(0) , AND  PHI(O). 

C 

S=0. 

DO  20  1=1, NP 
20  S=S*Y(I) 

R(1)=S/XNP 

C 

C COMPUTE  ACI),  B(I),  Cai,  PHICI)  FOR  1.L£*I.LE. K/2-1 

C 

DO  40  T=l,NH 
S=0. 

Sl=0. 

DO  30  J=1,NP 

ANGLE=ALPHA*FL0AT( (J-1)»I) 

S=SfY(J)*OCOS(ANGLE> 

30  S1=S1+Y(J)*0SIN(ANGLE) 

11=5^ I+l 
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R(I1) =S*C? 

R(Il*l) -S1*C2 

R(I1*2)=0S0RT (Cl 
R(Il*4)=R(Il+2) 

C 

C 00  NOT  COMPUTE  PHI(I)  IF  A{I}=B(IlaO. 

C 

if(R(ii).ne.o..op.r(ii*i)«ne.o.)Go  to  32 

WRITE (6,31)1,1 

31  FORMAT (///16H  TO  COHPUT ( PHK  I4,46H),IT  MOULD  BE  NECESSARY  TO  COM 
IPUTE  ATAN2(0,0)./48H  THESE  ARGUMENTS  HILL  NOT  BE  ACCEPTABLE,  SO  PH 
21 ( I4,23H)  HILL  NOT  dE  COMFUTEC.///) 

GO  TO  41 

32  R(Il«-3)  =C1*ATAN?  (R(I1>  ,R(I1+1)I 
40  CONTINUE 

C 

C COMPUTE  A(K/2),  b(K/2)=0,  C(K/2)  , PHKK/2)  IF  H IS  EVEN. 

C 00  NOT  COMPUTF  TWEP  IF  H IS  000  OR  A FULL-POINT  ANALYSIS  IS 

C NOT  COMPUTtn. 

C 

IF(M00(NP,2).NE.C)G0  TO  60 

IF(  (2*NH)  .f.E.  (2»  (NP/2)  ))G0  TO  60 

AI=-1. 

S-1  • 

00  50  1=1, NP 
AI=-AI 

5C  S=S4-AI*Y(I) 

NH1=5*NH+1 
R(NH1)=S/XNP 
R(NHI*2)=ABS(R(NH1) ) 

R(NHlf4)=R(NHl4-2) 

R(NHH-1I  = 0. 

R(NH1*3)=SIGN(90,,P(NH1I ) 

C 

C COMPUTE  CMAX. 

C 

60  C=0. 

NHl=NH«-l 
00  70  1=1, NHl 

70  C=JMAX1(C,08LE  (R(5*l)l  ) 

C 

C COMPUTE  Cdl/CMAX  FOR  O.LE.I.LE.NH. 

C 

DO  80  1=1, NHl 
11=5*1 

80  R(I1)=P(I1)/C 
C 

C COMPUT‘D  THE  CHECKSUMS  IF  A FULL-POINT  ANALYSIS  HAS  COMPUTED. 

C 

IF((2*NH) .NE. (2* (NP/21 IIRETURN 
I1=5*NH1+1 
DO  90  T=1,NH1 
90  R(Il)=R(Il)fR(5*I-4) 

11=114-1 

IF(HOO(NP,2l .EQ.0INHI=NH1-1 
DO  100  1=2, NHl 
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r 


1 

i 


m 100  R(Iil^RtIl)«-R<5*I-3)*0SIN{FL0Am-l**ALPHA» 

RdlJ  »2.*R{I1» 

RETURN 

END 


I 


REFERENCES 


Lerbs,  H.W.,  "Contra-Rotating  Optiraura  Propellers  Operating 
in  a Radially  Non-Uniform  Wake#"  David  Taylor  Model 
Basin  Report  941,  May  1955 

Morgan,  William  B.  and  Wrench,  J.W.,  Jr.,  "Some  Computational 
Aspects  of  Propeller  design,"  Methods  in  Computational 
Physics,  Vol.  4,  Academic  Press  Inc.,  New  York,  p 301- 
331,  1965 

Morgan,  W.B.,  "The  Design  of  Contrarotating  Propellers 
Using  Lerbs'  Theory,"  Transactions  of  the  Society  of 
Naval  Architects  and  Marine  Engineers,  Vol.  68,  p 6- 
38,  1960 

Nelson,  D.M.,  "A  computer  Program  Package  for  Designing 
Wake-adapted  counterrotating  Propellers:  A Users  Manual," 
Naval  Undersea  Center,  Fleet  Engineering  Department 
Report  NUC-TP-494,  December  1975 

Lerbs,  H.W.,  "Moderately  Loaded  Propellers  with  a Finite 
Number  of  Blades  and  an  Arbitrary  Distribution  of  Circulation, 
Transactions  of  the  Society  of  Naval  Architects  and 
Marine  Engineers,  Vol.  60,  p 73-117,  1952 

Hough,  G.R.  and  Ordway,  D.E.,  "The  Generalized  Actuator 
Disk,"  Advanced  Research  Report  TAR-TR-6401,  Therm, 

Inc.,  January  1964 

Caster,  E.B.,  Diskin,  J.A.,  and  LaFone,  T.A.,  "A  Lifting 
Line  Computer  Program  for  the  Preliminary  Design  of 
Propellers,"  David  Taylor  Naval  Ship  Research  and  Development 
Center,  Ship  Performance  Department  Report  SPD-595-01, 

November  1975 

Denny,  Stephen  B. , "Comparisons  of  Experimentally 
Determined  and  Theoretically  Predicted  Pressures  in 
the  Vicinity  of  a Marine  Propeller,"  Naval  Ship  Research 
and  Development  Center  Report  2349,  May  1967 

Morgan,  W.B.,  Silovic,  Vladimir,  and  Denny,  S.B.,  "Propeller 
Lifting  Surface  corrections,"  transactions  of  the  Society 
of  Naval  Architects  and  Marine  Engineers,  Vol.  76, 
p 309-347,  1968 


139 


10.  Bckhart,  M.K.  and  Morgan/  K.B.,  *A  Prqpeller  Design  Method/" 
Transactions  of  the  Society  of  Naval  Architects  and 

Marine  engineers/  Vol.  63/  p 305-370/  1955 

11.  Keller/  J.  Auf'm/  "Snige  aspecten  bij  het  ontwerpen 
van  Scheepsachroeven/*  Schip  en  werf/  No.  24/  p 658- 
662/  1966 

12.  Burrill/  L.C.  and  Baerson/  A./  "Propeller  Cavitation: 

Further  Tests  on  16-Inch  Propeller  Models  in  the  Kings 
college  Cavitation  Tunnel/"  Transactions  of  the  North 
East  Coast  Institution  of  Engineers  and  Ship  Builders/ 

Vol.  78,  p 295-320,  1963-64 

13.  Brockett,  Terry,  "Minimum  Pressure  Envelopes  for  Modified 
NACA  66  Sections  with  NACA  a«0.8  Camber  and  BUSHIPS 

Type  I and  Type  II  Sections,"  David  Taylor  Model 
Basin  Report  1780,  1966 

14.  Cumming,  R.A.,  Dictionary  of  Ship  Hydrodynamics  - 
Propeller  Section,  14th  International  Towing  Tank 
Conference  1975,  Report  of  Presentation  Committee, 

Appendix  VII,  1975 

15.  Tsakonas,  J.  and  Jacobs,  W.R.,  "Counterrotating  and 
Tandem  Propellers  Operating  in  Spacially  Varying, 
Three-Dimensional  Plow  fields,"  Davidson  Laboratory, 

Stevens  Institute  of  Technology,  Report  1335, 

September  1968 

16.  Abbot,  Ira  H.  and  Von  Doehoff,  Albert  E. , "Theory  of 
Wing  Sections  Including  a Summary  of  Airfoil  Data," 

Dover  Publications,  Inc.,  New  York,  Library  of  Congress 
Catalog  Card  No.  60-1601 

17.  Hoerner,  S.P.,  "Fluid-Dynamic  Drag,"  Published  by  the 
author.  Midland  Park,  New  Jersey,  1965 

18.  Lerbs,  H.W.  and  Rader,  H.P.,  "Uber  der  Auftriebsgradienten 
von  Profilen  ira  Propeller  Verband,"  Schiffstechnik, 

Vol.  9,  No.  48,  p 178-180,  1962 

19.  Kerwin,  J.E.,  "The  Solution  of  Propeller  Lifting-Surface 
Problems  by  Vortex  Lattice  Methods,"  Department  of 
Naval  Architecture  and  Marine  Engineering,  Massachusetts 
Institute  of  Technology,  Cambridge,  Massachusetts, 

1961 


140 


20.  Ketwini  J.E.  and  Leopold,  R. , "Propeller  Incidence 
Correction  Due  to  Blade  thickness,"  Journal  of  Ship 
Research,  Vol.  7,  No.  2,  p 1-6,  1963 

21.  Cheng,  H.M.,  "Hydotydna«lc  Aspects  of  Propeller  Design 
Based  On  Lifting-Surface  Theory:  Part  1 - Uniform  Chord 
wise  Load  Distribution,"  David  Taylor  Model  Basin 
Report  1802,  1964 

22.  Cheng,  H.M.,  "Hydrodynamic  Aspects  of  Propeller  Design 
Based  On  Lifting-Surface  Theory:  Part  II  - Arbitrary 
chordwise  Load  Distribution,"  David  Taylor  Model  Basin 
Report  1803,  1965 

23.  Hill,  J.G.,  “The  Design  of  Propellers,"  Transactions 

of  the  Society  of  Naval  Architects  and  Marine  Engineers 
Vol.  57,  1949 


DTNSROC  ISSUES  THREE  TYPES  OF  REPORTS 


111  DTNSROC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER. 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA- 
TION OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION, 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS.  NOT  FOR  GENERAL  DISTRIBUTION. 


